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HAUPTAUFSÄTZE 


Plastisches Knicken der Wandung von Hohlzylindern 
und einige andere Faltungserscheinungen 
an Schalen und Blechen.') 


Von J. W. GECKELER in Jena. 


1. Einleitung. Ein Hohlzylinder, auf dessen Endquerschnitte gleichmäßig verteilte 
Druckkräfte in Richtung der Mantellinien wirken, ist im allgemeinen einer zweifachen Knick- 
gefahr ausgesetzt. Wenn er hinreichend lang ist im Verbältnis zu den Q@uerschnitts- 
abmessungen, und die Last ein gewisses Maß übersteigt, biegt er als Stab in einem 
Bogen aus. Die Knicklast wird in diesem Falle, sofern sie noch in den Bereich der 
vollkommenen Elastizität fällt, durch die Eulersche Formel angegeben. Unter anderen 
Verhältnissen, nämlich bei kürzeren und dünnwandigen Zylindern, ist aber noch ein 
gänzlich anderer Vorgang möglich, der ebenfalls zu den Knickerscheinungen zählt, und 
darin besteht, daß sich bei geradlinig bleibender Achse in der Zylinderwandung ringsum 
laufende Wulste oder Falten bilden, was bei Fortdauer der Belastung ebenfalls zum Zu- 
sammenbruch führt (Abb. ı). Dieser Knickfall wurde 1908 von R. Lorenz?) und von 
S. Timoschenko?) erstmals theoretisch untersucht‘). 

Es lohnt sich aus verschiedenen Gründen dieses Problem nochmals aufzugreifen: 
l. sind in den genannten Arbeiten Randbedingungen vorausgesetzt, die praktisch nicht, 
öder allenfalls nur mit komplizierten Anordnungen zu verwirklichen sind und daher zu 
Ergebnissen führen, die in den gewöhnlichen Fällen mit der Beobachtung nicht über- 
einstimmen, 2. ist es wünschenswert die Untersuchung auf das Gebiet bleibender Form- 
änderungen auszudehnen, weil aus Gründen, die später leicht einzusehen sind, gerade 


!) Vortrag, zehalten gelegentlich der Werkstoffschau, in der Gesellschaft fi. angew. Mathematik 
und Mechanik in Berlin, am 28. Oktober 1927. 

?) R. Lorenz, Zeitschr. des VDI, Bd. 52 (1908), S. 1706. 

3) S, Timoschenko, Mitt. der Techn. Hochschule zu Kiew, 1907/08; Zeitschr. für Math. u. 
Phys., Bd. 58 (1910), S. 337. 

*) Auf das nicht achsensymmetrische Knicken, bei dem die Querschnitte nicht Kreise bleiben, 
sondern zu Vielecken ausknicken, soll hier vorerst nicht eingegangen werden. Eine Untersnchung dieses 
Falles hat ebenfalls R. Lorenz gegeben (Physik. Zeitsehr., 12. Jahrg. (1911), 8. 247. 
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dieses Gebiet von praktischem Interesse ist, 3. kann dieser einfache Fall zur Erläuterung 
einer auch für die Untersuchung verwickelterer Knickvorgänge brauchbaren Methode 
dienen. 


2. Differentialgleichung für das Knicken 
eines Rohres im plastischen Bereich. Für die 
Ableitung der auch jenseits der Elastizitätsgrenze 
gültigen Differentialgleichung ist es zweckmäßig, 
sich aus der zylindrischen Wandung einen Längs- 
streifen von der Breite adg ausgeschnitten zu denken 
und die Gleichung für diesen Stab aufzustellen. Auf 
diesen wirken folgende Lasten: 


l. In axialer Richtung der auf diesen Streifen 
trefiende Anteil an der Drucklast, also T’adg, wenn 
a T der auf die Längeneinheit des Querschnittskreises 
treiiende Betrag ist. 
2. Auf die länder kann ein Moment M und eine (Querkraft () wirken. 


3. Jede Ausweichung w der Erzeugenden — sie soll positiv gerechnet werden, 
wenn sie nach außen geht — ruft tangentiale, also Umfangsspannungen hervor. Deren 
über die Wanddicke genommene Resultante ist auf die Längeneinheit der Mantellinie 

T; = w. 


a 


Diese Kräfte, die man sich anstelle der weggeschnvittenen Zylinderteile längs des 
Stabes angebracht zu denken hat, haben eine in die Richtung der Flächennormale fallende 
Komponente vom Betrag 


T,dy wdgy. 


Auf die Längeneinheit der Stabbreite, also auf die Flächeneinheit des Mantels, 
trifft dann eine Kraft 


a 
Sie stellt einen der Ausbiegung proportionalen Widerstand dar, der das Ausbiegen, 


also das Knicken des Stabes zu hindern sucht. Der Ausdruck e spielt dabei dieselbe 


Rolle, wie in der Theorie des Balkens auf elastischer Unterlage — etwa der Eisenbahn- 
schwelle — die Bettungsziffer. Er soll hier als Widerstandszifier 5 eingeführt werden 
Eö 


Das Fehlen oder Vorhandensein eines solchen elastischen Widerstandes teilt die 
Knickerscheinungen in zwei wesentlich verschiedene Klassen. Zur ersten Klasse gehören 
die elementaren Knickfälle des geraden und gebogenen Stabes, der ebenen Platte und 
das Knicken solcher Schalen, bei denen eine Verbiegung ohne Dehnung der Mittelfläche 
möglich ist. Zur zweiten Klasse gehören die vorausgehenden Fälle dann, wenn ein 
äußerer Widerstand stützend wirkt, also z. B. der Stab mit elastischer Querstützung, 
sowie das Knicken solcher Schalen, deren geometrische Form und Randbedingungen 
keine Verbiegung ohne eine wesentliche Dehnung der Mittelfläche zulassen. Diese 
Dehnung ist in der Wirkung einem äußeren elastischen Widerstand gleichwertig. Während 
bei den Knickfällen der ersten Klasse die elastische Linie oder Fläche im ganzen einen 
Bogen bildet, besteht das Knicken im zweiten Falle in der Bildung meist zahlreicher 
Wellen, Falten oder Beulen. Weil die Schale auf diese Weise in mehrere kleinere 
Bereiche aufgeteilt erscheint, wird dieser Vorgang auch als »lokales Knicken der Wandung« 
bezeichnet '). Alle hier erörterten Fälle gehören dieser zweiten Klasse an. 


Die Diiierentialgleichung des in der oben bezeichneten Weise belasteten Stabes lautet: 


x 
a? 
B“ Tady-w-+| x+Madg=0 (3). 


ı) R. Zoelly, Ueber ein Knickungsproblem an der Kugelschale, Promotionsarbeit, Zürich 1915. 
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Dabei ist B die Biegungssteifigkeit, d.h. das Produkt E/ aus Elastizitätmodul und Träg- 
heitsmoment des Stabquerschnitts. Nach zweimaliger Differentiation und Einsetzen des Wertes 
Bw für g(x) erhält man dv Tau 

wo B sich nun auf die Breite 1 bezieht, so daß 

Eö° Eö 

a? 

Der Faktor — ; (r = Querzahl, — '/,) stellt eine kleine Korrektur dar, die anzubringen 


ist, weil der Querkontraktion wegen der Spannungszustand des Stabes, der sich im Zu- 
sammenhang mit den anschließenden Zylinderteilen befindet, nicht mehr einachsig ist. 
Bis jetzt war vorausgesetzt, daß sich alle Vorgänge im Bereich vollkommener 
Elastizität abspielen. Nur unter dieser Voraussetzung gibt die Lösung der Gl. (4) die 
Vorgänge richtig wieder. Für den geraden Stab haben nun Engesser') und Kärmän’) 
eine über die Elastizitätsgrenze hinaus gültige Verallgemeinerung der Eulerschen Glei- 
chung entwickelt. Diese Methode kann ohne weiteres auch auf den aus dem Rohr aus- 
geschnittenen Elementarstreifen und damit auf dieses selbst angewandt werden. Für die 
Anwendung dieser Theorie ist es notwendig, daß für den betreffenden Stoff das Form- 
änderungsgesetz über die Elastizitätsgrenze hinaus durch Druckversuche an Stäben er- 
mittelt ist. Vorausgesetzt wird sodann: 1. daß sich die Längenänderungen der Fasern 
eines wenig ausgebogenen Stabes aus der Annahme über das Ebenbleiben der Querschnitte 
genügend genau berechnen lassen und 2., daß in den Fasern eines solcben Stabes der 
Zusammenhang zwischen Dehnungen und Spannungen derselbe ist, wie beim einfachen 
Zug-Druckversuch und zwar sowohl im elastischen, wie im plastischen Gebiet. Diese 
Annahmen erweisen sich als gerechtfertigt‘). Abb. 2 zeigt für einen Querschnitt des Elementar- 
streifens die Verteilung der durch die Ausbiegung hinzukommenden Dehnungen und 
Spannungen im Augenblick des Ausknickens. Auf der einen, der konkaven, Seite einer 
sog. neutralen Faser, die aber nicht mehr durch den Schwerpunkt geht, werden die Druck- 


spannungen infolge des Hinzutretens der Biegung noch weiter erhöht. Dort gilt ein Modul 
=, der aus Druckversuchen ermittelt is. Auf der konvexen Seite findet eine 
Entlastung statt, wobei jedoch nur die rein elastischen Förmänderungen rückgängig 
gemacht werden. Hierfür gilt ein Modul Es, der nur wenig von dem gewöhnlichen 
Modul E verschieden ist. Engesser 
und Kärmän haben unter Zugrunde- 
legung einer solchen Spannungsverteilung 
durch einfache Gleichgewichtsbetrach- 
tungen, die ganz denen der elementaren 
Biegungstheorie entsprechen, die Diffe- |, ! 
rentialgleichung der elastischen Linie im | nahe ip 
plastischen Bereich abgeleitet, die sich [KRs«72) 

von der gewöhnlichen für das elastische Abb. 2. 

Gebiet nur dadurch unterscheidet, daß 

an Stelle des Moduls X£ ein Ausdruck K tritt, der die eben eingeführten Größen E, und E3 
in einer von der Form des Querschnitts abhängigen Weise enthält. Für den hier in Be- 
tracht kommenden rechteckigen (Querschnitt ist 


4 Eı 
Dieser sog. Knickmodul ist auf diese Weise eine Funktion der Belastung und nimmt ab, 
wenn letztere wächst. Diese Ueberlegungen bleiben für den aus dem Zylinder ausge- 
schnitten gedachten Längsstreifen gültig. Setzt man also in der Differentialgleichung (4) 
in dem Ausdruck für die Biegungssteifigkeit B anstelle von E diesen Modul X ein, also 


K 03 (6) 


« Dehnungen Spahnungen 


neutr : Faser 


I) F. Engesser, Schweizer. Bauzeitg.. Bd. 26, S. 24 (1895). 

2), Th. v. Kärmän, Dissert. Göttingen, Mitt. üb. Forschungsarb. a- d. Geb. d. Ingen.-Wesens, 
H. 81 (1910); vergl. auch W. Gehler, Verhandl. des 2. Internat. Kongresses für tech. Mechanik. 
Zürich 1926 (Verl. Orell Füssli). 

3) E. Meyer, Zeitsch. des VDI., Bd. 52, S. 167 (1908). 
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so braucht man bei der Auflösung der Gleichung zunächst gar keinen Unterschied zu 
machen, ob es sich um das rein elastische oder das plastische (sebiet handelt. Im ersteren 
Fall geht K in E über. Dagegen darf nicht etwa auch in 8 an die Stelle von EK ein- 
gesetzt werden. Vielmehr kann mit guter Annäherung E unverändert stehen bleiben, 
weil die Umfangsspannungen 7, vor dem Eintritt des Knickens keine hohen Werte an- 
nehmen und selbst eine bereits eingetretene Stauchung in axialer Richtung den Dehnungs- 
modul quer dazu nicht erheblich ändert. Allenfalls müßte dies durch genauere Versuche 
noch mehr geklärt und ein verbesserter Dehnungsmodul ermittelt werden. 


3. Lösung der Differentialgleichung. Diese lautet: 
w— Cı + 09 eh? +6 eh? 
wo ? die Wurzeln der Gleichung 


sind, nämlich 
4 r 
Die reelle Form der Lösung für ww ist verschieden, je nachdem 
m< 


Es wird sich zeigen, daß dieser Ausdruck den Kritischen Wert für die Belastung dar- 
stellt. Setzt man in diesem noch die Werte für B und ß ein, so erhält man dafür 
2% 
— v?) a (1 — v?) Aa 
Das Verhältnis der augenblicklich aufgebrachten Last zur kritischen soll weiterhin mit z 
bezeichnet werden: 


Vor dem Erreichen des kritischen Wertes kann man die Lösung in der Form schreiben 
w— Ch . . (10), 
wo Wı, Wa Integrationskonstanten, und 
4 
.V.? 


Wenn 7T=0, also z= (0 und damit «= y, erhält man die Formeln für die reine Biegung. 
w verläuft dann in einer stark gedämpften Sinuslinie. Aufeinanderfolgende Amplituden 
steherf im Verhältnis e”: 1, also etwa 23,14:1. Die Wellenlänge ist 


2n 4B 

2 - 2a)". . . . (12). 
ß 

Bringt man eine allmählich zunehmende Belastung auf, so nimmt die Wellenlänge dabei 

ab und erreicht bei der kritischen Belastung einen im Verhältnis 1: '/a kleineren Wert, 


nämlich 
4 9 4, 
B K 


Der Dämpfungsexponent hat dabei auf Null abgenommen. Die Lösung ist beim Errei- 
chen der kritischen Belastung in die Form übergegangen: 


Bei weiterer Steigerung der Last, wenn also 2 > 1, lautet die Lösung 


wo Yı und 9» die durch ? dividierten rein imaginären Wurzeln (7) sind. 


Diese Lösungen sind nun den Randbedingungen anzupassen. Dabei ist fol- 
gendes zu beachten: Die Druckbelastung ruft infolge der Querdehnung eine Vergrößerung 
des Zylinderdurchmessers um den Betrag ra; hervor. Im allgemeinen werden bei 
Versuchen in der lestigkeitsmaschine die Druckplatten, bei Konstruktionen die anschließen- 
den Konstruktionsteile die Räder ganz oder teilweise daran hindern, dieser Dehnung zu 


A 
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folgen (Abb. 3). Man hat also von Anfang an mit einer Verbiegung an den Rändern zu 
rechnen und hat für die Ränder bei @=0 und xz=/! nicht wie in der Lorenzschen 


Untersuchung w = 0, sondern w — — Ar zu fordern, da w —= 0 

die unverbogene Mantellinie ist. Auf diesen Umstand hat auch Ir 

L. Föppl!) hingewiesen. Solange man der Knicklast nicht zu =; 
nahe kommt, klingen diese von den Rändern ausgehenden Ver- [ \ 


biegungssysteme rasch ab und die an einem Ende hervor- 
gerufenen Spannungen haben keine merkliche Wirkung auf das 


andere, wenn der Zylinder nicht allzukurz ist (l> etwa 4,8 | aß). 
Bei oder nahe an der kritischen Last ist dies jedoch anders. Es 
zeigte sich, daß in diesem Falle der Dämpfungsexponent ganz 
verschwindet, und wie das Auftreten des Faktors & in der 
Lösung (14) erkennen läßt, spielt nunmehr die Länge des Zy- Abb. 3. 
linders eine gewisse Rolle. Die Ränder wirken gegenseitig ver- 


steifend aufeinander. Es ist nicht mehr möglich, die Randbedingungen an beiden Enden 
getrennt zu erfüllen. 


4. Besprechung der Lösung. Wenn die Rohrenden nicht eingespannt sind, und 
die Robrlänge ein gerades Vielfaches der Wellenlänge A, nach Gl. (13) ist, wird z. B. die 
Form zu der sich die Mantellinie unter der kritischen Belastung ver- 
biegt, dargestellt durch 


Diese Kurve ist für I-8”- und Z—16- in Abb. 4 dargestellt. 


Wenn die Rohrlänge, wie in diesen Beispielen, so groß ist, daß sich 
mehrere Wellen ausbilden können, kann man mit guter Annäherung 


als größte Ausweichung Wax Jene ansehen, die zu = gehört. 


Diese ist dann 


Aehnliche Formen ergeben sich auch für andere Verhältnisse zwischen 
Wellenlänge und Rohrlänge. Wird an der Stelle höchster Bean- 
spruchung (d. i. immer die dem Rande nächste Ausbiegung) die 
Elastizitätsgrenze überschritten, so geht die Formänderung an dieser 
durch die Ausbiegung geschwächten Stelle rasch weiter, während sich 
die übrigen Wellen wieder zurückbilden können. Auf diese Weise Abb. 4. 
erhält man bei Versuchen in der Festigkeits- 

maschine zuerst meist nur eine Welle von der 7” 


| 
typischen, immer nach außen gebogenen Form 24 
(Abb. 5). Ob diese Welle an beiden Enden, »_. HH IHRRLIR 
oder nur an einem, und an welchem Ende sie | 
dann auftritt, hängt von Zufälligkeiten des 
| 
| 
„111 
| 


234567890 N 
ze) 2 


Abb. 5. Abb. 6. 


) L. Föppl, Sitzungsber. d. bayr. Akad. d. Wiss,, math.-phys. Kl., 1926, 
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Materials und der Form ab. Unter Umständen (dünnwandige Rohre, i— (2 n+-,) Re 


schlagartige Belastung) erhält man aber auch eine längs des ganzen Rohres sichtbare 


Wellenbildung. Im allgemeinen bleiben die w für die kritische Last endlich, nur wenn 
die Rohrlänge !=(2n+1) "ist, ergeben sich unendlich große Werte. Stellt man (immer 
Y 


für z—= 1) Wax als Funktion der Rohrlänge dar, so erhält man Abb. 6. 

Zur Kennzeichnung der Verhältnisse bei der Annäherung an die kritische Belastung 
ist es nützlich, die maximale Ausbiegung als Funktion von z darzustellen. Die Grund- 
lage liefert Gl. (10). Unter Zulassung der obigen Vereinfachung für die Berechnung von 


Inn 
Wnax erhält man, wenn != 
(yi-r) 
„”„1-—e 
[77 a 
yı-:»‘ 
7 
und wenn 
a 
2 
U'nex == Ar (18b). 
yı „2 a 
| | | 1—e 
Ir 


| | Die beiden Kurven sind für n=4 in Abb. 7 dar- 
. | | gestellt. Die Kurven, die den dazwischen liegenden 
RItT I Verhältnissen entsprechen, verlaufen in dem Gebiet 
zwischen beiden. Man sieht, daß die Formände- 
rungen erst 5 bis 15 vH unter der kritischen Be- 
lastung stark anwachsen. Praktisch ist auch der 
Unterschied, ob die Rohrlänge ein gerades oder un- 
gerades Vielfaches der Wellenlänge oder ein da- 
zwischen liegender Wert ist, nicht sehr erheblich. 
Die Gl. (1Sb) liefert, wie oben schon gesagt für 
Gl. (15a) dagegen nimmt zunächst 
eine unbestimmte Form an, die aber nach Berech- 


| | | 1 nung des Grenzwertes 
00 085 090 095 100 P 


77) | 


- 


yı 


4 
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die Gl. (17) bestätigt. Zur Berechnung der w über z=1 hinaus 
hat man von Gl. (15) auszugehen. Eigentlich müßte auch immer 
berücksichtigt werden, daß y und damit //A für ein bestimmtes 
Rohr nicht einen festen Wert hat, sondern Funktion von z ist. Die 
Veränderlichkeit ist aber in der Nähe der kritischen Belastung 


nur gering und kann für eine allgemeine Erörterung außer acht 
gelassen werden. 


Abb. S gibt je ein Beispiel für die Verformung der Mantel- 
linien bei beidseitiger und einseitiger fester Einspannung. Im 


ersteren Falle ergibt sich, wenn / ein gerades Vielfaches von g- 


ist, die Cosinuslinie als Lösung, für = (?n+ wie in der 


Abbildung, erbält man: 


Im Falle fester Einspannung bleibt immer w < Jr, solange 
T< T,, darüber hinaus nehmen aber die Formänderungen eben- 
falls rasch zu, 


Abb. 8. 


) 
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Man erkennt, daß es sich unter Beobachtung der hier vorausgesetzten Randbe- 
dingungen im Gegensatz zu der Lorenzschen Lösung nicht eigentlich um ein Stabilitäts- 
problem handelt: eine Verzweigung des Gleichgewichtszustandes kommt an keiner Stelle 
in Frage, weil die triviale Lösung w—=0 von Anfang an der Randverbiegung wegen 
ausgeschlossen ist. Es handelt sich streng genommen um ein reines Spannungsproblem, 
ähnlich wie beim Stab mit zusätzlicher Biegungsbeanspruchung. Immerhin ist das gefähr- 
liche Anwachsen der Formänderunpgen und Spannungen so deutlich auf die nächste Um- 


gebung der kritischen Belastung beschränkt, daß es für die praktische Beurteilung der, 


»Knickgefahr« vollständig genügt, den kritischen Wert 7), zu kennen. Allerdings ist dann 
zu beachten, daß unter Umständen der Zusammenbruch schon bei einer um etwa 10 oder 
15 vH geringeren Last eintreten kann. 


Die kritische Spannung ist gegeben durch 
de 
I+VE 
do 


(wenn . . (20). 


Diese liegt nur dann unter der Elastizitätsgrenze, wenn es sich um außerordentlich 
weite dünnwandige Rohre handelt. Bei Eisen muß z. B. der Halbmesser mindestens 
500mal so groß sein wie die Wanddicke. Dann gehen die Gl. (8) und (20) in die von 
Lorenz und Timoschenko abgeleiteten Werte über. Praktisch werden aber solche 
Rohre kaum verwendet. Bei dick- 
wandigeren Hohlzylindern gelten 
dann die allgemeineren Formeln (8) —- 
und (20), die man etwa in Parallele Kalke \pannung \ | 
stellen kann zu der Tetmever- 5000 \-— 
schen Formel für das unelastische | I\ | | | | | 
Knicken des geraden Stabes. 5000 
Wie erwähnt, muß der Mo- | | | 

dul X ein für allemal für das in Be- #000, 
tracht kommende Material als Funk- 
tion der Spannung ermittelt sein. 3000 
Legt man z. B. die Versuche von Pe | 
Kärmän') mit Siemens-Martin-Stahl- 2000,— 
stäben zugrunde und stellt dann die | | | 

kritische Spannung als Funktion des 700074 ! | — 
Verhältnisses Rohrdurchmesser zu | | | | rein elastisches Knıcken 
Wandstärke dar, so erhält man den Ä 
in Abb. 9 aufgezeichneten Verlauf. 
Bis zu r/ö = 500 gilt die Lorenz- Abb. 9. 
sche Formel, welche sich als Hyperbel 
darstellt. Im plastischen Gebiet bleibt 


| 

| | 

_leßgrenze_| __| 

_ Elastizitätsgrenze 


me) 700 200 300 400 500 600 700 800 9007 1000 


dann die kritische Spannung erheblich 
hinter dieser Hyperbel zurück, um erst 7 rs Ä 
jenseits der Streckgrenze wieder stärker 6 
anzusteigen. Für die Praxis müßte eine 50 78 N 


solche Kurve für jeden Werkstoff ein- #0 
mal aufgezeichnet werden, um dann in 30 9 
jedem Fall ohne weiteres die kritische 20 Ya 
Spannung entnehmen zu können. Stellt 7017 
man die Tragfähigkeit eines Rohres von 44.0) 700 200 300 
bestimmter Wandstärke, z. B. 1 mm als 

Funktion des Halbmessers dar (Abb. 10), Abb. 10. 

so zeigt sich, daß im elastischen Gebiet 

die Knicklast vom Durchmesser ganz unabhängig ist, daß jedoch im plastischen Bereich 
die Knicklast mit dem Durchmesser stark abnimmt und diesem schließlich mehr und mehr 
proportional wird, wie es dem einfachen Stauchen entspricht. Auf diese Weise entsteht 
ein ganz stetiger Uebergang vom Knicken zum Stauchen. 
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I) Anm. 2, S. 343 unten. 
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Alle diese Ergebnisse werden durch Versuche gut bestätigt. 
Auf Erörterungen über den weiteren Verlauf des Vorganges (Bil- 
dung mehrerer Knickwellen nacheinander), über die Abgrenzung 
gegen das Eulersche Knicken im plastischen Bereich usw. soll 
nicht weiter eingegangen werden. Erwähnt sei nur, daß unschwer 
auch Zylinder mit nicht kreisförmigem Querschnitt in entsprechender 
Weise behandelt werden können. Für den elliptischen Zylinder 
ergibt sich z. B., daß die Breite der Knickfalte, deren. Bildung 
an der Stelle geringster Krümmung zuerst beginnt, ringsum pro- 
portional mit der Wurzel aus dem Kıiümmungsradius veränderlich 
ist, was ebenfalls durch den Versuch sehr gut bestätist wird. 
Abb. ı1 zeigt die Knickform für einen Zylinder mit rechteckigem 
Querschnitt. Ferner werden dadurch, daß man nur die allgemeinen 
Begrifte Biegungssteifigkeit, Dehnungssteifigkeit (#d), Widerstands- 
ziifer eingefiihrt, auch Zylinder mit komplizierterer Wandkonstruk- 
Abb. 11. tion (Netzwerkkonstruktionen, Zylinder aus Wellblech, solche mit 
Versteifungsrippen u. dergl.) der Untersuchung zugänglich. 


5. Faltenbildung an Behältern und Kesseln unter innerem Ueberdruck. 
Wie schon einleited gesagt, läßt sich auf diese einfache Weise nicht nur dieses Knick- 
probem des Zylinders, sondern noch eine größere Zahl anderer Aufgaben über das Knicken 
von gewölbten Schalen und ebenen Blecben lösen. Man schneidet aus dem auf Stabilität 
zu untersuchenden Objekt in geeigneter Weise einen 
Streifen aus, und ersetzt den Widerstand, der der Ver- 
biegung, solarge der Zusammenhang noch gewahrt ist, 
aus geometrischen Gründen entgegenwirkt, durch äußere 
elastische Kräfte, die sich dann wie eine elastische 
Bettung auswirken. 

Es möge z B. an ein abgeplattetes Rotations- 
ellipsoid (oder überhaupt an einen ähnlich geformten 
einfach zusammenhbängenden Behälter) gedacht werden 
(Abb. 12). Ein innerer Ueberdruck ruft in einem solchen 
längs eines Streifens von bestimmter Breite zu beiden 
Seiten des Aequators Druckspannungen 7% hervor, wenn 


Abb. 12. das Ellipsoid flacher ist, als dem Achsenverhältnis Ya:ı 

entspricht (oder allgemein, wenn der Meridiankrümmungs- 

radius kleiner ist als die Hälfte des 2. Krümmungshalbmessers Z,. d. h. der Normale bis 
zur Achse). Der im Aequator erreichte Höchstwert dieser Druckspannungen ist 


ap a 


wo e der Meridiankrüimmungsradius am Aequator ist. 

Ein sehr schmaler Streifen zu beiden Seiten des Aequators befindet sich nun in 
derselben Lage, wie ein auf Druck beanspruchter Kreisring, dessen Ausknicken in seiner 
Ebene durch einen radialen elastischen Widerstand behindert wird. Dieser Widerstand 
rührt von den mit dem Ausknicken notwendig verbundenen Dehnungen der Schalen- 
mittelfläche her. Es ist leicht zu zeigen, daß auch hier die Widerstandsziiier 


Eö 


ist. Erreichen die Druckspannungen in einem solchen Ring einen bestimmten kritischen 
Wert, so kniekt er, indem er sich längs seines ganzen Umfanges in Wellen legt. Das- 
selbe gilt fiir die Aequatorzone des Ellipsoids (Abb. 13). Die Rechnung entspricht ganz 
der für den geraden Sıab durchgeführten, nur hat man von der für den krummen Stab 
gültigen Gleichung der elastischen Linie 


B w ) 
wi == M 
ag? + 


auszugehen. Man findet für die kritische Spannung 


(22 


a: 


| 
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In diesem Falle handelt es sich übrigens um ein wirkliches Stabilitätsproblem. 
Bei Schalen, die so dünnwandig sind, daß ein Knicken in Frage kommt, ist 
1 
38T 67a 


ein so kleiner Wert, daß man einen nur sehr geringen Fehler macht, wenn man ihn in 
Gl. (23) gegenüber 1 vernachlässigt und wieder schreibt 


Diese Vereinfachung bedeutet, daß man für ein Stick des Aeqnatorkreises von der 
Länge einer Welle die Krümmung vernachlässigt. Die Wellenlänge ergibt sich wieder zu 


4 4, 
/ß 


Setzt man schließlich noch für 
B und P die Werte eın, und benutzt 
man noch die Gl. (21), so erhält man 
für den kritischen Ueberdruck 


2 1 ES? 


Diese Formel wurde durch eine 
erößere Reihe von Versuchen mit 
dünnwandigen, aus Messingblech ge- 
drückten Ellıpsoiden gut bestätigt. 
Allerdings findet, ähnlich wie auch bei 
anderen Knickversuchen, das Knicken 
infolge schwer zu vermeidender Fehler 
in Form, Wanddicke und Werkstoff- 
beschaffes,heit meist schon bei einem 
etwas geringeren Druck statt. Häufig 
tritt dann auch die Wellenbildung 
nicht längs des ganzen Umfanges in 
Erscheinung, sondern es bildet sich 
nur eine kräftige Falte an der 
schwächsten Stelle. Die Wellenlänge 
stimmt mit der aus Gl. (24) gerech- 
neten in allen Fällen sehr gut überein. 

Bemerkenswert ist, daß grund- 
sätzlich auch die Krempe eines Kessel- 
bodens in dieser Weise knicken kann 
(Abb. 14). Mittels der oben abge- 
leiteten Gleichungen und der früher 
für diesen Fall entwickelten Theorie ') 
ist es nicht schwer, den kritischen 
Druck zu berechnen. Bei den ge- 
bräuchlichen Kesselformen ist eine 
Knickgefahr im allgemeinen nicht vor- 
handen, wenn man auch bei einem 
zu geringen Bördelradius dem kriti- 
schen Wert nahekommen kann. Das 
Kesselmodell, Abb. 14, müßte mit 
a=150 mm, ö=0,28 mm nach Gt!. (24') 
n= 42 Wellen zeigen, während der Versuch 44 ergibt. 


Mit a=eg— R geht GI. (25) in die von Zoelly?) und Schwerin‘) für das Knicken 
der Kugelschale abgeleitete über. 


Pr = 


(25). 


Abb. 14. 


I) Forschungsarbeiten auf dem Gebiet des Ingenieurwesens, Heft 276. 
?2) Anm. 1, S 342 unten. 


3) E. Schwerin, Zeitsch. f. angew. Math. u. Mech., Bd. 2, S. 81, 1922. 
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6. Faltenbildungen an ebenen Blechen, insbesondere beim Ziehprozeß, 
Beachtenswert ist schließlich noch ein anderer Grenzfall: das Doppelblatt und das ebene 
Blech. Das Doppelblatt kann man sich aus dem abgeplatteten Drehellipsoid durch all- 
mähliche Verkleinerung der Drehachse bis auf Null entstanden denken. Wenn man ein 
solches, etwa aus zwei kreisrunden, längs des Randes aneinander geklebten Blättern 
hergestelltes Doppelblatt aufbläst, bilden sich in einem äußeren Riog radial laufende 


a 
(Abb. 15). Verwandt hiermit sind die Faltungserscheinungen, wie sie oft bei Arbeiten 
an der Ziehpresse zu beobachten sind. Wird z. B. ein ebenes Blech mittels des Zieh- 
stempels in eine tellerförmige Matrize 
gedrückt (oder umgekehrt, wie bei der 
folgenden Abb. 18), und auf diese 
Weise radial nach der Mitte zu kon- 
trahiert, so entstehen in der über- 
stehenden Randzone tangentiale Druck- 
spannungen, die ein wellenförmiges 
Ausknicken verursachen (Abb. 16). 
Zur Untersuchung dieser Er- 
scheinungen geht man zweckmäßig 


Falten, die so weit reichen, wie die Druckspannungen, nämlich von r=abisr = 


[RAses 215] 
Abb. 15. Abb, 16. 


von der Betrachtung der Stabilität eines radial gedrückten Kreisringes gegen räumliche 
Verzerrung aus. Dem Umstand, daß der Ring nicht frei beweglich, sondern längs seines 
Innenrandes festgehalten ist, kann wieder durch Annahme einer elastischen Querstützung 
(normal zur Ringebene) Rechnung getragen werden. Jeder Verschiebung w eines Punktes 
normal zur Ringebene wirkt eine Kraft g= Pw pro Flächeneinheit entgegen, die unten 
näher untersucht wird. Der auf die Länge 1 der Mittellinie des betrachteten Ringes 
(r = a) bezogene radiale Druck sei p (kg/cm). Labilität tritt ein, wenn 


B 
n? + 


_ 

n stellt die Anzahl der Vollwellen auf den Umfang dar. 5 ist die Breite des 

Ringes, B die Biegungssteifigkeit pro Breite 1 gegenüber einer Verbiegung des Ringes 

aus seiner Ebene heraus, also B= — C ist die Torsionssteifigkeit. Sie wird für flache 


Rechtecke (Ringquerschnitt), wenn b > 36 ist, genügend genau durch die Näherungsformel 


15 
angegeben. Das in Gl. (26) vorkommende Verhältnis B/C bleibt demnach meist kleiner 
als 1, und sogar für das Quadrat kleiner als 1,5, so daß man es mindestens für den 
Zweck einer allgemeinen Erörterung gegen n? vernachlässigen kann. Die Torsion spielt 
keine merkliche Rolle, wenn es sich um eine einigermaßen große Wellenzahl handelt. 
Setzt man noch n? -- I mit n? gleich, so wird 


a’ n? 


Der Wert n, der zur kleinsten kritischen Belastung (dem kleinsten Eigenwert) 


führt, wird aus der Bedingung sr — 0 gefunden zu 


ww. 
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Dies ist die Wellenzahl auf den ganzen Umfang. Der zugehörige kritische Radial- 
druck ist 


b 
Eine Wellung des Ringes tritt also ein, wenn die in seinem (uerschnitt über- 
tragene Druckresultante (66) T=2VBß . (29) 


ist. 

Es ist nun zu überlegen, welchen Wert man für die Widerstandsziiier 5 einzusetzen hat. 

Bei einer Ausbiegung des Ringes aus seiner Ebene heraus hat man zunächst an 
den Biegungswiderstand als widerstebende Kraft zu denken, und zwar ist jener Biegungs- 
widerstand gemeint, den der längs seines Innenrandes zwischen Stempel und Matritze 
fest eingespannte Ring einer achsensymmetrischen Aufbiegung entgegensetzt. Diesen 
Widerstand kann man sich von gleichmäßig über die Ringoberfläche verteilten Kräften y 
geleistet denken, die normal zur Ringfläche wirken, der Ausbiegung proportional sind 
(»Bettungsdruck«), und für sich allein die entgegengesetzt gleiche Durchbiegung der ring- 
förmigen Platte ergeben würden. Dem entspricht eine in radialer Richtung veränderliche 
Widerstandsziffer, die am Außenrand den kleinsten Wert hat; dort ist die Knickgefahr 
am größten. 

Der Zusammenhang zwischen Belastung 7 und größter Ausbiegung w am Rande 


kann in der Form IT 
8 B 


dargestellt werden, wo # ein vom Verhältnis 7 abhängiger, aber von I nur wenig ver- 


schiedener Faktor ist. Nimmt z.B. der Rand '/; des Gesamthalbmessers ein, so ist 
“= 1,102, ist der Rand nur '/‚o des Halbmessers, so ist = 1,036. Mit #— 1 wäre die 


Formel identisch mit jener für einen Kragbalken von der Länge b. Ausg= folgt, 
daß der Ausdruck — die Rolle der Widerstandszahl am Außenrande spielt, also 
sB 
(30). 
Damit gehen die Gl. (27) und (29) über in 
a 8 . 4B 2 
Setzt man für x einen mittleren Wert, etwa 1,07 ein, so wird 


z. B. erhält man mit a=54mm, b= 12,5 mm 
nach dieser Formel eine Wellenzahl n—= 17,1, 
was durch den Versuch (Abb. 17) bestätigt wird. 
Praktisch kommen natürlich nur ganze Zahlen in 
Betracht Es ist bemerkenswert, daß bei diesem 
Vorgang die Wellerzahl von der Blechdicke ganz 
unabhängig ist. 

Bei der Fabrikation müssen natürlich Abb. 17. 
diese Faltenbildungen vermieden werden. Man 
erreicht dies, indem man einen sogenannten Niederhalter mit pneumatischer oder 
mit Federkraft gegen den ausknickenden Rand preßt. Abb. 18 ist eine schematische 
Darstellung der Ziehpresse, wobei, entsprechend dem oben behandelten Vorgang, der 
Niederhalter noch nicht in Wirksamkeit getreten ist. Dagegen ist dies in Abb. 19 bereits 
der Fall. Dann wird natürlich die Widerstandszifier über den Betrag der Gl. (30) hinaus 
erhöht, und damit die Wellenzahl größer, als die Gl. (31) angibt. Der weitere Verlauf 
des Vorganges hängt dann von Einzelheiten der Zieheinrichtung und der Ziehform ab. 
Häufig spielen bei stärkerer Durchbiegung, oder nachdem der Niederhalter die weitere 
Vergrößerung der Biegungswellen hemmt, oder wenn die Form zu starken kegelförmigen 
Aufbiegungen des Randes zwingt, die tangentialen Dehnungsspannungen die Hauptrolle. 


Im letzteren Falle wird Eö 
= . . D . . . (32). 
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Dem entspricht 


( E 


Hier ist dann auch die Blechdicke wieder von Einfluß auf die Zahl der Wellen. Letztere 


ist größer, als die von Gl. (31) angegebene, wenn 5> 1,06 |ad. Für das in Abb. 16 
gezeigte Beispiel liefert die Formel (33) mit «= 150 mm (mittlerer Randhalbmesser) und 
ö=1imm die Wellenzahl n = 19,1, während tatsächlich 20 Wellen — in der Nähe 


Abb. 18. Abb. 19. 
des inneren Randes — zu zählen sind. Die gegen den Außenrand hin beobachtete 


Gabelung einzelner Wellen in zwei ist hauptsächlich auf eine durch den Niederhalter 
verursachte spätere Knickung des flach gedrückten Wellenscheitels zurückzuführen. 
Genauer soll auf diese Einzelheiten hier nicht eingegangen werden. 

In diesem und dem vorausgehenden Abschnitt ist auf die Verhältnisse nach dem 
Ueberschreiten der Elastizitätsgrenze nicht mehr besonders Rücksicht genommen. Dies 
kann selbstverständlich in derselben Weise wie beim Zylinder geschehen, indem man in 
3 oder $ oder in beiden je nach den besonderen Umständen den Modul X oder E, an- 
stelle von E einführt. s64 


Über die Stabilität der Kärmänschen Wirbelstraße 


gegenüber beliebigen Störungen in drei Dimensionen.’ 
Von KARL SCHLAYER in Berlin. 


ie Frage nach einer stabilen Anordnung von Wirbeln in einer doppelten Reihe 
(»Wirbelstraße«) ist erstmalig von Kärmän?) im Jahre 1910 behandelt worden 
zum Zweck einer Erklärung des quadratischen Widerstandsgesetzes in Flüssig- 
keiten und Gasen. Der Geltungsbereich dieses quadratischen Gesetzes (im Gegensatz 
zum linearen Stokesschen) ist bekanntlich der, wo man die Flüssigkeit als reibungslos 
ansehen kann, was um so genauer zutrifft, je kleiner die Zähigkeitskonstante bzw. je 
größer die Geschwindigkeit des durchgeschleppten Körpers ist. Da weder die sogenannte 
»stetige« Potentialströmung noch die durch Annahme eines »toten« Wassers hinter dem 
Körper entstehende »unstetige« Potentialströmung diesen Grenzfall richtig darzustellen 


I) Dissertation, Universität München; Referenten: Prof. Sommerfeld und Prof. W. Wien. 
Januar 1927. 

) Kärmän, Göttinger Nachrichten 1911, S. 509 ff. und 1912, S. 547 ff., sowie Kärmän und 
Rubach, Phys. Zeitschr., 13. Jahrg. (1912), S. 49 ff. Vergl. auch Lamb, Hydrodynamics, 4. Auflage, 
Kap. VII, Art. 156. 
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vermögen, da ferner die primitivste Beobachtung bereits deutlich die Tendenz der Flüssig- 
keit im »Kielwasser« erkennen läßt, einzelne diskrete Wirbelfäden zu bilden, lag es nahe, 
nach möglichen stabilen Wirbel-Anordnungen zu suchen, die eventuell hinter dem Körper 
entstehen könnten (»als Endprodukt der aufgelösten Wirbelschicht«)'). Gibt es eine der- 
artige stabile Wirbelanordnung, die vor den übrigen ausgezeichnet ist, so gestattet sie, 
in eindeutiger Weise die im Mittel pro Zeiteinheit beim Schleppen des Körpers durch 
die Flüssigkeit erzeugte Impulsmenrge und damit den zu überwindenden Widerstand zu 
berechnen. Und es zeigt sich, daß der so gewonnene Ausdruck für den Widerstand in 
der Tat auch quantitativ gut mit den Beobachtungen übereinstimmt. 

Die Existenz einer derartigen ausgezeichneten Wirbelkonstellation zu erweisen, 
bildet den (regenstand der angeführen Kärmänschen Arbeiten. Die Wirbel werden 
dort als streng linear angenommen, so daß eine Verschiebung von Wirbeln längs ihrer 
eigenen Achse an der Konstellation nichts ändert und das Problem als zweidimensional 
behandelt werden kann. Auf Grund dieser Idealisierung läßt sich die Aufgabe in ge- 
schlossener Form bewältigen, bekanntlich mit dem sehr bemerkenswerten Ergebnis, daß 
die »symmetrische« Anordnung (wo je zwei Wirbel sich direkt gegenüberstehen) instabil, 
aber auch die »unsymmetrische« (s. Abb. 2, S. 356) nur dann stabil ist, wenn der Quer- 
abstand zum Längsabstand in einem ganz bestimmten zahlenmäßigen Verhältnis steht, 
nämlich 0,281:1. Dies Resultat kommt dadurch zustande, daß bei gegebener »Form« 
der Straße (nämlich: »alternierende Doppelreihe«) durch Längs- und Querabstand der 
Wirbel die Einwirkung die ein Wirbel erleidet, nach Größe und Richtung eindeutig be- 
stimmt ist, und daß es ein und nur ein bestimmtes Verhältnis der beiden genannten 
Abstände gibt, für welches ein etwas verschobenes (»gestörtes«) Wirbelelement eine Ge- 
schwindigkeit genau in Richtung auf seinen »Ruhepunkt« zu bekommt. 

Es ist nun interessant zu untersuchen, ob dieses Ergebnis der zweidimensionalen 
Behandlung auch dann noch bestehen bleibt, wenn man die Wirbel nicht mehr als streng 
linear auffaßt, sondern auch verbiegende Störungen zuläßt, wie sie dadurch zustande 
kommen, daß die störenden Impulse nicht jeden Wirbel in seiner ganzen Länge mit 
gleicher Größe und Richtung erfassen, sondern das einzelne Element für sich in beliebiger 
Weise. Dadurch wird das Problem offenbar zu einem dreidimensionalen, und es fragt 
sich, ob auch noch hier, angesichts solcher allgemeinster Störungen, die von Kärmän 
aufgestellte Bedingung ausreicht, um die Stabilität des Wirbelsystems sicherzustellen. 

Was die Behandlung des so verallgemeinerten Problems von der des zZweidimen- 
sionalen wesentlich unterscheidet, ist die Tatsache, daß bei einem verbogenen Wirbel 
nicht mehr der ganze Wirbel durch ein einzelnes Element symbolisiert werden kann, 
sondern daß jedes Element eines jeden Wirbels für sich hinsichtlich seiner Bewegungen 
betrachtet werden muß. Dabei wirkt jedes Element des gesamten Systems auf jedes 
andere ein, insbesondere auch innerhalb eines und desselben Wirbels, welch letztere 
Wirkung als »Selbstinduktion« bezeichnet sei. Bei der Selbstinduktion ergibt sich ins- 
berondere die Schwierigkeit, daß unmittelbar benachbarte Elemente unendlich große Ge- 
schwindigkeiten aneinander hervorrufen, solange man den »Wirbelquerschnitt«, wie beim 
zweidimensionalen Problem, als verschwindend annimmt. Diese Schwierigkeit gilt es zu- 
nächst zu überwinden, d. h. es ist festzustellen, welche Gestaltänderungen ein einzelner 
unendlich schwach verbogener Wirbel an sich selbst induziert, und ob er ein stabiles 
Gebilde ist, oder etwa im Laufe der Zeit sich auflöst und aufhört ein zusammenhängender 
Wirbel zu sein. Erst wenn man erkannt hat, daß auch noch angesichts dreidimensionaler 
Störungen sinnvoll von einem »Wirbel« gesprochen werden darf, kann man den Einfluß 
mehrerer derartiger Wirbel aufeinander untersuchen. 


1. Einzelwirbel mit Selbstinduktion. Es sei also betrachtet ein einzelner — 
zunächst unendlich langer — Wirbel, dessen verschiedene Elemente von der Ruhelage, 
in welcher der Wirbel streng linear wäre, um unendlich kleine Strecken nach beliebigen Rich- 
tungen entfernt sind. In allgemeiner Weise ist die Frage nach der Gestaltänderung eines solchen 
Wirbels mit der Zeit bereits von Lord Kelvin behandelt worden. (»Mathematical and 
Physical Papers«, Vol. IV, p. 152 fi. sowie 202 ff.)’),. Von dessen Ergebnissen gilt es, 
diejenigen auszusondern und in eine für uns passende Form umzuschreiben, welche für 
den Ansatz unseres speziellen Problems benötigt werden, d.h. diejenigen, welche einer 


!) Sorgfältige Beobachtungen über derartige Wirbelbildungen waren früher von H. Beward 
veröffentlicht worden (Comptes rendus, Paris 1908). 


2) S.a. Hicks, Transact. of the Royal Soeiety, Vol. 175 (1884), p. 161 ft. 
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Verbiegung des Wirbels als ganzen entsprechen, nicht aber einer bloßen gegenseitigen 
Verlagerung der verschiedenen Wirbelemente, aus welchen sich der Wirbelkern zusammen- 
setzt. Eine derartige Verlagerung wie die letztgenannte kommt für die Wirkung des 
Wirbels nach außen nicht in Betracht, ist also für die Wechselwirkung verschiedener 
Wirbel aufeinander und somit für die Stabilität der Straße bedeutungslos. Die Methode 
durch welche Kelvin zu seinen Resultaten gelangt, ist folgende: es wird betrachtet 
der Fall einer inkompressiblen llüssigkeit, welche sich in einer annähernd kreisförmigen 
Bewegung um eine gerade Achse (die »Wirbelachse«) befindet. Zu diesem Zwecke werden 
die Eulerschen Bewegungsgleichungen dieser Flüssigkeit in ihrer auf Zylinderkoordinaten 
©, z transformierten Gestalt angesetzt. 

Die Forderung annähernd kreisförmiger Flüssigkeitsbewegung um die z Achse läßt 
sich erfüllen, wenn man für diese Gleichungen ein Lösungssystem ansetzt in der Form 
(p bezeichne den Druck): 


r—=ocosmz-sin(nt —i%), — T+rceosme:cos(nt — i®) 
z= wsinmz-sin(nt — id), p=P+mwcosmz:cos(nt — 
mit "T?dr 
r 


wobei eingeführt sind die »Störungsamplituden« E, 7, w und ®, welche offenbar Fank- 
tionen von r sind und voraussetzungsgemäß sehr klein gegen T, die ungestörte Bahn- 
geschwindigkeit eines Flüssigkeitsteilchens. Es sind ferner eingeführt die drei Parameter 
m, n und i, welche eine einfache geometrische Bedeutung besitzen: 


m bedeutet (bis auf den konstanten Faktor 2) die reziproke Wellenlänge der 
Störung längs der z-Achse; 

n bedeutet die Kreisfrequenz der Störung bei ihrer Öszillation ||z und || vr; 

’ ist eine beliebige ganze Zahl zwischen 0 und », und bedeutet die Anzahl der 
»Sinuswellen«, die man als ursprüngliche Deformation (zur Zeit /—= 0) jeder 
einzelnen der Kreisbahnen aufzwingen will, welche in ihrer Gesamtheit unseren 
Wirbel darstellen. Zeichnet man eine dieser Kreisbahnen an, so lassen sich 
die in der z-Ebene gelegenen Komponenten der beispielsweise den beiden 
Werten ’— 1 und ?= 6 entsprechenden Deformation durch die folgende Ab- 
bildung veranschaulichen. 


Dem Falle ©— 0 entspricht ersichtlich 
eine undeformierte Kreisbahn, so daß hier 
die Störungsbewegung keinen Umlauf von 
»Ausbuchtungen« um die Wirbelachse dar- 
stellt, sondern lediglich eine Pulsation der 
Kreisbahnen als Ganzes, d. h. eine suk- 
zessive infinitesimale Vergrößerung und 
Verkleinerung des Bahnradius, verbunden 
mit einer synchronen ÖOszillation aller einem 

Abb. 1. 2 = Schnitt zugehörigen Teilchen parallel 
der z= Achse. 

In dem allgemeinen Lösungsansatz (I) ist die ungestörte Bahngeschwindigkeit 7 
der Flüssigkeitsteilchen zunächst in willkürlicher Weise von r abhängig. Die Speziali- 
sierung auf einen eigentlichen Wirbel, d. h. ein Gebilde mit einem als fester Körper 
rotierenden Wirbelkern und einer sich nach außen hin (stetig!) anschließenden rotationslosen 
Flüssigkeitsbewegung (Geschwindigkeit abnehmend nach dem 1/r-Gesetz) gestattet es 
vermöge der so gewonnenen Randbedingungen, ganz bestimmte Werte für die Störungs- 
frequenz (bzw. Kreisfrequenz) n zu errechnen (s. Kelvin, ]l. c. p. 161 u. 162 sowie 155/56). 
Diese ist naturgemäß verschieden, je nachdem welche infinitesimale Gestaltänderung man 
dem Wirbel anfänglich erteilt, d.h. also je nach der Größe der ganzen Zahl ;. 

Die Fälle, die Kelvin näher untersucht, sind gegeben durch ?—= 0 undi=1, und 
dies sind auch — wie später gezeigt werden wird «(s. im 3. Abschnitt) — die einzigen, 
die für uns Bedeutung besitzen. Man erhält übersichtliche Resultate unter der Voraus- 
setzung, daß m-a << 1 sein soll, d. h. physikalisch gesprochen, daß man nur solche 
Störungsformen betrachtet, bei welchen die Wellenlänge der Störung längs 


z sehr groß ist verglichen mit dem Radius des Wirbelkerns. Machen wir diese 
Einschränkung, so gilt: 


i=7 


A 
3 
- 
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1. für 
ma 


q 
wo qg die Wurzeln der Gleichung Ju (q) = 0 bedeutet, wenn Joa die Besselsche 
Funktion erster Art nullter Ordnung. Da deren Anzahl unendlich groß ist, so 
ist n unendlich vieldeutig. Der kleinste Wurzelwert ist angenähert q — 2,4; 
2. für ©=1 als kleinster Wert unter den auch in diesem Falle an sich unend- 
lich vielen: 


"mia! (In 0,366) 

2 m a 
d.h. der Umlaufsinn der Störung ist in diesem Fall dem des Wirbelkerns ent- 
gegengesetzt, da es sich ja hier, im Gegensatz zum Fall ?=o, tatsächlich um 


eine die Wirbelachse umkreisende Störung handelt. 


Anm.: Daß wir im Falle =] in der Tat nur den kleinsten der n-Werte zu 
betrachten brauchen, liegt daran, daß hier dieser Wert der Größenordnung nach sehr 
klein ist gegen alle übrigen Werte, und daß darum bei der Aufspaltung in unendlich 
viele »Teilwellen« — entsprechend den » vielen n-Werten — (welche erfolgt, wenn wir 
den Wirbel in die durch Abb. 1 angedeutete Gestalt deformieren und dann sich selbst 
überlassen), derjenigen Welle, welche die Achse mit der Winkelgeschwindigkeit n-i. um- 
kreist, eine sehr viel größere Amplitude zukommt, als allen ürigen Teilwellen'). Durch 
die entsprechende Ueberlegung folgt dagegen für ’— 0, daßin diesem Falle die Amplituden 
der den »ersten« (d. h. hier größten) n-Werten entsprechenden Teilschwingungen unter 
sich der Größenordnung nach gleich sind und zwar abnehmend mit abnehmendem n. 

An Hand der von Kelvin ].c. gegebenen Formeln läßt sich leicht die kinematische 
Bedeutung der beiden letztgenannten Störungstypen ermitteln. Es zeigt sich, daß (unter 
der Voraussetzung m-a X 1) abgesehen von kleinen Größen zweiter und höherer 
Ordnung im Falle i= 0 eine reine ÖOszillation der Flüssigkeitsteilchen || z vorliegt, im 
Falle ©—= 1 dagegen eine reine Rotation des verbogenen Wirbels um seine ungestörte 
Ache herum; und zwar beides so, daß sowohl die Dicke als auch der kreisförmige Quer- 
schnitt des Wirbels dabei erhalten bleibt. Physikalisch bedeutsam ist dieses Ergebnis 
darum, weil daraus hervorgeht, daß ein anfänglich durch Stauchung oder Dehnung ( || z) 
gestörter Wirbel dauernd nur ÖOszillationen ||z ausführt (also Bewegungen ohne x- und 
y-Komponente), und ebenso ein ohne z2-Verschiebung verbogener Wirbel dauernd nur 
um seine ursprüngliche Achse rotiert ohne in der z-Richtung zu schwingen. Somit sind 
Longitudinal- und Transversal-Störungen des Wirbels voneinander unabhängig, und 
dieser Umstand ist es, der weiter unten (Abschn. 3) die allgemeinste Störung der Wirbel- 
straße durch die beiden eben besprochenen Störungstypen Kelvins darzustellen erlaubt. 

Ebenso wie für ©=0 und i= 1 läßt sich übrigens auch für alle anderen :-Werte 
zeigen, daß imaginäre n-Werte niemals auftreten können, wodurch sich also der gestörte 
Einzelwirbel ganz allgemein als ein stabiles Gebilde erweist. 

Zum Schluß möge noch bemerkt werden, daß aus den obigen »harmonischen « 
Lösungen (Ansatz I) sich in Fourierscher Weise ganz beliebige Störungsformen aufbauen 
lassen, so daß die Kelvinschen Ansätze alle überhaupt möglichen infinitesimalen Störungen 
umfassen. 


2. Die allgemein gestörte Wirbelstraße. Einführung der Störungen 
Nachdem so die Stabilität des dreidimensional gestörten Einzelwirbels gesichert und die 
Geschwindigkeit bekannt ist, die seine einzelnen Elemente infolge der Selbstinduktion 
erhalten, können wir an den Kernpunkt unseres Problems herantreten, die Wirkung der- 
artiger allgemeiner Störungen auf eine ganze Straße von linearen Wirbeln zu untersuchen. 

Wir betrachten ein bestimmtes, dem gten Wirbel angehöriges Wirbelelement (» Auf- 
punkt« oder »Aufelement«) unter der Einwirkung aller anderen Wirbel sowie aller übrigen 
Elemente desselben Wirbels. Dabei sei die Straße als unendlich lang vorausgesetzt, 
der Querschnittsradius jedes Wirbels sehr klein verglichen mit dem gegenseitigen Ab- 


a. zweier Wirbel der Straße. Ferner seien folgende Bezeichnungen eingeführt 
8. Abb. 2): 


!) Siehe Kelvin Il. e., Bemerkungen auf S. 158 untem u. 159 oben, sowie 161 oben. Diese Be- 


merkungen beziehen sich zwar zunächst auf den Fall des Hohlwirbels, wo nur 2 mögliche n-Werte 
existieren, lassen sich aber ohne weiteres auf unseren Fall übertragen. x 
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h für den Abstand der beiden Wirbelreihen im ungestörten Zustand, 

für den Abstand zweier aufeinanderfolgender Wirbel derselben Reihe, 

d für den Abstand der beiden die Flüssigkeit begrenzenden Platten (Ebenen 
z— 5 d/2) voneinander, | 

p (und g) für die Nummer eines betrachteten Wirbels (von — » bis + ©), 

“ für die Wirbelstärke oder »Zirkulation«, definiert durch die Gleichung: x = wo 
= Winkelgeschwindigkeit X Fläche des Wirbelkerns. 


Die x-Achse liege parallel zu den beiden Wirbelreihen und zwar so, daß + x 
in der Richtung wachsender positiver Wirbelnummern wächst. Die y-Achse sei so ge- 
richtet, daß der »ersten« Reihe negative, der »zweiten« Reihe positive, dem Betrag nach 
gleiche y-Werte zukommen, wie dies auch aus der folgenden Abbildung ersichtlich. 

Die Wirbelstärke # sei in der 
kz ersten Reihe positiv, in der zwei- 


| ten Reihe negativ angenommen. 
oO) | O) 
| Indem wir nun die ganz allge- 
meine störende Verrückung in drei 


x Dimensionen, die wir jedem Wirbel- 
element eines jeden Wirbels der un- 
y4 endlich langen Straße erteilen wollen, 
wirklich einführen, gelangen wir zu- 
nächst zu einer unendlicheu Anzahl 
von Gleichungen mit unendlich vielen 
d Unbekannten. Jedoch gelingt es durch 
|? einen Ansatz, der eine direkte Verall- 


m gemeinerung des zum gleichen Zwecke 
A von Kärmän beim zweidimensionalen 
7 Reihe Problem gemachten darstellt, die Ver- 


m 7 rückung eines jeden beliebigen der 
h unendlich vielen Elemente durch die- 
jenige einer endlichen Zahl von Ele- 
menten auszudrücken. Da es vorteil- 
haft ist, die Stabilitätsbedingung (wie ja auch Kärmän es tut) in komplexer Form auf- 
zustellen, so müssen die Störangen selbst von vornherein komplex definiert werden. 
Dies geschieht durch den Ansatz: 


Abb. 2. 


+ ii - de 


wobei der Index 1 sich auf die »erste«, der Index 2 auf die »zweite« Wirbelreihe 
bezieht und die Nummerierung in der Straße so vorgenommen gedacht ist, daß die 
geraden (ganzen) Zahlen p, den Wirbeln der ersten, die ungeraden Zahlen p, den- 
jenigen der zweiten Reihe zukommen. Ferner bedeuten « und «@ willkürliche Para- 
meter; und zwar regelt « (von 0 bis ©) die Wellenlänge der (sin-förmigen) Verbiegung 
eines jeden Wirbels längs seiner (zur z- Achse parallelen) »Ruh«-Achse, (von 0 bis z) 
die Aenderung der Amplitude dieser Verbiegung von Wirbel zu Wirbel, oder — wie man 
auch sagen kann — die Wellenlänge der Störung längs der Straße. Schließlich sind 
die Doppelvorzeichen in den Exponenten so zu verstehen, daß jeweils alle vier mög- 
lichen Kombinationen zugelassen sind, dagegen das Doppelvorzeichen vor {jo und [20 80 
zu wählen ist, daß für Kombination — + und — — im Exponenten das negative, für 
Kombination + — und — + das positive Vorzeichen gilt. Der Grund hierfür liegt darin, 
daß unser Ansatz imstande sein muß, den Randbedingungen an den beiden die Flüssig- 
keit begrenzenden Ebenen 2 = —+ d/2 zu genügen. 

Diese Randbedingungen besagen, daß keine Flüssigkeit durch die Begrenzungs- 
platten ein- oder austreten darf; und dies wird dann und nur dann erreicht, wenn die 
Wirbelachsen in ihren Durchstoßungspunkten auf den Ebenen z2= + d/2 senkrecht 
stehen und durch diese Ebenen hindurch nach beiden Seiten hin ins Unendliche fort- 
gesetzt gedacht werden. Infolge der so erreichten spiegelbildlichen Symmetrie beiderseits 
der Platten ist ohne weiteres klar, daß der Grenzbedingung hinreichend genügt wird. 
Notwendig jedoch ist das Senkrechtstehen auf alle Fälle darum, weil nur so die im 


> 
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Wirbelkern nach Art eines festen Körpers rotierenden Teilchen daran verhindert werden, 
eine Komponente der Geschwindigkeit normal zur Begrenzungsplatte zu besitzen. Wenn 
wir also wieder darauf ausgehen, die Störungen in »harmonischer« Form anzusetzen, mit 
dem Vorbehalt, daraus dann nach Fourier die allgemeinsten möglichen zusammenzusetzen 


(wie dies schon beim Einzelwirbel geschah), so werden wir mit den Grenzbedingungen 
im Einklang bleiben, wenn wir schreiben: 


— &, cos (up): cos | 


(up): sin (@z) 


(gültig sowohl für die erste wie für die zweite Reihe, wenn die entsprechenden Indizes 


ı und g bzw. 2 und « angefügt werden) — und außerdem den zulässigen Werten des 


Parameters « die Beschränkung auferlegen: 
d 


wo n eine beliebige ganze Zahl von 0 bis ®. 


Diese Beschränkung bedeutet offenbar anschaulich, daß nur ganze Vielfache einer 
vollen Wellenlänge zwischen den Platten liegen dürfen, wie dies aus obiger Forderung 
des Senkrechtstehens ohne weiteres folgt und aus unserm Ansatz (Ila) für den Fall, daß 
— wie es bei uns angenommen ist — die Ebene 2— 0 die Mediane zwischen den Be- 
grenzungsplatten darstellt. 

In komplexer Form geschrieben, ist nun offenbar der Ansatz (lla) mit dem obigen 
Ansatz (II) identisch einschließlich der dort gegebenen Zeichen-Vorschriften. 


3. Bedeutung der Störungen für die Selbstinduktion. Bevor wir daran 
sehen, das vollständige Diiierentialgleichungssystem unseres Problems anzugeben, mag 
noch kurz darauf hingewiesen werden, daß die im 1. Abschnitt behandelten Kelvinschen 
Störungstypen ?=0 und ?—1 tatsächlich die einzigen für die Straße bedeutungsvollen 


sind, insofern als mit ihrer Hilfe die im vorigen Abschnitt charakterisierte Störung der 
Straße vollständig dargestellt werden kann. 


Es handelt sich hierbei, wie wir gesehen haben, um eine Störung, die — abgesehen 
von der Amplitudenänderung von Wirbel zu Wirbel — einen jeden Wirbel vermöge der 
transversalen Störungen und schlangenförmig verbiegt, außerdem aber, vermöge 
der longitudinalen Störung [, ihn auch noch abwechselnd »staucht« und »dehnt«, und 
zwar gerade dort maximal, wo die transversale Elongation ihre Minima hat, was daraus 
hervorgeht, daß bei { der sin steht, wo bei & und 7 der cos. Diese Störungsform, aus 
der in Fourierscher Weise jede beliebige iiberhaupt mögliche Störung der Straße auf- 
gebaut werden kann, ist so angesetzt, als bestünden die Wirbel nur aus ihren Achsen, 
d.h. als wäre der Radius des Wirbelkerns Null. Dies darf, wie wir bereits in der Ein- 
leitung bemerkten, bei verbogenen Wirbeln nicht mehr angenommen werden, wenn unsere 
Ausdrücke endlich bleiben sollen; wohl aber dürfen wir annehmen, daß von den Nachbar- 
wirbeln aus gesehen, ein jeder Wirbel auf seine Achse zusammengeschrumpft erscheint. 
Und daher werden wir zwar für den Wirbel, dem das »Aufelement« angehört, eine end- 
liche Dicke annehmen, aber nur solche Verbiegurgen als für die Straße wichtig betrachten, 
bei welchen die kreisförmige Gestalt des Kernquerschnitts (mindestens in erster 
Näherung) erhalten bleibt. Man sieht nun unmittelbar, daß die so charakterisierten 
Störungstypen gerade mit denjenigen übereinstimmen, welche im 1. Abschnitt (nach 
Kelvin) durch /=0 und i= 1 gegeben waren. Denn einerseits bleibt bei diesen und 
nur bei diesen der Kreisquerschnitt der Wirbel in erster Näherung erhalten, andererseits 
liefert der Fall@2—0 die longitudinale Komponente ({) unserer Störung (d.h. die 
Stauchung und Dehnung der Wirbelfäden in Richtung ihrer Achsen), der Fall @—=1 die 
transversale Komponente (d. h. die eigentliche Verbiegung; E und 7). Alle übrigen 
Störungsformen bestehen in einer »Riffelung« der Oberfläche des Wirbelkerns und ver- 
schwinden im Mittel über eine volle Umdrehung desselben (mit Winkelgeschw. ®); sie 


sind also für die Straße bedeutungslos, da ja» sehr groß ist gegen alle durch die 
Störung bewirkten Frequenzen. 


Wie wir bereits sahen, besteht der eine Störungstypus hinsichtlich seines zeit- 
lichen Verlaufes (der Wirbel ist zunächst noch isoliert gedacht) in einer Oszillation 
eines jeden Wirbelelements (speziell also des »Aufelements«) parallel z mit der Fre- 
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juenz n (0); der andere Störungstypus in einer Rotation des verbogenen Wirbels um 
seine »Ruh«-Achse herum mit der Winkelgeschwindigkeit nr. - ı). 

Diese beiden Ausdrücke n;-.0, und n«-—ı) wollen wir hier zum Schluß noch auf 
unsere (im vorigen Abschnitt eingeführten) Bezeichnungen umschreiben, da diese von 
den Kelvinschen abweichen. 


Wir nennen: den Radius des Wirbelkerns ou (Kelvin a) 
den Parameter der Störung längs z-Achse « (Kelvin m), 
die Winkelgeschwindigkeit des Kerns ® (unverändert gegen Kelvin), 
die Frequenz der z-Osziallation A, (Kelvin nv-—0,), 
die Winkelgeschwindigkeit des verbogenen Wirbels um seine Ruh- 
Achse A, (Kelvin 1). 
Somit haben wir endgültig für die beiden Größen, welche den von der Selbst- 
induktion herrührenden Bestandteil der Störungsfirequenzen eines der gestörten Straße an- 
gehörigen Wirbelelements bestimmen, die Ausdrücke: 


20009 


A= (g= 2,4 usw., siehe 1. Abschn.!) | 
1 

Aı= — (In + 0,366) \ 
2 0, 

4. Das Differentialgleichungssystem des Problems. Nach diesen Vorbereitungen 
gehen wir daran, das Differentialgleichungssvstem unseres Problems aufzustellen. Wir 
betrachten — wie bereits erwähnt — ein dem ten Wirbel angehöriges Element unter 
dem Einfluß sämtlicher übrigen Wirbelelemente der ganzen Straße. Dabei sei für dieses 
»Aufelement« der z-Wert z—= 0 angenommen, was olienbar keine Beschränkung der All- 
gemeinheit bedeutet, da man jederzeit durch bloße Vertauschung der Integrationsvariabeln 
z gegen (z — z,), wo 24 = des Aufpunkts] zu einem beliebigen anderen Wert übergehen 
kann, ohne daß sich dadurch unsere folgenden Gleichungen ändern. 

Wir setzen an das Biot-Savartsche Gesetz der Wirbelbewegung, dem gemäß 
die in einem Aufpunkte A induzierte Geschwindigkeit v, sich darstellt durch die Gleichung: 

"fr, as] 

4 

zu integrieren über alle an der Induktion beteiligten Wirbelachsen-Elemente d®, d.h. in 
unserem speziellen Falle, wo eine unendliche Anzahl diskreter Wirbel existiert, einer- 
seits zu integrieren über die ganze Länge jedes einzelnen Wirbels, andererseits zu 
summieren über alle vorhandenen Wirbel. Dabei können wir — wie im 2. Abschnitt 
erläutert wurde — um den Randbedingungen an den Begrenzungsplatten zu genügen, 
die Integration zwischen den Grenzen + d/2 und — d/2 ersetzen durch eine solche 
zwischen unendlichen Grenzen, und dafür die Beschränkung des die Deformation des 
Wirbels regelnden Parameters auf eine diskrete Reihe erlaubter Werte einführen [s. Gl. (1II)]. 
Bezeichnen wir nun noch der Kürze halber mit y=e*ivr,#ir: = e*ivp, 
die Funktion, welche die Zurückführung der unendlich vielen unabhängigen Störungs- 
größen auf 6 unabhängige Variable leistet [vergl. Gl. (II) und (Ila)], ferner mit dem In- 
dex A die Zugehörigkeit zum Aufpunkt, mit dem Index 1 und 2 diejenige zur 1. bzw. 2. 
Wirbelreihe, (wobei zu beachten ist, daß 214 = 224 — 0 nach unseren obigen Festsetzungen), 
und erteilen wir jedem Wirbelelement die ihm zukommende infinitesimale Verrückung 
(£n{), so schreiben sich die Bewegungsgleichungen des » Aufelementes« wie sie das Biot- 
Savartsche Gesetz liefert: 


+ % . 09 p “ 
in, (rı + J rı)? 
| (V) 
( ( 
+ (: +2 — [E14— (29p + Eop)] 
u PA 
(rıa + r2)° \ | 
analog je eine Gleichung für (v,)ı« und (v,)ı, durch entsprechende zyklische Vertauschung 
der Größen &, y, z und S, », {, wobei jedoch die Sonderrolle der z Koordinate bewirkt, 
08 


geht und umgekehrt; dazu kommen noch die 3 weiteren Gleichungen für 
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welche aus den drei ersten erhalten werden durch Vertauschung von Index »1« ge- 


sen Index »2« sowie Summationsindex »p,« gegen »p.« und umgekehrt. Dabei ist 
als Abkürzung eingeführt: 


+4 rı) = +4 = + + &,,)) + yı4) — (yı,» 


sowie: 


+ Ira)=(rsı + — (0, + + + 114) — + 
Anm.: [In diesen Gleichungen (V) ist rein formal auch der von der Selbst- 
induktion herrührende Bestandteil der Störungsgeschwindigkeiten mit enthalten in Ge- 
stalt desjenigen Summengliedes, für das der Summationsindex gerade den Wert annimmt, 
welcher dem den Aufpunkt enthaltenden Wirbel selbst zukommt. In diesem Gliede wird 
im Nenner r=0, weshalb bekanntlich die Behandlung nach dem Biot-Savart- Gesetz 


in diesem Falle mit Schwierigkeiten verknüpft ist. Wir wollen statt dessen auf die im 
3. Abschn. erhaltenen Resultate zurückgreifen (s. S. 358) | 


Entsprechend dem Sinne der infinitesimalen Störung seien nun die Ausdrücke 


auf der rechten Seite des Systems (V) nach Taylor entwickelt und nur die ersten 
Potenzen der £7{ berücksichtigt. 


Bei (v.)ı4 und (v,)4 tritt je ein gleichlautendes, von diesen Störungen unabhän- 
giges Glied auf, die gemeinsame Fortschreitungsgeschwindigkeit der gesamten Straße, 
naturgemäß von genau derselben Größe wie beim zweidimensionalen Fall (denn wir gin- 


gen ja von genau der gleichen ungestörten Grundkonstellation aus), nämlich in Kärmän- 
scher Schreibweise: 


Das entsprechende Glied bei (v,)ı4 und (v,)»a wird Null (s. Kärmän ]l. e.), und bei (v.)ı4 
und (v.)s4 treten überhaupt von vornherein nur Glieder von der Ordnung $n{ auf (weil 
ja eine Geschwindigkeit parallel z erst durch die störende Verbiegung erzeugt wird). 


Was nach Abspaltung dieser endlichen Bestandteile auf der linken Seite übrig 


bleibt, sind die reinen infinitesimalen Störungsgeschwindigkeiten des Aufpunktes in beiden 
Reihen: 


dt dt dt 
Auf der rechten Seite stehen die Koeffizienten der Taylor-Entwieklung multipliziert 
mit den Größen S$ı bis (. Es ist aber dabei wie wir sahen (s. oben, Anm.) bei allen Summa- 
tionen das vom »Aufwirbel« herrührende Glied auszulassen und durch denjenigen Be- 


standteil der Störungsgeschwindigkeiten zu ersetzen, welchen wir aus Kelvins Ansätzen 
für die Selbstinduktion erhielten. 


Dieser Anteil der Selbstinduktion an der Störungsbewegung des Auipunkts 
besteht, wie wir im 3. Abschn. sahen, in einer Rotation um die ungestörte Wirbelachse 
(in der &y-Ebene) mit der Winkelgeschwindigkeit Aı [siehe Gl. (IV)]| und einer davon 
unabhängigen Oszillation parallel zur z-Achse mit der Frequenz Au. Wir beschreiben 


diese Verhältnisse durch die — für den isoliert gedachten Aufwirbel gültigen! — Glei- 
chungen: 
+ Aı na } 
(VI). 
Cu 


Die übrigen Taylor-Koeffizienten sind gewisse Integral-Summen, welche lediglich 
von der Verbiegung und gegenseitigen Entfernung der Wirbel abhängen. Es zeigt sich, 
daß man sie auf 13 unter sich verschiedene Größen ZR, S,... reduzieren kann, die, um 
den Gedankengarng nicht zu unterbrechen, im Anhang [unter a)| zusammengestellt sind. 


Indem wir also den durch Gl. (VI) beschriebenen Effekt mit dem von den übrigen 
Wirbeln der Straße herrührenden superponieren, erhalten wir das endgültige Difierential- 
gleichungssystem unseres Problems: 


24* 
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- — +0% 
+Plo + (8, 530 + 
— T — (Sı +0%o . 
— (5, + Tnıo — QSıo . 
Z — Ynıo — Xn20 


Von den Koeffizienten dieser Gleichung sind reell: 81,83: U,V; X,Y; P,Q; 
rein imaginär: 7, O0, Z [siehe Anhang a)). 


5, Berechnung des Wirbelkern-Radius. Für die Stabilıtäts-Diskussion, d.h. 
die Entscheidung der Frage, unter welchen Bedingungen die charakteristische Gleichung 
des Systems (VII) imaginäre Wurzeln besitzt, ist es offenbar wesentlich, nicht nur die 
Koeffizienten Rı, A, usw. bis Z unter sich, sondern auch mit 4, und A, numerisch zu 
vergleichen für alle möglichen Werte der Störungsparameter ““ und «. In der erstge- 
nannten Gruppe von Koeffizienten kommt nun der Radius 0, des Wirbelkerns nicht vor, 
da sie nur von den Wirkungen der Wirbel aufeinander herrühren, wohl aber ist 0, in 
Ao und A, enthalten, weil letztere Größen ja durch die Wirkung des »Auf-Wirbels« auf 
sich selbst gegeben sind. Solange nun über die Größe des Wirbelkern-Radius im Ver- 
gleich zu den übrigen Längenabmessungen der Straße (d.h. also den Wirbel-Abständen 
bzw. h und /) nichts bekannt ist, können die beiden genannten Koeffizienten-Gruppen in 
keinerlei Größenbeziehung zueinander gesetzt werden. Es muß daher versucht werden, 
0, aus den pbysikalischen Bestimmungsstücken der Wirbelstraße zu berechnen. 


Der Gedankengang welcher der Berechnung des Wirbelkern-Radius zugrundeliegt, 
ist folgender. Es kann einerseits nach den von Kärmän') durchgeführten Reehnungen 
der Widerstand, den ein Körper in einer Flüssigkeit erfährt, vollständig ausgedrückt 
werden durch Größen, die auf die Wirbelstraße Bezug haben, und zwar tritt in diesen 
Ausdrücken der Radius o, des Wirbelkerns nicht auf. Den letzteren Umstand kann 
man sich plausibel machen, indem man sich vergegenwärtigt, daß dieser Widerstand ge- 
wonnen wird als propoıtional der pro Zeiteinheit beim Schleppen des Körpers erzeugten 
(und eben in der Wirbelstraße invertierten) Impulsmenge, und daß ja der Impuls eines 
einzelnen Wirbelpaares, etwa in der x-Richtung, naturgemäß ganz unabhängig ist von 
der »Dicke« der Wirbel, da er ja (als einseitig gerichtete Größe!) seinen Wert sogar dann 
beibehält, wenn bei gegebenem Wirbelmoment die Wirbel unendlich dinn werden, und 
somit »in deren Achse« die Teilchen unendlich rasch rotieren. Auf Grund dieser Wider- 
standsformel kann man weiterhin die Arbeit, die in gegebener Zeit zum Schleppen des 
Körpers geleistet werden muß, wiederum ausdrücken als Funktion der charakteristischen 
Konstanten unserer Wirbelkonfiguration. Diese Arbeit muß nun offenbar gleich 
sein der kinetischen Energie, welche das in der gleichen Zeit entstandene Stück 
der Wirbelstraße besitzt. Denn energetisch kann der Schlepp-Prozeß in unserer reibungslosen 
Flüssigkeit ja garnicht anders interpretiert werden, als daß einerseits Arbeit beim Schleppen 
verbraucht wird, andererseits dafür kinetische Energie der Wirbelstraße erscheint; ein 
weiterer Effekt tritt nirgends auf. Nun verhält sich aber hinsichtlich der Abhängigkeit 
von der Wirbeldicke die kinetische Energie eines Wirbelsystems prinzipiell anders als 
dessen Impuls. Denn es ist unmittelbar evident, daß unendliche Rotationsgeschwindig- 
keiten »in«e der Achse (d.h. eben für 9 = 0) auch die kinetische Energie unendlich 
machen müssen (da sich ja hier unendlich große, aber entgegengesetzt gerichtete Geschwin- 
digkeiten nicht, wie bei Bildung des Impulses, zerstören). Es muß also in den Aus- 
druck für die kinetische Energie der Radius 0., des Wirbelkerns durchaus eingehen. 
Daraus folgt schließlich, daß die geforderte Gleichheit von geleisteter Arbeit und erzeugter 
Energie nur dann bestehen kann, wenn 2, als Funktion der charakteristischen Konstanten 


der Straße (eventl. auch noch der Schleppgeschwindigkeit) einen ganz bestimmten Wert 
hat. Diesen gilt es nun zu ermitteln. 


Zunächst schreiben wir den Ausdruck für die pro Zeiteinheit geleistete Arbeit an, 
welcher aus den Kärmänschen Formeln (l. ec.) ohne weiteres folgt. Es bezeichne: 


') Phys. Zeitschr. 13. Jahrg. 1912, S. 49 bis 59. 


- 
== 
!o = 
ao = 
= 


ch. 
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U die Geschwindigkeit ( x) mit der der Körper durch die Flüssigkeit gezogen 
wird; 

o die Dichte der Flüssigkeit; 

u die Fortschreitungsgeschwindigkeit der Straße; 

k das Verhältnis A:l; 

W den Widerstand beim Schleppen mit der Geschwindigkeit U. 


Dann ist (A = Arbeit pro Zeiteinheit): 
in 


Die kinetische Energie des pro Zeiteinheit abgelösten Stücks der Straße ge- 
winnen wir durch folgende Ueberlegung. Kennt man für ein gegebenes Wirbelsystem 
die »Stromfunktion«, so läßt sich durch letztere die kinetische Energie des Wirbel- 
systems ausdrücken gemäß der Formel '): 


wo unter %; die Wirbelstärke des ’-ten Wirbels, unter ”', die Stromfunktion, genommen 
»am Örte« des Wirbels, zu verstehen ist. 

Für die Wirbelstraße ist diese Stromfunktion aus dem von Kärmän (l. ce.) ange- 
gebenen »komplexen Potential«: 


2n sin IT 4 Y» 


also +9), n—:=|- e)+il, 


„als imaginärer Teil zu errechnen. Wir erhalten: 


Sin? - cos? Xı + Yı - sin? X, 
Sin‘ Y,- cos? + Y,- sin? 


(X) 


Am 


indem wir abkürzend einführten: 


x): all; Y=(- + y) all; y) al. 


Die durch (X) gegebene Größe w ist — wie die Ableitung bei Kärmän l.c. zeigt — 
die Stromfunktion bezogen auf ein mit der Wirbelstraße mitbewegtes, d.h. mit der Ge- 
schwindigkeit u fortschreitendes Koordinatensystem. Um deren Wert in bezug auf ein 
(gegen die unendlich weit vor dem durchgezogenen Körper : befindliche Flüssigkeit) 
»ruhendes« Koordinatensystem zu erhalten, haben wir dem % noch ein additives Glied 
hinzuzufügen derart, daß dessen Ableitung nach y gerade die negative ©-Komponente der 
Geschwindigkeit im ruhenden System, also gerade — u liefert. Derjenige Wert unserer 
Stromfunktion, den wir (am Orte der Wirbel gebildet!) in Gl. (IX) einzuführen haben, 
ıy8 Sin? eos? + Bo]? Y3 - sin? X, 


in 
Exakt an einer Stelle gebildet wo ein Wirbel sitzt, wird das zweite Glied von (Xa) 

unendlich, wie es sein muß, da ja dann, wie wir sahen, die kinetische Energie bereits 
für einen einzelnen Wirbel unendlich ist. Nun ist aber der Sinn des Energieansatzes 
(IX) — wie aus seiner Ableitung (s. Lamb, |. c.) hervorgeht — der, daß die gesamte 
Energie der mit Wirbeln durchsetzten Flüssigkeit (inclusive der Energie der Wirbelgebiete) 
durch eine Integration über das wirbelfreie Gebiet dargestellt wird; und in der Tat 
gibt ja die Funktion ® die Geschwindigkeit nur im wirbelfreien Gebiet. Unsere Frage 
nach der Größe des Kernradius o, kann also auch so formuliert werden: wie weit hat 
man das wirbelfreie Gebiet bis an die Wirbelachsen heran zu erstrecken, damit die 
kinetische Energie des Systems gerade den »richtigen« (durch Gleichsetzung mit (VIII) 
zu erhaltenden) Wert bekommt? Wir dürfen somit (ebenso wie unsere vorhin erwähnte 
Integration, so auch) unsere »Entfernungen« soweit sie in den Formeln für % auftreten, 
nur bis an die Grenze des wirbelfreien Gebietes erstrecken, haben also für die »Punkte«, 
in denen Wirbel sitzen, anstatt der exakten Werte 


1; (wo ] bel. ganze Zahlen) 


I) Vergl. z. B. Lamb, Hydrodynamies, 4. Aufl. (1916). Kap. VII, Art. 157. 
?) Lamb, Hydrodynamies, 4. Aufl. Kap. VII. Art. 156. 
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die »ungefähren« Werte 
+, 
einzusetzen, welche eigentlich keine »Punkte« mehr, sondern eben kleine Gebiete vom 
Radius 0 kennzeichnen. Es sei natürlich wieder 0% <& h bezw. X! angenommen wie 
bisher stets. Damit ist die kinetische Energie in ihrer Abhängigkeit von », dargestellt. 
Wenn wir auf der rechten Seite von Gln. (X) Zähler und Nenner nach dem so einge- 
führten ou entwickeln und im Nenner nur die erste Potenz berücksichtigen, während wir 
im Zähler die mit (o,/l) multiplizierten Glieder gegen die mit (Ah/l) multiplizierten streichen, 
so erhalten wir: 
x Sol? (nk); 


== In 


dies an der betreffenden Stelle (als 2. Glied) in (Xa) eingeführt und alsdann das 1. Glied 


für die Orte gebildet, an denen Wirbel sitzen, gibt uns schließlich den Wert von w,, 
weleher in unsere Energieformel (IX) einzusetzen ist, als: 


Um die in der Zeiteinheit erzeugte kinetische Energie zu erhalten, haben wir die 


Summe (IX) über die Anzahl N der pro Zeiteinheit abgelösten Wirbel zu erstrecken; diese 
aber ist im Mittel: 


2 (U-ıu) 
Daher bekommen wir schließlich (Z#.. — Energie pro Zeiteinheit): 
U—u 0%” 2nh 


ıYs 


FR 
Gleichsetzen von (XI) mit (VIII) ergibt: 


2; (Si? 7 \ 
ı ys 7109 ) 2 

|\ 


= (7 k) ! u) Zrl ı) . . (XI). 


Nun setzen wir speziell den Kärmänschen Wert für das Verhältnis k=h:l ein, 
nämlich 


und daraus: 


k=0283, 


ferner den (s. Kärmiän |l.c) empirisch bestimmten Wert 


u 


und erhalten so 


/ 


2 und U-l- 141; 


daraus folgt dann hinreichend genau: 


1 
I . . 0,01 ] 
7T 


und: 


Die Geschwindigkeit des geschleppten Körpers hat sich also ganz 
herausgehoben, der Kernradius der in einer Kärmänschen Straße zu- 
sammengestellten Wirbel erscheint als eindeutig gegeben durch die Längen- 
dimensionen dieser Straße, d.h. also infolge der (h/!)-Bedingung, eindeutig 
gegeben durch die Breite A des Körpers, 


ein, 


Oktober 1928 Schlayer, Ueber die Stabilität der Kärmänschen Wirbelstraße 363 


6. Die Stabilitäts-Diskussion. Wir haben uns nun die Frage vorzulegen, unter 
welchen Bedingungen die durch das Gleichungssystem (VII) dargestellte Bewegung der 
gestörten Wirbelstraße eine »stabile«, d.h. eine periodische Bewegung ist. Die Lösungen 
eines solchen Systems homogener linearer Differentialgleichungen erster Ordnung, lassen 
sich darstellen als lineare Funktionen von e’‘, wo A; die Wurzeln der »charakteristischen 
Gleichung«, und die Stabilitätsuntersuchung läuft hinaus auf eine Untersuchung, unter 
welchen Bedingungen nur periodische Lösungen auftreten, d. h. sämtliche }, rein- 
imaginär werden. 


Führen wir die Abkürzungen ein: 
[2 + (U+ A)]=Kı 


. (VIla), 
+V)=L;: 
so lautet die charakteristische Gleichung unseres Systems (VII): 
Kı —T’ +0 
| -K +-P +26 +T 
\ 
— 7 +0 (— 4) — Kı 
— X (—iAo 4) 


Ausgeführt gibt dies eine Gleichung 6. Grades für A, welche wir in der Form 
schreiben wollen: 


dı —- dı —- da 2? + da 2 = . . . (XIV). 


Die Koeifizienten dieser Gleichung setzen sich in komplizierter Weise aus den 
Elementen der Determinante (XIII) zusammen. Ihre genaue Kenntnis ist natürlich not- 
wendig, wenn exakt berechnet werden soll, welche stabile bezw. instabile Bewegung von 
irgend einer gegebenen Störung ausgelöst wird. Eine solche Berechnung ist jedoch über- 
flüssig, da man sich ohne Gefährdung des prinzipiellen Ergebnisses auf Größenordnungs- 
Betrachtungen beschränken kann. Immerhin erweist es sich als vorteilhaft, über die 
roheste erste Annäherung hinaus (die sich direkt aus der unaufgelösten Determinante 
(XIII) ergibt), noch zu einer zweiten Näherung der Gl. (XIV) fortzuschreiten, da erst 
hier die schwache Instabilität zutage tritt, welche dem dreidimensionalen Problem 
eigentümlich ist. 


Für diese Rechnung in zweiter Näherung braucht man nun ebenfalls eine ziemlich 
genaue Kenntnis der Koeffizienten der Gl. (XIV). Da es sich aber um recht unüber- 
sichtliche Ausdrücke handelt, so sollen sie erst im Anhang zusammengestellt werden. 
Dort finden sich auch die Einzelheiten der Stabilitäts-Diskussion angegeben, soweit sie 
zur Kontrolle der Ergebnisse erforderlich sind. Die eigentliche Rechnung in extenso zu 
entwickeln, wäre viel zu weitläufig und ein durch das Ergebnis nicht gerechtfertigter 
Aufwand, da im Wesentlichen nur eine Bestätigung der Resultate Kärmäns heraus- 
kommt. Hier im Text soll nur der Gedankengang dargelegt werden, der zur ersten 
bzw. zweiten Näherung führt. 


Betrachten wir die numerische Darstellung der Größen A, Aı; Rı: A, usw. bis Z 
(Anhang, Abb. 3 und 4!), welche die Elemente der Determinante (XIII) bilden, so er- 
kennen wir, daß es zwei Bereiche von Störungen „ibt, für welche sich besagte Deter- 
minante in erster Näherung erheblich vereinfachen läßt. Dies sind erstens der Bereich 
»Kärmän-ähnlicher« Störungen, d. h. schwacher Verbiegung (kleine Werte des Para- 
meters X), zweitens der Bereich sehr starker Verbiegungen. Im ersteren Falle nämlich 
sind alle von der &- und 4 Komponente des Wirbelmoments herrührenden Koeffizienten 
im Gleichungs-System (VII) sehr klein, verglichen mit den Koeffizienten %, 5 und T, 
welche bereits im zweidimensionalem Problem vorkommen'), so daß wir anstatt (XIII) 
schreiben Können: 


) Bei Kärmän l.ce. durch A, B, C, bei Lamb durch «, ß, y bezeichnet. 
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(—A) R . T NS 
R (— +-T 
+7 — BR (— A) 
oder: 
R — T 
R 
(— 4) — R 
| — R (— 4) 
was die Wurzeln liefert: —= + und T-+ -S?, d. h., da 7’ rein- 


imaginär und alles ürige reell ist, Stabilität, wofern nur die Kärmänsche Bedingung: 
erfüllt ist. 

Im entgegengesetzten Grenzfall sehr starker Verbiegungen (k’ > 1) kann die 
Determinarte in erster Näherung sogar noch einfacher geschrieben werden, da hier 4, 


und 4A, weit über alle anderen Elemente überwiegen (s. Kurven I u. 2 in Abb. 3 des 
Anhangs!) nämlich als: 


(+ 4 NIIb), 


was die Wurzeln liefert: 


= + und hiss = + | 
d.h. Stabilität unter allen Umständen. 

Physikalisch besagt dies, daß bei schwacher Verbiegung der Wirbel im wesent- 
lichen ihre gegenseitige Beeinflussung bemerkbar ist, welche bewirkt, daß Stabilität 
zwar möglich ist, aber nur für einen bestimmten Wert des Verhältnisses h:!; bei starker 
Verbiegung dagegen überwiegt die Wirkung der Selbstinduktion; d.h. die Präzession 
und Oszillation, die jeder verbogene Wirbel an sich selbst hervorruft, wird so rasch, daß 
sie durch die viel langsamere Einwirkung benachbarter Wirbel nicht mehr merklich 
beeinflußt werden kann, die Stabilität also von der Anordnung der Wirbel in der Straße 
in erster Näherung unabhängig wird. 

Zwischen diesen beiden Grenzgebieten aber liegt das kritische Gebiet wo die gegen- 
seitire Induktion mit der Selbstinduktion in Resonanz treten kann, weil die von 
beiden erzeugten Oszillations- und Rotations-Frequenzen von gleicher Größenordnung 
sind, d. h. mathematisch gesprochen, weil sämtliche Elemente der Determinante (XIII) dem 
gleichen Größenordnungsbereich angehören. Hier darf also die charakteristiche Gleichung 
nicht so radikal vereinfacht werden wie dies durch Streichung der meisten Glieder in 
den beiden obigen Fällen geschah. Trotzdem aber wird man vermuten können daß auch 
hier noch die Ausdrücke für die Störungs-Freqjuenzen 4, von ähnlicher Form sind wie 
die aus Gleichungen (XIIla und b) folgenden. 

Streicht man nämlich in der Determinante (XIII) alle Elemente außer As; Kı,K%; 
L,,I12 und 7, so erreicht man einerseits, daß die so gewonnene Gleichung in den beiden 
Grenzfällen schwacher und starker Verbiegung in die Ausdrücke (XIlIa) und (XIIIb) 
übergeht, andererseits daß für jede beliebige Störungsform Größen stehen bleiben, welche 


dem Betrag nach größer sind als alle gestrichenen (siehe Abb 3 und 4 des Anhangs!); 
und man erhält: 


— t 4°; bise T- — 


Diese Ausdrücke für die Störungsfrequenzen erweisen sich in der Tat als in 
grober Näherung ausreichend für die meisten Störungstypen. Es gibt jedoch einen 
beschränkten Größenbereich des Verbiegungsparameters k', für welchen sie den Tatsachen 
nicht mehr entsprechen, wo vielmehr Selbstinduktion und gegenseitige Induktion derart 
in Resonanz treten, daß eine (allerdings äußert langsame) instabile Bewegung resultiert. 


4 


eh. 


ie 
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Der Beweis, dessen arithmetische Einzelheiten im Anhang zusammengestellt sind, 
relingt auf folgendem Wege: Zunächst wird eine Störungsiorm betrachtet, für welche die 
Determinante (XIII) sich erheblich Meregen ohne daß irgendwelche Vernachlässigungen 
nötig sind. Es ist dies der Fall «= 7/2 (d. h. eine Störung bei der die Wirbel der ersten 
Reihe alle mit gleicher Amplitude verbogen sind — und zwar so, daß ihre Spurpunkte 
in der Ebene z= 0 abwechselnd in positiver und negativer «- und %- Richtung elongiert 
sind, während diejenigen der zweiten Reihe alle in ihrer Ruhelage bleiben), welcher die 
Gleichung 6. Grades in A (XIV) auf eine solche dritten Grades in A” reduziert. Die 
Wurzeln dieser kubischen Gleichung lassen sich angeben und erweisen sich als durchgehends 
reinimaginär für alle Grade der Verbiegung (Parameter k') außer für den Fall wo 
die Selbstinduktion auf die gegenseitige Induktion in besonderer Weise abgestimmt ist, (nämlich 
so, daß annähernd 4° = T”, siehe Anhang unter e)'). Hier werden die }, teilweise 
komplex, also die Bewegung instabil. Aber diese Instabilität ist äußerst schwach, 
die betreffende aperiodische Bewegung sehr langsam im Vergleich zu den gleichzeitig (oder 
auch in den durchaus stabilen Nachbarbereichen) stattfindenden stabilen Oszillationen. 
Ueberdies wird die Bedeutung dieser Instabilität noch dadurch erheblich herabgesetzt, daß 
sie ja nur in zwei engen Bereichen innerhalb der unendlichen Menge möglicher Störungs- 
formen überhaupt auftreten kann. 

Der Uebergang auf Störungen mit beliebigen w-Werten gelingt durch eine Nähe- 
rungsmethode, die darauf beruht, daß sich die allgemeine Gleichung 6. Grades (XIV) 
auffassen läßt als eine Modifikation der speziellen kubischen Gleichung. Der Nachweis 
der Berechtigung dieses Verfahrens sowie dessen wichtigste Schritte finden sich im 
Anhang zusammengestellt (unter d)). Es ergibt sich jedoch gegenüber dem Spezialfall 
«= n/2 nichts wesentlich neues, nämlich auch hier wieder Stabilität außer im Bereiche 
Au” ” T” usw, wo eine äußerst schwache Instabilität hinzukommt. 


7. Prinzipielles zur neuen Stabilitätsbedingung. In den Ergebnissen des 
vorangehenden Abschnittes sind zwei Stabilitätsbedingungen enthalten. Die eine ist die 
alte Kirmänsche, welche aus der Störung ohne Verbiegung (k’—= 0) folgt und, als von 
vornherein notwendige Voraussetzung jeglicher Stabilität, bereits der Berechnung des 
Wirbelkerns (Abschn. 5) sowie derjenigen der Koeffizienten unseres Differentialgleichungs- 
Systems (VII) — s. Anhang, Abb. 3 und 4 — zugrunde lag. Neu hinzu kommt eine 
zweite Bedingung, die den Plattenabstand d betrifit. Es läßt sich jedoch beweisen, daß 
diese neue Bedingung praktisch bedeutungslos ist, und wegen der prinzipiellen Wichtig- 
keit der Sache soll dieser Beweis hier kurz angegeben werden. 

Besagte Bedingung ist aus den angegebenen Ausdrücken für die Störungsfrequenzen 
darum nicht ohne weiteres ablesbar, weil in der obigen Diskussion immer nur vom Ver- 
biegungsparameter k' und nicht vom Plattenabstand d die Rede war. In der Tat ist die 
Verbiegung durchaus das Primäre, denn sie allein erlaubt es, die Länge der Wirbelfäden 
überhaupt ins Problem hereinzubringen, d. h. zwei in der z Richtung voneinander ent- 
fernte Punkte zu unterscheiden und ihren Abstand zu messen. Von der Verbiegung aus 
eelangt man über die Randbedingungen zum Abstand der Begrenzungsplatten als einer 
Strecke, die nur ganze Vielfache einer Wellenlänge der sinusförmig verbogenen Wirbel 
betragen darf. Darum ist der Plattenabstand zunächst durchaus nicht eindeutig in unseren 
Ausdrücken enthalten, denn Verdoppelung, Verdreifachung usw. desselben bei gleichzeitiger 
Verdoppelung, Verdreifachung usw. der zwischen den Platten gelegenen Anzahl von 
Störungs Wellenlängen führt zu lauter hinsichtlich ihrer mathematischen Behandlung 
völlig identischen Problemen. Trotzdem wird durch einen gegebenen Plattenabstand 
unter sämtlichen zulässigen Wellenlängen eine ganz bestimmte als die größte ausge- 
sondert (nämlich diejenige, die gleich dem Plattenabstand ist, d.h. n—=1 in @). (II); 
s. a. Abb. 3 des Anhangs!) und auf diese Weise die Eindeutigkeit wiederhergestellt. Diese 
Verbiegung von gıößter Wellenlänge besitzt gleichzeitig die schwächste Krümmung, 
welche ein Wirbel zwischen den Platten annehmen kann. Nun besteht — wie im vorigen 
Abschnitt gezeigt — für gewisse Grade der Wirbelkrümmung eine schwache Instabilität, 
und unter diesen instabilen Krümmungsgraden gibt es einen größten, jenseits dessen die 
Wirbel absolut stabil werden Um absolute Stabilität zu erzielen ist es daher notwendig 


!) Auch die Abstimmung der longitudina'en (d. h. z2-) auf die transverslae (d. h. x&- und y-) 
Komponente der Selbstinduktion (derart, daß Ay ” ,Aı ), welche bei starker Verbiegung eintritt, bewirkt 
durch Vermittlung der übrigen Wirbel eine schwache Instabilität. Diese besitzt, wie im folgenden 


Abschnitt erläutert wird, eine gewisse prinzipielle Bedeutung (s. Anh., Abb. 3) Pr 
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(und hinreichend), die Platten soweit zusammenzurücken, daß die schwächste mögliche 
Wirbelkrümmung bereits im Gebiet absoluter Stabilität liegt. Diese neue Stabilitätsbe- 
dingung ist aber darum praktisch bedeutungslos, weil die »letzte« Instabilität (erzeugt 
durch die Schwingung Asonax) bei so hohen Krümmungen liegt, daß ein Plattenabstand, 
welcher alle schwächeren Krümmungen ausschlösse, bereits von der Größenordnung des 
Wirbelkernradius o, wäre. (Nämlich !:d” 10*°, wie aus Abb. 3 folgt). 

Dieses Ergebnis wird aus der exakten Behandlung des Spezialfalls «—= r/2 (wo 
die charakteristische Gl. eine kubische ist) gewonnen und kann durch den Uebergang 
auf allgemeinste Störungen nur noch ergänzt werden, ungeändert aber bleibt die prin- 
zipiell wesentliche Aussage, daß von den größten Plattenabständen bis zu den kleinsten, 
für die unsere Ansätze noch gelten, eine sehr schwache Instabilität auftritt, sobald die 
Wirbel verbogen werden. Da übrigens die Wirbel außer mit der Frequenz Aomax noch 
mit einer (unendlichen) Reihe von kleineren Frequenzen A, parallel zur z-Achse schwin- 
gen, 80 sind die Instabilitäten fast gleichmäßig über die Verbiegungsgrade verteilt, und 
es erfolgt kein plötzliches Aufhören der Stabilität von einem bestimmten 
Plattenabstand an, sondern eine allmähliche Zunahme der Instabilität je 
weiter die Platten auseinanderrücken, wie dies im folgenden Abschnitt näher erläutert 
ist. In allen Fällen aber bleibt die Instabilität von höherer Ordnung klein. 


8. Zusammenfassung und Ergebnis. Ein kurz zusammenfassender Riückblick 
auf den Weg der zur Lösung unseres Stabilitätsproblems führte, und eine kritische Be- 
leuchtung der Resultate ergibt folgendes Bild. 

Zunächst mußte die Schwierigkeit überwunden werden, die darin liegt, daß ein 
dreidimensional gestörter d.h. verbogener Wirbel nicht mehr als unendlich dünn behandelt 
werden darf. Dies geschah mit Hilfe gewisser Ansätze von Lord Kelvin, welche er- 
laubten, die Störungsbewegung eines dreidimensional gestörten isolierten Wirbels in Ab- 
hängigkeit von seiner »Dicke« d.h. seinem Kernradius, anzugeben. Nach Aussonderung 
der beiden einzigen (unter den unendlich vielen möglichen) Störungstypen des Einzel- 
wirbels, welche für die Wirbelstraße von Bedeutung sind, geschah die eigentliche »Koppe- 
lung« dieses Einzelwirbels mit den übrigen Wirbeln der Straße durch die Berechnung 
seines Kernradius aus den (geometrischen) Bestimmungsstücken der Straße. Durch 
diese Berechnung, deren Grundgedanke von Kärmän herrührt, wurde der nach Kelvins 
Methode behandelte »Auf«-Wirbel aus einem beliebigen Einzelwirbel zu einem notwendig 
der Kärmänschen Wirbelstraße angehörigen: seine Störungsbewegung war nun (bei ge- 
gebener Störung»form) durch die Abmessungen der Straße eindeutig festgelegt. Dadurch 
war es schließlich möglich, zu eindeutigen (in Abschn. c) und d) des Anhangs gegebenen) 
Ausdrücken für die durch Zusammenwirken von »Selbstinduktion« und »Fremdinduktion« 
sich ergebenden Störungsfrequenzen eines gestörten Straßenwirbels zu gelangen, welche 
die Stabilitätsbedingung naturgemäß implizieren. 

Es ergab sich, daß die Kärmänsche Sıabilitätsbedingung in allen Fällen ausreicht, 
um die Stabilität auch angesichts der allgemeinsten Störungen zu sichern; abgesehen von 
einer schwachen Instabilität, welche unter allen Umständen auftritt, solange der Abstand 
der Begrenzungsplatten (d) größer ist als etwa !/,oo der Sıraßenbreite (h) oder Breite des 
durch die Flüssigkeit geschleppten Körpers, was natürlich praktisch stets der Fall. Es 
wäre denkbar, daß in dieser schwachen Instabilität ein Teil der Erklärung dafür liegt, 
daß ja in der Tat die Wirbelstraße sich nie weiter als etwa bis zum achten Wirbel ver- 
folgen läßt. 

Eine Abhängigkeit der Stabilität vom Plattenabstand konnte sich — außer der eben 
erwähnten praktisch bedeutungslosen — nicht ergeben. Man könnte aber immerhin eine 
solche in folgendem Umstande erblicken. Durch die Grenzbedingungen an den Platten 
wird, wie wir sahen, unter den Fourier-Komponenten, in die wir eine beliebige Wirbel- 
verbiegung zu zerlegen haben, eine »erste« zulässige Komponente festgesetzt, nämlich 
eben diejenige Sinus-Störung, welche eine volle Wellenlänge zwischen den Platten besitzt. 
Die Krümmung dieser schwächsten zulässigen Verbiegung ist nun naturgemäß umso 
erößer, je kleiner der Plattenabstand; d. h. die scherenden Störungstendenzen in der 
Flüssigkeit, die eben eine Abweichung von der bloßen zweidimensionalen Verlagerung 
des ganzen Wirbels bedingen, müssen umso stärker, die Flüssigkeit umso »unruhiger« 
sein, je näher die Platten aneinander stehen, damit überhaupt eine (»zulässige«) Verbiegunng 
zustandekommen kann. Wenn man nun annehmen darf, daßunter gegebenen Experimentier- 
bedingungen ein bestimmter, stets annähernd gleichbleibender Grad von »Ruhe« der 
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Flüssigkeit erzielt werden kann, so muß man schließen, daß die Wahrscheinlichkeit 
für das Auftreten dreidimensionaler, nicht — Kärmänscher Störungstypen, und somit 
unserer (wenn auch schwachen) Instabilität. um so geringer ist, je kleiner der Plattenabstand. 
Zum Schluß mag noch über die Bedeutung der von uns im Abschn. 5 durchge- 
fübrten Berechnung des Kernradius eines der Straße angehürigen Wirbels einiges gesagt 
werden. Zunächst ist klar, daß die Kärmänsche Wirbelstraße die einzige der be- 
kannten stabilen Konstellationen linearer Wirbel ist, welche diese Berechnung zuläßt, 
denn nur in diesem Falle kann die Arbeit, die bei der Erzeugung des Wirbelsystems zu 
leisten ist, durch physikalisch meßbare Größen (h, U, w) ausgedrückt werden. Ferner 
gilt die ganze Rechnung natürlich nur auf der Basis der Abstraktion: inkompressible, 
reibungslose Flüssigkeit. Und nur insoweit als man zuläßt, daß diese Abstraktion auch 
noeh in Dimensionen anwendbar ist und zu »richtigen: Resultaten führt, wo die Kapillar- 
kräfte bereits wirksam sind, wird man ihrem Ergebnis physikalische »Realität« zu- 
trauen. Umgekehrt aber könnte durch eine experimentelle Feststellung des Kernradius 
in Abhängigkeit von der Breite der Wirbelstraße der Beweis erbracht werden, daß die 
obengenannte Abstraktion ausreichend ist. Freilich sind solche Experimente m. W. bisher 
nicht gemacht worden und wohl auch mit großen Schwierigkeiten verknüpft. 


Anhang. 


a) Darstellung der Koeffizienten des Differentialgleichungssystems (VII): 


Zum Zwecke der numerischen Auswertung seien drei dimensionslose Parameter 
eingeführt: 


k—h!l: 


Innl 


(Maß der Wirbelverbiegung!), wo n infolge der Rand- 
bedingungen beliebige ganze Zahl zwischen 0 und » (siehe GI. III); 


Es seien ferner beachtet die Beziehungen: 
(yıa—Yı,) = 0 
ml 


sowie die Symmetrieverhältnisse bei den Summationen und Intgrationen. 

Dann lassen sich die im Text genannten (s. Abschn. 4) 13 verschiedenen Koeffizienten 
der Taylor-Entwicklung der Gl. (V), eben die Koeffizienten S,... des Gleichungs- 
systems (VII), durch folgende Ausdrücke beschreiben: | 


ro 
0 


0 Po? + =?) ®/a 


1 3%? 
— Pu ( o.3/, 2 =. ö/ ) dz: 


Rs — wie #, bis auf Vorzeichenumkehr; ferner ein unter der ersten Summe 
hinzutretendes Glied: 


+ 2 ( 


1--cos(up,)' 
+ 2°) 
schließlich unter der zweiten Summe den Zähler °/, p.? anstatt 3 k°?. 


Sı = wie zweite Summe von #ı multipliziert mit 
» » Ra ) [— cos cos (k' z)| 
+» 
4 4 sin (up )] eos (k' 3) 
T= — ) dz 


Diese 5 Koeffizienten sind dadurch vor den übrigen ausgezeichnet, daß sie bei 
verschwindender Wirbelverbiegung, also im Grenzfall X —= 0 nicht verschwinden, sondern 
in die bei Kärmän l.c. mit A, B, C bezeichneten übergehen, da sich in diesem Fall die 
Integration in geschlossener Form ausführen läßt (wobei Aı = AR, und $Sı — $8, wird). Alle 
übrigen Koeffizienten dagegen verschwinden im Falle X — 0, so daß hier das Gleichungs- 
system (VII) direkt in das Kärmänsche übergeht. Um nun die Rolle der Kärmän-Bedingung 
auch für kX' 4:0 bequem verfolgen zu können, sei 


Rı (k=0) = R; (k’=0) = Rx und (k’=(0) = S, (k’=0) Sk 
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in unseren Formeln beibehalten und demgemäß geschrieben (vor allem für die numerische 
Darstellung in den Abb. 3 und 4!): Rı= Rx — (Rx — Rı) usw,, wobei also nur jeweils 
das zweite Glied vom Parameter /’ abhängt, während das erste davon unabhängig ist. 

Die restlichen 8 Koeffizienten hängen also sämtlich vom Verbiegungsparameter X’ 
ab (meist im Sinne direkter Proportionalität). Sie sind ähnlich aufgebaut wie die 5 ange- 
sebenen und halten sich größenordnungsmäßig (wie aus Abb. 3 und 4 abzulesen) in der 
Nähe der von k’ abhängigen Bestandteile von ZA und $S, sind also im ganzen Wertebereich 


70 


70° 
70° 


703 


| | 
RAerzz3] 19 70-7 0° m’ 02 793 


Abb. 3. 


Erläuterung zu Abb. 3: 


Auf der Abszisse sind die Werte des Parameters k’' von 10”? bis 10% aufgetragen, 
auf der Ordinate die zugehörigen Werte der zu diskutierenden Koeffizienten. Und 
zwar wird durch: 


Kurve 1 dargestellt: Aomax ' (- ) Kurve 5 dargestellt: —) 


r Anl? (es ist hier das obere oder untere Vorzeichen 
(- ) zu wählen, je nachdem man in Gl, (II) das 
in? obere oder untere Vorzeichen wählt!) 
Kurve 6 dargestellt: (- — 
am? 
» » —(Rx—Rı) - 7), 
» 8 - —). 
Die an der Abszisse in Klammern zugefügten Zahlen geben den Wert von 
k’ 
r == r für n= 1 und gestatten offenbar für einen gegebenen /!/d-Wert direkt die 
an d 


zur ersten erlaubten Verbiegung gehörigen Kurvenwerte abzulesen: alle links von 
dem diesem (/d)-Wert zugeordneten Abszissenpunkt gelegenen Kurventeile sind durch 
die Grenzbedingung verboten, nur die rechts gelegenen erlaubt. (Stets erlaubt 
sind natürlich die Kurvenwerte für X=1017"%, was mit der zweidimensionalen 
Kärmänschen Störung identisch ist.) 

Bei Berechnung der Größen Aymax und A, ist der in Abschn. 5 errechnete 
09-Wert verwendet worden. Dadurch ergibt sich die Wellenlänge der Störung 
längs 2: = 2/k’ + 10°. 00; und dies ergibt als oberste Grenze von X’ für die unsere 
Ansätze ihren Sinn behalten, einen Wert von der Ordnung 10° denn allen unseren 
Betrachtungen lag ja die Voraussetzung zugrunde: 00 — Wellenlänge der Störung 
längs 


| 5 
Z 
/ | | 
| | 
70) 4 + 4 
/ | | d 
107? | | | 
(4 | | | | 
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‚he der Störungsparameter X und « sehr klein im Vergleich zu den jeweils größten der 
ile von k' unabhängigen Koeffizienten und der Selbstinduktionskoeitizienten A, und Aı. Diese 
ist. Eigenschaft ist für uns wesentlich, denn auf ihr beruhen die verschiedenen Vernach- 
K | lässigungen sowohl zur Erzielung der ersten Näherung (s. Abschn. 6 des Textes) als auch 
ge- der zweiten Näherung (s. Anhang unter c) und d)). 
der b) Die Koeffizienten der charakteristischen Gl. (XIV) lauten im einfachsten Sonder- 
ich fall u—=nr/2 (für welchen unter c) die Stabilitätsdiskussion im Einzelnen durchgeführt 
ist), wobei die Größen 83, V,X, P,O,Z verschwinden!): 
0° 7 T 
| | | 


0° 


N 


| 
| | 
| | | > 
| | | \ 
| \ | | | 
| | | | 
| 1 | + 
[RAarz7#| # 8 d 
Abb. 4 
Erläuterung zu Abb. 4: 
Fs stellt dar: 
Kurve 1 die Größe: TI — ) fürk =0|also: —-iTk- (- ) 


» >» » — ) » 10 


“ 
» 2» » : 8, (-- ) für =0| also: Sk'I— ) 
4 


» 3.» » v(- ) » kK=2n:10” 
» 4 » W”=2n+107 


Hinsichtlich der übrigen Koeffizienten der Gl. (VII) genügt es für unsere Zwecke, zu 

beachten, daß: 

So den gleichen Gang zeigt in Abhängigkeit von « wie Kurven 2 und g 
und an keiner Stelle der Größenordnung nach von der Kurve für $j 
abweicht; 
für größere Werte von k’ die gleiehe Kurvenform behält in Abhängigkeit 
von 4, lediglich immer weiter in die Gegend der Kurve 2 rückt, wenn 
k'-—> 10° geht; 

IX ebenfalls mit wachsendem X’ gleiche Kurvenform behält, nur das Maxi- 
mum höher hinaufrückt, so zwar, daß es stets annähernd mit dem 
Kurvenpunkt von für zusammenfällt: 

'Q der Größenordnung nach genau dieselben Kurven gibt wie Y;; 

Pl » » » » » » >» >» |X; 


Y 


Z. höchstens Werte von der Ordnung X max oder IY’Imas erreichen können 
bei gegebenem k'; 
'U' und 'V höchstens Werte von der Ordnung der Kurve 6, Abb. 3, erreichen 
für die ensprechenden K'-Werte; 
IRK—Rı und |R£k— Ra höchstens Werte von der Ordnung der Kurven 7 und 8 
Abb. 3, erreichen können für die entsprechenden k'-Werte. 


I) Es sind hierbei die in Gl. (VITa) Abschnitt 6 eingeführten Abkürzungen benutzt. 
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a = (4°’— 2 — 2 T?+2QY) 
a = + 2(2Kı KkQY-+ Kı 
az; —= A?(Kı Ka — T?)’ — T’Q?Y? 
Im allgemeinsten Fall bei beliebigen «-Werten ist in den drei Koeifizienten 
dy, dx zu ersetzen: KıKs3 durch — L,L.), 
QaY durch PX— 02). 
Ferner treten eine Reihe weiterer Glieder hinzu, die jedoch im ganzen Werte- 


bereich der Störungsparameter sehr viel kleiner bleiben als die oben von uns angeschrie- 
benen und daher im Einzelnen nicht interessieren. 


Die beiden Koeffizienten aı, ax lauten nun: 
+...) 
a = 
(folgen weitere Glieder gleicher Größenordnung). 
Von den Koeffizienten sind reell; as; rein-imaginär: aı, 


c) Stabilitätsdiskussion im Spezialfall «—/2. Hier reduziert sich die G]. (VID) 
auf eine kubische Gleichung in 4°, schreibbar in der Form: 


uU (mit und a a; = As, gesetzt). 
Wurzeln der kubischen Gleichung sind: 


wo € und &” zwei konjugiert komplexe dritte Einheitswurzeln bedeuten, ferner abkürzend 
eingeführt ist: 


yı ebenso, aber mit negativer Quadratwurzel und: 

Die gesuchten Wurzeln A, der charakteristischen Gl. (XIV) für den Spezialfall 

4 — n/2 sind daher: 
und unsere Aufgabe besteht darin, zu ermitteln, unter welchen Bedingungen alle drei 
Ausdrücke (k = 1,2,5) reipimaginär werden. Dazu ist notwendig undhinreichend: 
u, reell und negativ 


und dies impliziert als jedenfalls notwendige Bedingung: 
+ 


reinimaginär. 


Führen wir nun die Abkürzungen ein: 
b=(4KhK T’+24’QY), 
KRT-+...; 


(letzteres unter Vernachlässigung einer Reihe kleinerer Glieder), so schreibt sich der 
Radikand der Quadratwurzel: 


ba? by? 1 | N 272 3 1 2 
was wegen der in Abb. 3 dargestellten Größenordnungs-Verhältnisse näherungsweise über- 
geht in: 
P 4 
(2c—ab) und weiterhin = — Kı 


gültig für alle möglichen Werte des Parameters k' mit Ausnahme zweier kleiner Bereiche 
für welche angenähert @«— 0 wird, das eine Mal, weil bei X’ 10°? (verhältnismäßig 
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schwache Wirbelkrümmung) > 7” wird (7? ist notwendig negativ!) während 
und @ Y dagegen sehr klein, das andere Mal bei X’ = 10*° (sehr starke Krümmung) weil 
dann — , die anderen Glieder wiederum sehr klein. 

Die Kärmänsche Bedingung k = 0,283 aber besagt für uns, daß (A, K,) stets 
negativ, für alle Werte des Parameters %. Führen wir sie ein, so genügen wir also 
unserer notwendigen Stabilitätsbedingung (»reinimaginäre (Juadratwurzel«) überall dort, 
wo der obige Näherungswert für den Radikanden ausreicht. 

Im Ausnahmefall «= 0 dagegen haben wir in Näherung: 


bg3 b3?\ _ 1 

und dieser Ausdruck ist jedenfalls positiv, da c reell und b wegen Ueberwiegens des 
Gliedes 2 Au’Q@ Y und wegen der Ungleichnamigkeit von Q und Y (siehe Erläuterung zu 
Abb. 3) notwendig negativ. In diesem Falle besteht also jedenfalls Instabilität. 


Daß in allen Fällen, in denen > + ; reinimaginär wird, auch die hinrei- 


chende Stabilitätsbedingung erfüllt ist, läßt sich leicht nachweisen. Wir erhalten außer- 
halb der Bereiche 0 und für “= n/2 die Wurzeln der charakteristischen Gl. (XIV) 
in Annäherung als: 


Ist jedoeh «= 0, so tritt in den Wurzeln A; bis A, zu dem imaginären [Bestandteil 
noch ein reeller hinzu, was das Auftreten einer Instabilität bedeutet. Jedoch ist der reelle 
Bestandteil stets sehr klein gegen den imaginären. 


d) Stabilitätsdiskussion für beliebige #-Werte. Es wird ein Näherungsverfahren 
angewandt, welches darauf beruht, daß die Lösungen der allgemeinen Gleichung 6. Grades 
(XIV) sich in erster Näherung durch diejenigen der speziellen kubischen Gleichung 
unter c) approximieren lassen. ls folgt nämlich aus der Größenordnung der Koeffizienten 


des Gliedes dritter und erster Ordnung (a, und «a, in G]. (XIV)), daß man diese in erster 
Näherung gleich Null setzen darf. Unter dieser Annahme löst man die verbleibende 
Gleichung 3. Grades nach dem unter c) dargestellten Verfahren und setzt die so erhaltenen 
Lösungen erster Näherung (}/ genannt) anstelle von A” und A? in Gl. (XIV) ein '). 

Die Wurzeln der so modifizierten kubischen Gleichung erweisen sich als nur um 
Kleines höherer Ordnung von denen der ersten Näherung verschieden und liefern so nach- 
träglich den Beweis für die Berechtigung des angewandten Verfahrens. 

Die ganze so modifizierte Stabilitätsbetrachtung weicht an keiner Stelle wesentlich 
von der unter c) durchgeführten speziellen ab. Wir haben lediglich überall X, X, durch 
(KıKa — LıL;) zu ersetzen (was für X —0 übergeht in Rx?— Sx?). 

Demnach erhalten wir als Näherungswurzeln der charakteristischen Gl. (XIV) für 
beliebige Werte der Parameter « und X (wiederum mit Ausnahme der Bereiche von 
k', wo angenähert «= 0 wird!) die Ausdrücke: 


Diese Wurzeln sind imaginär, falls die Kärmänsche Bedingung erfüllt ist. Für 
die beiden Grenzfälle X = 0 und » gehen sie über in 


bezw. in 4° und 

Die kleine Abweichung von den im Text (Gl. XIlIa und (XIllb) für diese Grenz- 
fälle gegebenen Ausdrücken rührt daher, daß die von uns gewählte Approximation für 
schwache Verbiegungen schlechter ist als für starke. 

Im Ausnahmefall «+0 tritt auch hier zu den imaginären /-Werten ein reelles 
(»instabiles«) Glied hinzu. Ueber dessen Größenordnung und Bedeutung ist genau das 
unter c) bezüglich des dort auftretenden instabilen Gliedes Gesagte zu wiederholen. 


I, Dabei ist ein Punkt zu beachten. Bezeichnen wir das Lösungssystem erster Näherunge durch 
hy’ bis Ag’, so erhalten wir jede der Lösungen zweiter Näherung A, bis /; dureh Einsetzen jeweils der 
mit der gleichen Nummer versehenen gestrichenen Größe. Dies ist in unserem Fall darum von 
Wichtigkeit, weil in einem sewissen k-Bereich die Wurzeln }ıa von anderer (Größenordnung sind als 
und Ase. 
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Anmerkung: Man sieht übrigens leicht, daß das in A, und A, eingehende Ver- 
hältnis 00 :! nicht gerade den im 5. Abschnitt errechneten und in Abb. 3 benutzten Wert 
haben muß, damit unsere obigen Stabilitätsbetrachtungen und deren Ergebnisse gültig 
bleiben. Im Gegenteil ist ihre Empfindlichkeit gegen Aenderungen von @/l sehr gering: 
eine Multiplikation des in Abb. 3 verwendeten Wertes mit einer Zehnerpotenz würde sich 
noch an keiner Stelle ändernd bemerkbar machen. Diese Feststellung ist von Bedeutung 
für den Fall, daß sich etwa experimentell ein anderer Wert für @/l! bzw. o,/h ergeben 
würde als der theoretisch ermittelte (siehe auch Abschnitt 8, Schlußabsatz). 872 


Das Förderhöhenverhältnis radialer Kreiselpumpen 
mit logarithmisch-spiraligen Schaufeln.' 


Von A. BUSEMANN in Göttingen. 


(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


ie gleichnamige Arbeit von W. Schulz stellt den Versuch dar, zur Lösung des 

gesamten Problems mit einer vereinfachten Abbildungsfunktion vorzudringen und 

die hierdurch benötigten Berichtigungsfunktionen durch Anschluß an die bekannte 
Förderhöhe bei unendlicher Schauielzahl aufzufinden. Trotzdem bleibt dies nur eine 
Näherungslösung, die in rechnerisch schwer zu erfassender Weise von der wirklichen 
Lösung abweicht. Unabhängig davon ist die Arbeit jedoch vorbildlich in der einfachen 
llerausarbeitung der verlustlosen Kennlinie. Außerdem sind die Ergebnisse wirklich alle 
durchgerechnet und mitgeteilt. Es ist daher anzunehmen, daß vielfach der Wunsch be- 
steht, doch auch die genaue Lösung wenigstens als Vergleich für einige Werte zu kennen. 
Aus diesem Wunsch ist die vorliegende Arbeit entstanden und es erscheint daher zweck- 
mäßig, sie in engem Anschluß an den Schulzschen Text abzufassen. Sie soll jedoch 
auch für sich allein verständlich bleiben, zumal sich zeigen wird, daß sich die genaue 
Lösung in sehr einfacher und bandlicher Form darstellen läßt. 


1. Rechnungsgrundlagen und Voraussetzungen. Die Voraussetzungen 
zweidimensionaler idealer Flüssigkeitsströmung, die sich mit den Methoden der konformen 
Abbildung behandeln läßt, sind dieselben wie bei Schulz, ebenso die Bezeichnungen: 
Bezeichnungen: 

z — Liaufrad-Schaufelzahl, 
P = 7/2 — y = Laufradschauielwinkel, 
r — Laaufradschaufelradien, 
Indizes 1, 2 für Ein- bzw. Austritt, Index s für die Schaufelkontur, 
«@ — 7,2 — @— Schaufelwinkel eines etwaigen Leitrades unendlicher Schaufelzahl ’), 
theoretische Förderhöhe bei endlicher und unvendlicher Schaufelzahl, 
Y = sekundliche Fördermenge, 
Winkelgeschwindigkeit, 
» — u + ir — Koordinaten des Gitters, 
Koordinaten der Schaufelradebene im Maßstab n, = 1, 
ww — u 42er Koordinaten bei Abbildung der Schaufeln auf den Einheitskreis, 
— Längen und Winkel in der -Ebene (Abb. 1), 
S = V/z = (uellstärke je Schaufel, 
J = Zirkulation um jede Schaufel, 
— Umfangsgeschwindigkeit auf dem Kreis r;, 
— mittlere radiale Durchflußgeschwindigkeit auf dem Kreis 
h, /to, /t,, — Koordinaten und Bestimmungsstücke der Kennlinien (Abb. 2), 
h., ec. = Koordinaten des Punktes stoßfreien Eintritts (Abb. 2). 


Abkürzungen: Mathematische Bezeichnungen: 
a—1lw, a=1/w a —= konjugiert komplex zu a, 
vergl. Gl. (12), N (a) = (a+a)—reeller Bestandteil von a, 
S,8: (17), (a) (a — a) = imaginärer Bestandteil 
N (24), von a, 
P.i (30), | \,., 
‚a == !a»-a= absoluter Betrag von a. 


I) Ergänzung zum gleichnamigen Aufsatz von W. Schulz, diese Zeitschr. Bd. 8 (1928), S. 10. 


?), Da das Leitrad sehr klein anzunehmen ist (punktförmig), kann «’ auch als effektiver Schaufel- 
winkel bei endlicher Leitschaufelzahl anzesehen werden. 


DW. Band 8, Heft 5 
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er- 2. Vorbereitende Untersuchungen. Die z Schaufeln werden durch Potenzieren 
ert mit z auf eine ebenfalls logarithmiseh spiralige Schaufel gleichen Schaufelwinkels aber 
ug verkleinertem Innenradius abgebildet. Diese wird abweichend von Schulz mit der be- 
5: kannten Königschen!) Methode auf eine Gerade abgebildet. Es wird jedoch zum besseren 
ich Vergleich mit Schulz die Königsche Gerade gleich auf einen Kreis abgebildet und 
ng nur das Gitter zur Orientierung mit angegeben. Die Abbildungsfanktion einer logarith- 
en 


misch spiraligen Schaufel auf den Kreis entspringt aus dem Vergleich der Wirbelquell- 
strömung aus dem Mittelpunkt der [-Ebene, deren logarithmisch- spiralige Stromlinien die 
Radien unter dem Winkel y schneiden und somit die Schaufel störungsfrei umströmen, 
mit der Strömung einer ebensolehen Wirbelquelle im Punkt wo, der w-Ebene, die den 
Einheitskreis als Hindernis umströmt. Die Singularität der wirklichen Strömung im 
Mittelpunkt der Ü-Ebene wird daher nur beim Potenzieren entsprechend der Schaufelzahl 
für jedes Riemannsche Blatt verringert und bleibt sonst beim Uebergang auf die 
w-Ebene dieselbe. 


du | 3. Die Abbildungsfunktionen. Die vorher erwähnten Abbildungsfunktionen 
and | lauten (Abb. 1): | | 

ine wo — 1 wo — | 

hen 


Da der Leser die Richtigkeit gegebener Abbildungsfunktion jederzeit leicht prüfen kann, 


hen werden hier gleich weitere Beziehungen angegeben, die auch bei einer solchen Prüfung 
alle von Vorteil sind. 

be- 

Ien. IV 

:ck- 


E= sin 


T, 
1°), Abb. 1. Abbildungsebenen. 
‚ahl, 
4. Abbildungskonstanten, $Schaufelgleichung. Durch Diiferentiieren von (1) 
ergibt sich für die Linienelemente folgendes: 
> wo — I/w w 
Die Gleichung der Schaufelkontur ist die Gleichung des Einheitskreises der w-Ebene. 
Diese lautet: ws — e'® 
daraus: 
wo— erd\l/z (wo — t wo — ei ® wo — 
— 1 — 1 2 wo l uw — 1 
Die äußere Schaufelkante hat die Koordinaten: 
G=1, 02.0.4. 
Die innere Schaufelkante liege bei d, = — 28. Dann ergibt sich: 
wo + e2ie\l/z (wo + 
‚dteil wo — 1 (5) 
ER wg + wn + e2i: t | 
- )+etin( ) = 2-cosyIn‘“' —4siny-e | 
a. wo —1 | 
3. 10. I) E. Könie, Potentialströmung durch Gitter, diese Zeitschr. Bd. 2 (1922), 8. 122. Vergl. 
aufel- auch E. Sörensen, Potentialströmungen dureh rotierende Kreiselräder, diese Zeitschr, Bd. 7 (1927), S. 89 
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Daraus folgt das Radienverbältnis vr, zu rz: 


_ . 
0 2 


01 Wo) E cos Yy, Wo *C0OSE—+ 


und 0 sind keine Vektoren). 


Gl. (6) ermöglicht bei gegebenem r,, 7/2? — y und z den Punkt der w-Ebene 
zu bestimmen!). Längs der Schaufelkontur gelten außerdem die vereinfachten Beziehungen: 


=) 
er 


- 6 wg — em! 
= ( ) . | 
1 . 


5. Das komplexe S$Strömungspotential. Das Strömungspotential werde als 
Superposition der drei Teile a) Durebflußströmung bei ruhendem Schaufelstern, b) Ver- 
drängungsströmung des rotierenden Schaufelsterns allein, e) Zirkulation um die Schaufeln 
zur Ermöglichung glatten Abflusses behandelt. Nachdem die Abbildungsfunktionen fest- 
liegen, läßt sich die Strömung in jeder Ebene weiter behandeln. Für die Beachtung der 
Grenzbedingungen bei a) und c) eignet sich am besten die w-Ebene, wo die Schaufel 


ein Kreis ist, nur die Verdrängungsströmung des rotierenden Schaufelsterns hat in der 
s-Ebene einfachere Grenzbedingungen. 


a) Reine Durchflußströmung. Da die Beschränkung auf drallfreien Wasser- 
eintritt keinerlei Vereinfachung der Rechnung bringt, sei hier abweichend von Schulz 
der allgemeine Fall behandelt. Hierfür lautet das Potential, wenn die Quellstärke 5 und 
der Winkel der Stromlinien der Wirbelquelle gegenüber den Radien « ist: 


(1 — ite «)1In (w — 


oder mit Umströmung des Einheitskreises: 


Za —itze)In(e — + (1-+3% tg @)1In ( 
zn w 
Die Geschwindigkeiten der Durchflußströmung in der w-Ebene sind: 
1 \ 
(1 —itga) (1—itg«) / 
Ss +itene s +itga)w w 


Für die Schaufelkanten lauten diese Geschwindigkeiten unter Berücksichtigung von (4) 
und (5) und der Beziehungen: 


Si2cosy 


analog‘) 


ay 


2nwı 


Der Fall drallfreien Eintritts (Schulz) ist durch «=0 gegeben. Für den Fall «=}Y 
bilden die Schaufeln in der Ü-Ebene kein Hindernis, in der w-Ebene werden dann 
und 0, Staupunkte (dies ist die »Abbildungsströmung«, vergl. Ziff. 2). 


b) Verdrängungsströmung. Die Rotation des Schaufelsterns wird durch An- 
bringung von @uellen auf der Vorderseite und Senken auf der Hinterseite der Schaufel 


I) Siehe Ziff. 8, G1. (32) und (32a). 
?) Die Analogie ist vollständig bei Verzicht auf die unsymmetrischen Festsetzunzen: wa = 1 und 
+. Das verlangt "Car Statt cag und statt y in (9. 


. N . . . (7). 


cew, 
lech. 


als 
Ver- 
ıfeln 
fest- 
der 
rufel 
der 


sser- 
ıulz 
und 


n (4) 


‚aufel 


1 und 
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dargestellt, deren Ergiebigkeit so zu wählen ist, daß die Geschwindigkeiten an der 
Schaufelkontur normal zur Schaufeloberfläche gleich der Normalkomponente der Geschwin- 
digkeiten bei der Rotation der Schaufel werden. Das Potential dieser sorenannten Ver- 
drängungsströmung kann in der w-Ebene durch eine Laurent-Entwicklung') angesetzt 
werden: 
w w“ k—=1 W 

die mit komplexen a, für das Aeußere des Kreises die allgemeinste Quell-Senken-Vertei- 
lung auf der Schaufelkontur darstellt. Zur Bestimmung der a, dient die Umfangsgeschwin- 
digkeit der rotierenden Schaufel (die wegen der Annahme r, — I aus Dimensionsgründen 
den Faktor hat): 


und die Strömungsgeschwindigkeit in der [-Ebene: 
dw 
w 
Dreht man ein Schaufelelement so, daß seine Normale in die Richtung der reellen Achse 


zeigt, so muß für dieses Element die Relativgeschwindigkeit (ex — c,) rein imaginär sein wegen 
der Undurchdringbarkeit der Schaufel. Daraus ijolgt für die ganze Schaufelkontur: 


Hieraus bestimmt man die a bzw a. Durch Einsetzen der Werte aus Gleichung (5) 
und (7) erkennt man, daß in dem 2. Faktor in der Klammer auf der linken Seite von 
(11a) ein Ausdruck zusammengefaßt ist, der rein reell ist und daher auf die reehte Seite 
gebracht werden kann. Der Rest auf der linken Seite wird durch Diiierentiieren von 
Gleichung (10) gewonnen (er stellt die Normalkomponente der (reschwindigkeit an der 


Schaufelkontur in der w-Ebene dar, die durch Multiplikation mit ©» reell gemacht ist und 
deren Größe durch die Rotation in der {-Ebene vorgeschrieben ist). Dann erhält man: 


e'Y 
wo — e! wo 
wobei die Abkürzung eingeführt ist: 
1 + 1+c0os2y isin?2 1 + e2'7 | 


Q 


u 


Durch Multiplikation mit e* '*®. dö/27 und Integrieren zwischen den Grenzen —z und +4 7 


lassen sich die a, bzw. a, links trennen: 
+r 


== — 2 .o-ik6.d Ö (13). 
ZT. 8 ( 1 ) 
uo — 6 w— ei 7 


Durch (13) sind die a; bestimmt. Die Summierung der a, läßt sich jedoch auch vor 
Auswertung des Integrals durchführen. Die gesuchte Geschwindigkeit c,;, an der Aus- 
trittskante der Schaufel lautet in der w-Ebene wegen un = 1: 


+r 

i’ cos 0 1 \ 
k—0 1 z ( 1 ( k—( 
wo — e! w— 6, 
') W.Spannhake, Konforme Abbildung; Hydraulische Probleme, herausgeg. vom Wiss. Beirat 


des VDI, Berlin 1926. Vergl. auch Anm. 1, S. 376. 
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0, ist rein imaginär wegen der Grenzbedingung (11). Und so kann man zur Erhöhung 
der Symmetrie schreiben: 


2 2 (wo — 1)": ywu'e wor el® 1 
(ei . . (15). 


Ersetzt man die Summen im letzten Faktor, die geometrische Reihen darstellen, durch 
ihren bekannten Summenwert 


so erhält man ein Integral '), das sich graphisch auswerten läßt und dem von Schulz aus- 
gewerteten entspricht und sogar für = 90°, y — 0° identisch mit seinem ist. 
Man geht jedoch zweckmäßig den umgekehrten Weg, daß man alle anderen Fak- 


toren des Integrals in Summen von Potenzen von e'? verwandelt. Dann erhält man mit 
l/wo =a als Abkürzung: 


4 
wo +eiY-al-am 
Da das Integral zwischen — z und + nur für die nullte Potenz von e’? einen Beitrag 


liefert, sind alle positiven Werte & zu ermitteln und durch /, m und n auszudrücken, für 
die der Exponent Null wird. Die Summe läßt sich dann vor das Integral ziehen und 
das Integral liefert nur den Faktor 27. Sieht man die Summe von endlichen und un- 
endlichen geometrischen Reihen als bekannt an und berücksichtigt 

1 


—a uo—]1 02 


so hat man nur noch zu beachten, daß die Binomialkoeffizienten folgende Beziehung erfüllen ?) 


um die Gleichung (15) auf folgende Gestalt zu bringen: 


(1 — a) (1 — a) 92 


(16). 


Das Integral aus Gleichung (15) ist hierdurch in einen Ausdruck mit zwei Summen verwandelt: 


') In dieser Form ist das Integral von der Laurent-Entwicklung unabhängig. Sie diente nur 
zur Ableitung des Satzes, daß bei einer Quellsenkenverteilung auf einem Kreis die Normalgeschwindig- 
keit am Punkte Ö des Kreises gleich der halben mittleren Quellstärke plus der halben Quellstiärke am 
Punkte Ö ist. Da hier die mittlere Quellstärke Null ist, stellt die linke Seite von (11b) die Normal- 
geschwindigkeit und zugleich die halbe Quellstärke längs des Kreises dar, was durch (13) ausgedrückt 
ist. (14) berechnet bei zegebener Quellstärke die Geschwindigkeit am Punkte wg = 1 als Vektorintegral, 
(15) als Komponentenintegral bei bekannter Richtung. Hierbei tritt der reziproke Wert der Entfernung 
— des Punktes von wa =1 als Summenwert der geometrischen Reihe auf. 


=», 
(1 a) N > -1 
l—a Mm m 


= (1 - a) ai. 


| 
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am 
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S;ı und 8, sind reell und lassen sich reell auswerten, wenn man stets die konjugiert kom- 
plexen Faktoren beieinander läßt'!),,. Außerdem sind sie miteinander verwandt, so daß 
die meiste Rechenarbeit für beide gilt. Sı konvergiert stets rasch, $, etwas weniger rasch 


für a-a nahe 1. Das Verhältnis S, zu Sı wird bei konstantem q für a-a— 1 am größten. 


Für diesen Fall (a yore l, rı/ra = 0) läßt sich die Summe $, unter Benutzung der Gamma- 
funktion schreiben: 


( 


. .(17a), 
I (gQ)- II (Q) 
die für reelle Werte g tabellarisiert ist*). ‚Sy bekommt dann den Wert 


Dies sind zugleich die größten Werte, die Sı ung 8, sowie S,/Sı erreichen können. 
Gleichung (16) läßt sich auch in den übrigen Faktoren reell schreiben unter Benutzung 


der Abb. 1 und der Gleichungen (6), (12) und (17) sowie (a-a)? — a. 


2 08 cos — cosy 
Sı + | "2.83 
ı 727 + cosE + a; cosy (16 ) 
” > 1+cos2y 2:-sin?2y 
analog’) ergibt sich c,.: 
2° a co8y COSE + a Cosy 
COSE— a COSE— a 
cos2y 2e'sin2y 
2° 
+ a cos y) 


c) Schaufelzirkulation. Der Grenzfall zäher Flüssigkeit mit verschwindender 
Zähigkeit unterscheidet sich von der idealen Flüssigkeit nur dadurch, daß an der Schaufel- 
austrittskante nur dann eine stabile Strömung besteht, wenn ein glatter Abfluß vorhanden 
und keine Umströmung der Kante nötig ist. Die bei der idealen Flüssigkeit unbestimmte 
Schaufelzirkulation erhält durch diese Bedingung einen festen Wert aus der Gleichung: 


Die Schaufelzirkulation hat in der w-Ebene das Potential 


Z=' "lkw. 


Hieraus berechnet sich die Geschwindigkeit für die w-Ebene 
d Z- J 
Ce = 
dw 2nw 
und für die Kanten: 


wı 


d) Resultierende Strömung. In Gl. (17) ist die Bedingung für glatten Abfluß 
angegeben, aus der sich J bestimmt 


2 
J= 
Die Umströmung der Eintrittskante hat dann die Grschwindigkeit 
a3 + 
wy 


Diese Geschwindigkeit bedeutet in der [-Ebene eine theoretisch unendliche Geschwindig- 
keit an der Eintrittskante. Die hierauf zurückzuführenden Verluste bezeichnet man als 


I) Siehe Zift. 8. 
2) Jahnke und Emde, Funktionentafeln, Leipzig 1909. 


3, vergl. Gl. (9). Jedoch sind $ı und Sy nicht zu vertauschen und ist rı? durch rg? nach Gl (6) 
ersetzt. = 9=— 


| 
j 


h. 
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»Stoßverluste«. »Stoßfreier« Eintritt ist erfüllt für c, = 0, und ergibt mit (9), (9a), (16a) 
und (16b) 


Ca? 1 


w; wı ei 02 
2 cosy 2 a 
4 + COS — ja) C0S + COS Y 
wy (1 — a)! -(1 — a)! 
Nach Gl. (6) ist aber: 
1 1 2 al cosy 2 la cCOSE+ a C0Sy COS & — C08 y 
cOSE— a cosy cosE + |a c08 cos ! a cos y cosE+ a 
Die Bedingung für stoßfreien Eintritt lautet daher: 
/ — 2.89 
(1 —a)!-(1 — 


6. Das Förderhöhenverhältnis. Für die Förderhöhe der Kreiselpumpe ist die 
Schaufelzirkulation ausschlaggebend. Dem Aultriebssatz geradlinig bewegter Körper von 
Kutta-Joukowsky entspricht der Momentensatz bei rotierenden Körpern: 

wobei z:J die Zirkulation um alle Schaufeln und 2:8 = V die Quellstärke aus dem 
Zentrum angibt. Die theoretische Förderhöhe /,, ist die Arbeit am Schaufelrad bezogen 
auf die Durchflußmenge V. Sie ergibt sich aus der Arbeitsgleichung: 


Mm — 


— 
2: J-a 
in der Form: Hh=——, e oder, da J aus (19) bekannt ist: 


Nach den Gleichungen (16a) und (9) ist bei gegebener Laufradschaufelung «3 


proportional ® und proportional —tge). Gl (21) ist daher bei gegebener 
Schaufel darstellbar in der Form: 


9 \ 
+h,: (te — tg > 4 (21a). 


r„»@ ist die Umfangsgeschwindigkeit «„ der rotierenden Schaufel an der Austrittskante. 


m. 


ist die mittlere radiale Geschwindigkeitskomponente ce, der Durchflußströmung auf 


dem Kreis 73. Auf diese Weise wird die Gleichung für F,. noch anschaulicher: 
. .» . .  (21b). 
dabei bedeuten die Größen %, und A.: 
19°: Ss-(tgy—tga) 09 a +cosy 


Die Eulersche Förderhöhe einer Pumpe, die für unendliche Schaufelzahl bei end- 
licher Länge gilt, hat, wie aus der elementaren Pumpentheorie bekannt ist, den Wert: 
= Ur‘ — U Wı 
wobei hier u, und u, die Umfangsgeschwindigkeiten, w, und wz die absoluten Ein- und 
Austrittsgeschwindigkeiten und «@ und P' die Winkel zwischen u, und wı sowie us und w; 
sind. Man findet daraus für den zweidimensionalen Fall (konstante Schaufelhöhe) mit 
den neuen Bezeichnungen: 
Es nehmen also beide Größen A, und A, im Eulerschen Fall den Wert 1 an. 
Durch Division mit c„? läßt sich (21b) dimensionslos schreiben: 
Hrn (r 
h= — A [> (tg «)) =M + h,-c, (210). 


Cu 


Durch Gl. (21e) ist die verlustlose Kennlinie der Pumpe in den dimensionslosen Koordi- 
naten c, und A gegeber, und zwar ist sie stets eine Gerade (Abb. 2). Jede Kennlinie 
ist daher durch zwei Bestimmungsstücke festgelegt. 
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7. Numerische Aus- 
weriung. Die numerische 
Auswertung soll sich auf die 
Bestimmungsstücke der ver 
lustlosen Kennlinie und auf i 
den Punkt stoßfreien Kin- u |, y. 
. r . . 
tritts auf dieser Kennlinie 
erstrecken. Um anschauliche 
Größen darzustellen, sind für 
die Auswertung folgende Be- VAN 
stiimmungsstücke in den di- | 
mensionslosen Koordinaten 4 
gewählt: die Förderhöhe ‚Be 
bei der Menge Null, c„ die | 
der Förderhöhe Null, sowie Y , usw‘ | 
die Fördermenge und A. 
die Förderhöhe bei stoß- Drehsinn (466.7) 
meinsame Nenner aller Be- 
stimmungsstücke soll durch | 
den Buchstaben N abgekürzt Abb. 2. Verlustlose Kennlinie. 
werden: 
4cos’7 2...sin?Y 
N=(1—-a):(1 = -e : 
Dann haben die angegebenen Bestimmungsstücke folgende Werte: 
2 a *cosy E — *C08Y 
. + - — Rp) 
N 
h) 2 ja| cosy ar 
Cn = (2 
hy N 
und für den Pankt stoßfreien Eintritts durch Einführung der neuen Bezeichnungen in (20a): 
—2'8 2.89 
N (26), (27) 


In den Abb. 5 bis 10 und 11 bis 16 sind die Größen /o, —c„ und A, für ver- 
schiedene Schaufelzahlen und Schaufelwinkel über den Radienverhältnissen aufgetragen. 
Auf die Darstellung von c, wurde verzichtet, weil aus der Proportion «,:c„— (Ro — h.):h, 
leicht zu finden ist. Die Werte c, und r, sind stets negativ in Abb. 2, d.h. man erreicht 
sie beim Drehsinn — ®» in Abb. 1. . 

Für die einzelnen Größen ergeben sich folgende Grenzen: 


= | für rn = 1 
h. =0 für rı/rn = 1 
! h. = Mo für rılra ze 0, 


Man erkennt daraus, daß die Kennlinien für positive Werte h stets unter der Eulerschen 
Geraden liegen, soweit sie nicht in die Eulersche Gerade selbst übergehen. Der Schnitt- 
punkt der Kennlinien mit der Eulerschen Geraden liegt also stets bei negativem h, wie 
dies auch schon in der Schulzschen Arbeit festgestellt ist. Es kommen auch tatsächlich 
alle negativen Werte A als Schnittpunkte vor, eingeschlossen die Grenzen Null (für rı/r3 — 1) 
und —» (für rı/ra = 0). 

Der fast unstetige Anstieg der Aho-Kurven bei kleinem % und kleinem z erklärt sich 
daraus, daß die Schaufeln die Strömung erst dann gut führen, wenn die Schaufellänge 
einen bestimmten Mindestwert überschritten hat (siehe Abb. 3 für f = 5"). 
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8. Anweisung für die Berech- 
nung eines Beispiels. Gegeben ist in 
praktischen Fällen rı, Den 
Leitschaufelwinkel « braucht man, sofern 
überhaupt ein Leitapparat vorgesehen ist, 
erst bei der Uebersetzung des dimensions- 
losen Diagramms in die Wirklichkeit. 
Hierin liegt der Vorteil, der durch Ein- 
führung der Größe c, als Abszisse er- 
zielt ist. Aus den beiden Bestimmungs- 
stücken - + (>) und findet 

sın 

man aus Abb. 4 die Entfernung «, und 
unter Berücksichtigung des Maßstabes 
auch die Zahl wu (in Abb. 4 ist der 
Radius des Einheitskreises gleich I cm 
gewählt), Der reziproke Wert von wo 
Abb. 3. Enge Spirale. ist a. Der spitze Winkel & berechnet 

sich aus: sine= a cosß. Mit dem ge- 
fundenen a berechnet man die Summen Sı und 8, (Gl. 17) unter Benutzung von 


a-a= a’, y=90’—P und: 


1? 2? 14 
Unter Verwendung des Nenners Gl. (24) ergeben sich die Größen ho, /., Cm und c, aus 
den Gleichungen (22a), (25), (26), (27) durch einiache Rechenoperationen. 

Diese Anweisung versagt, wenn die Größe a so nahe an 1 ist, daß der zugehörige 
Punkt z, aus Abb. 4 nicht mehr gut zu entnehmen ist. Dies ist immer der Fall, wenn 
out führende Schaufeln gewählt sind, wie das in der Praxis meist geschieht. Man kommt 
daher in allen praktischen Fällen mit den Entwicklungen der Bestimmungsstücke in der 


Nähe von a — I aus. Diese lauten mit der Nebenbedingung, daß IE 1, die be- 


+0° 20° 70° 0° 
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Abb. 4. Lage der Punkte wo = 1/a, 
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sart, daß die gegenüberliegenden Führungswände der Kanäle zwischen zwei Schaufeln 
keinen stumpfen Winkel miteinander bilden sollen: 


S 1 
ho = „= (ho) == ho max ( .) = 1 — In N 
7) l N ali==1 

2 

Re = Roman — ( ) 

rg 


wobei P folgendes Produkt bedeutet: 


4 cos? y 4 cos? y 4 cos? y 
4 y 4 cos? I=1 ] 4 y 
z / 
rı\2 
G.=- ho max + ( ) N P . 1). 


Diese Weite Aomax Sind in Abb. 17 dargestellt und lassen sich auch für andere Schaufel- 
zahlen interpolieren. Man braucht zu ihrer Berechnung nur die schnell konvergierende 


Summe Die Ableitung der Summe nach wird unendlich für in 


der Nähe von a = 1. Untersucht man den Charakter der Summe näher, wobei der 
Uebergang in die Eulersche Gleichung für A, bei z=» einen Fingerzeig gibt, so findet 
man das angegebene Produkt P, das sich als schnell konvergierend erweist. Sollten für 
einen Fall die Näherungen nicht genügen, so kann man stets aus Gl. (6) in der Form: 


r2 


a’— 


2 
( eos co8° y + _ ( Y.e2': sin? y 


oder für a nahe 1, &e nahe y: 


den Wert a’ genau bestimmen, indem man € aus Abb. 4 entnimmt und sukzessive ver- 
bessert aus sine— a -siny, und nun wirklich die Gleichungen (22a), (25), (26), (27) be- 
rechnen. 


9. Vergleich mit bekannten Arbeiten. Die Gedankengänge dieser Arbeit 
sind nicht neu, sondern gehen auf die Kucharskische‘) und auf die Spannhakesche’) 
Arbeit zurück, in denen der Fall # = 90° durchgerechnet ist. Die Uebereinstimmung der 
Ergebnisse ist vollkommen. Der Uebergang in die elementare Eulersche Pumpentheorie 
bei unendlich guter Führung (2=” oder 3 = 0°, letzteres wegen der guten Führung der 
enggewickelten Schaufel) ist erfüllt. Die Werte der Arbeit von König*) für n/rn = | 
(gerades Schaufelgitter) werden ebenfalls erreicht. Das ganze Gebiet umfassen nur die 
Arbeiten von Sörensen’) und Schulz‘). Schulz arbeitet in seiner Arbeit die verlust- 
lose Kennlinie heraus, allerdings mit einer Näherungsmethode, die die 3 Parameter r,/r3, 
z in 2 herabdrückt. Die Uebereinstimmung ist für 90° und = 0° exakt, bei 
nahe 90° gut, jedoch bei ungefähr %?= 20° weniger gut (siehe Abb. 17 und 18). Die 


q q q q 


4Jcos? 


al=1 Y + 


?, W. Kucharski, Strömungen einer reibungsfreien Flüssigkeit, München 1918, 
3) W. Spannhake, a.a.O0. 

E. König, a.a. 0. 

E. Sörensen, a.a.0. 

6) W, Schulz, a.a.O. 
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Arbeit von Sörensen unterscheidet sich von der vorliegenden, soweit sie den Punkt 
stoßfreien Eintritts betrifft, den Sörensen untersucht, nur dadurch, daß Sörensen statt 
der Summen Integrale auszuwerten hat. Dadurch blieb die Rechenarbeit noch zu groß, 


70 /3:90° x 70 3-60 oo 
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Abb. 5 bis 10. Ao-Kurven. 
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Abnh. 11 bis 16. —cm- und Ahs- Kurven. 
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um daß ganze Gebiet in Kurvenform darstellen zu können. Das von ihm angegebene 
Beispiel stimmt folgendermaßen überein: 


Es handelt sich um eine Pumpe mit z= 6, rı/r, = 0,455, 7 = 60°, daraus 1 — a’ 


—= 0,00294 nach Gl. (32a). Das Verhältnis der Drebzahlen für stoßfreien Eintritt bei 
endlicher und unendlicher Schaufelzahl ist 


| 
C T2 1 


— = = — 0,703 (bei Sörensen 0,715). 


Das Verhältnis der Drehmomente bei der Drehzahl entsprechend :v,' ist dann: 


2 
h, N Sı + 2° rı 2 
( N a\=]| 


72 


Steigert man die Schaufelzahl bei konstantem a, so erhält man 


lim 4 cos? y -log(1— a?) 
Do 


2:0 


10 10 
Pro 
08 08 
06 06 
60° 30° 0° 02 06 08 70 
\RA 8987 16 
Abb. 17. Ahomax-Kurven. 
Abb. 18. M= -Kurven 
Hrn ==() 


nach Schulz für = 20°. 


Die vorliegende Arbeit liefert bedeutende rechnerische Vereinfachungen gegenüber 
den bisherigen, soweit es sich um die äußeren Daten einer radialen zweidimensionalen 
Kreiselpumpe mit log.-spiraligen Schaufeln handelt. (Für die inneren Daten, Geschwin- 
digkeit- und Druckverteilungen längs der Schaufeln kommt eine Summe hinzu, die schwerer 
auszuwerten ist.) Damit erıwuchs dem Verfasser die Pflicht, das ganze Gebiet in Kurven- 
form darzustellen, was in den Abbildungen in groben Zügen geschehen ist. Zugleich 
muß aber gewarnt werden, den Mangel an theoretischen Untersuchungen bei anderen 
Schaufelformen damit beheben zu wollen, daß die Ergebnisse, insbesondere soweit 
sie sich auf den stoßfreien Eintritt beziehen, ohne weiteres verallgemeinert werden. 
Es sei nur noch bemerkt, daß die bei Pumpen nicht übliche Strömung von außen 
nach innen durch den Schaufelkranz bereits in den Formeln enthalten ist, und zwar 
für 7 > 90°. 898 
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Mathematische Theorie der Wärmespeicher. 
Von WERNER SCHMEIDLER in Breslau. 


inleitung. Die mathematische Theorie der Wärmespeicher, wie sie im Hüttenbetriebe 

bei den Hochöfen bzw. bei Siemens-Martinöfen zur Vorwärmung des Kaltwindes in 

Gebrauch sind, ist trots wertvoller Vorarbeiten!) noch nicht soweit entwickelt worden, 
daß sie für die Praxis direkt hätte nutzbar gemacht werden können. Da mir andererseits 
ein lebhaftes Interesse an einer solchen Entwicklung der Theorie von berufener Seite 
geäußert worden ist’), so soll es die Aufgabe der vorliegenden Zeilen sein, diese Lücke, 
soweit möglich, auszufüllen. 

Die strenge Durchführung dieser Aufgabe führt auf ein System von Volterraschen 
Integralgleichungen, deren direkte Lösung infolge der in den Lösungen der Natur der 
Sache nach auftretenden Unstetigkeiten mit erheblichen Weiterungen verbunden ist. Um 
diesen zu entgehen und trotzdem zu einem für praktische Zwecke brauchbaren Gleichungs- 
system zu gelangen, beschränke ich die Fragestellung auf einen zeitlichen Mittelwert der 
in Frage stehenden Temperaturen, der in erster Linie praktische Bedeutung hat, und 
gelange auf diese Weise mit Hilfe gewisser Approximationen von den Integralgleichungen 
zu einem System von gewöhnlichen Diifferentialgleichungen, die linear sind und deren 
Koeffizienten überdies nur ein wenig von konstanten Werten abweichen. Da auch die 
ursprünglich entstehenden Differentialgleichungen etwa mit der Runge-Kuttaschen Re- 
gel leicht numerisch behandelt werden können, dürfte der dem Problem entsprechende 
Grad der Einfachheit erreicht sein. Für grobe Annäherungen genügt die Annahme von 
konstanten Koeffizienten und damit die Approximation durch Exponentialfunktionen. 


1. Formulierung des Problems. Die Wärmespeicher, die hier behandelt wer- 
den sollen, beruhen in ihrer Wirkungsweise darauf, daß ein zylindrisch ’) gestaltetes 
»Gitterwerk« aus feuerfestem Steinmaterial zuerst in der sogenannten Heizperiode einem 
heißen Gasstrom ansgesetzt wird, der durch Verbrennung des Gases in dem »Brennschacht« 
erzeugt, dann in der Richtung der Zylinderachse durch das Gitterwerk geleitet wird und 
dort seine Wärme in erheblichem Maße abgibt. In der zweiten Halbperiode, der soge- 
nannten Windperiode, wird alsdann durch das aufgeheizte Gitterwerk in entgegengesetzter 
Richtung der Kaltwind hindurchgeleitet, der aus den aufgeheizten Steinen Wärme ent- 
nimmt und damit selbst erwärmt wird. Alsdann wiederholt sich der Vorgang von neuem, 
und so fort in periodischem Wechsel. Offenbar wird sich, wie auch die Erfahrung zeigt, 
nach einer Ansahl von Wiederholungen des Prozeßes ein periodischer Temperaturzustand 
einstellen, bei dem die Temperaturen des Gases und des Windes periodische Funk- 
tionen der Zeit werden. Dieser im wirklichen Betriebe letzten Endes sich ein- 
stellende Temperaturzustand soll mathematisch eriaßt werden, das ist das 
eigentliche Problem. 

Zur mathematischen Formulierung dieser Aufgabe gehen wir aus von einem 
zylindrisch gestalteten Körper von beliebigem Querschnitt, dessen Achse wir zur z-Achse 
unseres Koordinatensystems machen, und dessen Temperatur wir mit u=u(%,y,z,t) be- 
zeichnen. Seine spezifische Wärme sei c, sein spezifisches Gewicht y, seine Wärmeleit- 
zahl 4. Die Temperatur des an ihm in der positiven oder negativen z-Richtung vorbei- 
streichenden Gas — bzw. Luftstromes sei dı bzw. ös. Diese Größe wird in dem uns 
interessierenden Falle eine Funktion von z und f, dagegen von x und y als unabhängig 
vorausgesetzt. Die Wärmeübergangszahl des Körpers gegenüber dem Gas- bzw. Luft- 


I) vergl. z. B. Nusselt, Die Theorie des Winderhitzers, Zeitschrift des Vereins Deutscher In- 
genieure (1927), Nr. 3, 8. 85. 

*) Insbesondere von der Wärmestelle des Vereins Deutscher Eisenhüttenleute in Düsseldorf, in 
deren Mitteilung Nr. 82 vom 27. Januar 1926 wohl die zur Zeit besten Beobachtungen über die in 
Rede stehenden Vorgänge gemacht worden sind. — Uebrigens habe ich in einigen Vorträgen gelegent- 
lich eines Hochsehulkursus der T. H. Breslau unabhängig von den genannten Arbeiten im Januar 1927 
an Hand von Beobachtungen von Harnickell (Hochofen-Ausschuß Nr. 66) eine Berechnungsmethode 
für Wärmespeicher mit sogenannter Semmelsteinfüllung entwickelt, deren Weiterführung und Ausgestaltung 
auf Speicher mit Gitterwerk die vorliegende Arbeit darstellt. 

3) Dies ist lediglich so zu verstehen, daß der Querschnitt durch den Wärmespeicher in horizon- 
taler Ebene unabhängig davon ist, in welcher Höhe der (Querschschnitt gelegt wird. Im Normalfalle 


besteht das Gitterwerk aus quadratischen Schächten, die durch Steinplatten von bekannter Dicke von- 
einander getrennt werden. 
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die Windperiode. Die z-Achse werde, der Richtung des Gasstroms entsprechend, nach 
unten positiv gerechnet; bei Eintritt des (Grases sei z2=0. 

Eine erste, auf Grund der praktischen Erfahrungen zulässige Vereinfachung der 
Ansätze!) besteht nun darin, daß wir die Wärmeleitung in der --Richtung vernachlässigen 
wollen; unter dieser Voraussetzung bestehen dann folgende Relationen zwischen den 
eingeführten Funktionen: 


strom sei &, bzw. “; ein für allemal bezeichne der lIadex 1 die Heizperiode, der Index 2 


u 1 On / 
= l 


Das zweite Glied der rechten Seite von Gl. (2) soll dabei die Berücksichtigung der Wärme- 
verluste enthalten, die bei jedem Wärmespeicher durch Abgabe nach außen entstehen. 
Um die Vorstellung durchführen zu können, daß die Gas- und Windtemperaturen über 
den ganzen (Juerschnitt binweg dieselben, also von & und y unabhängig sein sollen, wird 
der Sitz der Eutstehung dieser Verluste von der äußeren Begrenzungsfläche hinweg und 
eleichmäßig über die ganze Berührungsfläche des Körpers mit dem Gase bzw. Winde 
hin verteilt gedacht und außerdem der jeweils herrschenden Gas- bzw. Windtemperatur 


proportional gesetzt. Es bedeutet ferner «u die Temperatur an der Steinoberfläche; n be- 
zeiehnet die nach dem Gase zu gerichtete Normale der Steinoberfläche. In dem Falle, 
daß die Gastemperatur konstant ist, läßt sich hieraus bekanntlich die Funktion vu berechnen, 
wenn noch die Anfangstemperatur des Steines zur Zeit /=0 bekannt ist; die Lösung 
schreibt sich in der Form: 
a j-1 
Hierin stellt Z, eine von x und y abhängige Eigenfunktion der Differentialgleichung (1) 
mit der Randbedingung (2) dar; %, bedeutet den zugehörigen Eigenwert. Die Koeifi- 


zienten RB, sind durch die Entwicklung der bekannten Funktion u — (1 —_ ) ” zur Zeit 


a 
t—=( in der Form (1 -) = B, Z; gegeben; selbstverständlich hängen ihre 
Werte ebenso wie die von Z; von der Form des Querschnitts ab. 

Ist nun ® von ? abhängig, so gilt eine etwas allgemeinere Formel, die wir nun 
herleiten wollen. 

Wir gehen zu diesem Zwecke aus von der Darstellung der Funktion 1 durch die 
Eigenfanktionen Z, die etwa die Gestalt 


j-1 
haben möge; indem wir durch Variation der Konstanten die gesuchte Temperatur « in 
der Form 
ansetzen, ergeben sich durch Einsetzen dieses Ausdruckes in die Difierentialgleichung (1) 
und durch Multiplikation mit (4) die Relationen: 
o 
2 A;(t).e (1 - - (1 -#) 


J 


[0 


N 71 ot [24 J 1 ot 
Hiernach genügt es, die Funktionen 4,(f) den folgenden Gleichungen gemäß zu wählen: 
4,(t)= (1 — 
Diese Gleichungen lielfern nun für die gesuchten Funktionen 4;(t) die Ausdrücke 
t 
\ 89 
A;,(t) (1— dr+B, 


I) Jeh verdanke diese Angaben wie auch die Hervorhebung der Wichtigkeit des periodischen 
Betriebszustandes Herrn Dr. Schaek von der Wärmestelle Düsseldorf 


. 


ılen: 


schen 
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so daß sich für die Temperatur « die Gleichung ergibt: 
t 


ot j=1 a IT 


Hierin sind jetzt wiederum die Koeffizienten B, durch den Anfangszustand bestimmt. 
Uebrigens sei bemerkt, daß diese Größen von z abhängen werden, wenn der Anfangs- 
zustand und ® von zZ abhängen, was für die Gültigkeit der Differentialgleichung (1) 
natürlich keinerlei Schwierigkeiten macht. 

Während bisher die Gastemperatur ® als gegeben zu betrachten war, wird dies in 
unserm eigentlichen Problem nicht der Fall sein; die Gas- bzw. Windtemperatur bestimmt 


sich vielmehr erst durch die vom Stein aufgenommene bzw. abgegebene Wärmemenge. 
Die Gleichung, die diesen Uebergang regelt, lautet: 


10% 
Hierin steht links nach Multiplikation mit /z die an der Stelle z von der Mantelfläche 
des Zylinders in Höhe von Jz Metern abgegebene, rechts die vom Gase in der ent- 


sprechenden Schicht aufgenommene Wärmemenge 0 bedeutet die Gesamtoberfläche, 


u die mittlere Öberflächentemperatur, / die Höhe des Zylinders, % die in der Zeiteinheit 


zugeführte Gasmenge mal spezilischer Wärme des Gases, »® die Geschwindigkeit des Gas- 
stromes. 


Die Einsetzung von (5) in (6) liefert: 


Oz vvVt 


Diese Gleichung dient zur Bestimmung der gesuchten Gastemperatur, wenn der Anfangs- 
zustand zur Zeit = (0 bekannt ist; eine entsprechende Formel würde für die Wind- 
temperatur gelten. Z; bedeutet den Mittelwert von Z, an der Oberfläche. 

Wir haben nun zum Ausdruck zu bringen, daß zwar der Anfangszustand unbe- 
kannt, dafür aber der periodische Betriebszustand eingetreten sein soll. Zu diesem 
Zwecke suchen wir in der Formel (7) den Anfangszustand für z2—= 0 zu ersetzen durch 
den zur Zeit {= — p, wo p die Länge der Halbperiode bedeutet, alsdann diesen durch 
den zur Zeit {= — 2p, und so fort in infinitum. Bei dieser Verschiebung des zeitlichen 
Anfangspunktes ist zu beachten, daß in den auftretenden Integralen der Integrand jeweils 
mit dem Index | oder 2 zu versehen ist, je nachdem es sich bei dem betreffenden Inter- 
vall um eine Gas- oder Windperiode handelt, und daß außerdem infolge der Verschieden- 
heit, die die Werte (0) und 9, (0) und ebenso (p) und (p) gegeneinander auf- 


weisen werden, jedesmal ein entsprechendes Zusatzglied hinzutritt. So ergibt sich bei- 
spielsweise für den ersten Schritt. die folgende Rechnung: 


Nach Gl. (5) gilt zunächst, wenn wir von o bis p eine Gasperiode annehmen, 


nn 


In etwas allgemeinerer Schreibweise lautet die Gl. (5): 
t 


J £ 
fo 


Dabei ist 


Speziell gilt demnach für = — p und t—0: 
f —1 [44 IT 
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Die Subtraktion der Gleichungen (5) und (11) liefert nun durch Vergleichung der 
Koeffizienten nach Multiplikation der linken Seite mit (4): 


0 
also | 


Setzen wir diese Werte in die Gl. (5) ein, so ist damit der erste Schritt vollzogen. 
Allgemein ergibt sich auf die geschilderte Weise das folgende Resultat: 
0 —p 


— 2p 

| 

(( =t=p), 
wobei die Integranden jeweils die Funktionen e"}?W = ‚_ mit dem Index 1 oder 2 be- 
UT 


deuten. In abgekürzter Schreibweise und nach Summation der geometrischen Reihen 
erhält der Ausdruck die Gestalt: 


o 
- 2ak}?} 
p). 


Auf ganz analoge Weise erhält man für die Windperiode: 
t 


[04 [73 


J 
— 


ak?2y 


(pz=t=2p). 
Die Einsetzung dieser Ausdrücke in die Gl. (6) liefert: 
und für die Windperiode: 
t 
( 
— da + -) C 
— 
— ak?p 
— H:=0*% 2 r u0 
p .e ak? t—p) (15). 


Die hiermit entwickelten Gleichungen unseres Problems lassen sich noch umformen, 
wenn man rechter Hand die Integrale durch partielle Integration umgestaltet; dabei fallen 
alsdann die Zusatzglieder mit den Differenzen (9, — #,):—-o und (9, — Vy):--„ gerade weg, 
so daß eine formale Vereinfachung entsteht. Durch Einführung der neuen Variablen 

- 
„—=t-— an Stelle von fund 9 (z,f)—%(z,n) wird ferner die linke Seite zu ß: 
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die nun entstehende gewöhnliche Differentialgleichung für ® in : läßt sich durch Quadratur 
integrieren, wenn der Anfangszustand von # für — 0 bekannt ist. Auf diese Weise 


entsteht das in der Einleitung erwähnte System von Volterraschen Integralgleichungen: 


151 

0 

al 
2 


Wir werden indessen für den folgenden Zweck von dieser an sich wohl eleganteren 


Formulierung des Problems keinen Gebrauch machen, sondern uns an die Gl. (14) 
und (15) halten. 


2. Näherungsweise Lösung der Gleichungen für die mittleren Gas- und 
Windtemperaturen (Periodenmittel). Um eine praktisch brauchbare Lösung der 
Gl. (14) und (15) zu gewinnen, integrieren wir diese nochmals über ? zwischen den Grenzen 
o bis p bezw. p bis 2p und dividieren durch p, bilden also die Mittelwerte für eine Gas — 
bezw. Windperiode. Bezeichnet wı (2) bezw. ws (z) die mittlere Gas — bezw. Wind- 
temperatur an der Stelle x des Winderhitzers, so entsteht links 


Pı (2) [9 (2, 9ı | bzw. E 2)+ 19 (,0)— 
) 


Auf der rechten Seite können die Differenzen (9 — :-—o und — ohne 


einen erheblichen Fehler durch zo (z) — ws (z) ersetzt werden. Die Integrale der rechten 
Seite zerlegen wir in die drei Bestandteile 


— ip 
I: — III: Je«k?r dr+I+...: 
ÖT Ort 
0 —D — 3p —2p —4p 
p 0 —2p 
OT OT OT 
— 2p —p — 3» 


Umformung der ersten Bestandteile mit Hilfe partieller Integration liefert 


I: (2,t) (2,0) — dr 


bzw. I: (2,t)e — (z,p) — fa dr. 


p 
Die Bestandteile II und III der Integrale berechnen wir näherungsweise auf 


folgende Art: Zunächst fassen wir die ersichtlich entstehenden geometrischen Reihen 
zusammen und erhalten 


1: dr: je dr: 
p 
pP 
o —2ak?2p 
e j 1 
III: | dr euk?r dt» ——, 
ö IT 1 bzw. OT Da 


—p 
Diese Integrale werden nun in der Weise approximiert, daß wir in jedem den monoton 
wachsenden Faktor e“*; ® herausziehen und durch seinen mittleren Wert 


bzw. — 
ersetzen; alsdann werden die Bestandteile II und III: 
— ak: 
1 
[93 (2, 0) 9 (2,p)] bew. 
ak;ı'p 
ak?p 
1 
— : (2,p) — 9ı (2, 0)] 
ak?p [ ı ( p) 1 ( 
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kr 1 2a kr p 


1 
(2,0) — 9 (2, p)) 


Da nun jedes der Integrale noch mit e=“*,‘ multipliziert war, so ergeben sich 
aus II und III für die in Frage stehende Mittelwertformel die Bestandteile: 


a u 


a0 [7 Ci; Z; (i ak, P) e k; 
(akt 
a/j=1 
Die hierin auftretenden Summen 


[9 (2,p) — (z,0)) bzw. 
[93 (2, 0) (2, 


j=1 (ak? 
sind beide kleiner als 1. 

Der Bestandteil I liefert, in ähnlicher Weise behandelt, zunächst (Vorausziehen des 
ebenfalls normalerweise monotonen Faktors 9 nach dem zweiten Mittelwertsatz und Er- 
setzung der dann im übrigbleibenden Integral entstehenden Zwischengrenze durch { selbst!). 

(2,8) (2,0) — (2, 1) (et? — (2,1) — (2,0) bzw. 


nach Multiplikation mit e “', ‘ und Bildung des Mittelwertes entsteht näherungsweise 


l ak?p 
| 
bzw. — (z) — (z, p)] - ) 
«a 


Wir setzen: 


- 
ji p 
auch diese Konstante ist kleiner als 1. 

Alle bisher näherungsweise berechneten Glieder sind zusammengenommen ver- 
hältnismäßig klein, da die dabei auftretenden Differenzen sich zum größten Teil gegen- 
seitig zerstören. Von erheblicher Bedeutung ist demgegenüber das durch Mittelbildung 
aus dem zweiten Gliede rechter Hand in (14) bzw. (15) entstehende Glied. Es liefert, 
wenn die Differenzen, wie erwähnt, durch 0, (*)— ws (2) ersetzt werden, in beiden Gleichungen 


(2) — (2)]. 


Wir setzen die hierin auftretende Konstante 


107 (1 ‘ ) 
+e 


Somit ergibt sich insgesamt ein Gleichungssystem der Gestalt: 


(19). 


Pa (2) —- (=) — wy (2)] 


ao 
on — (dl "p) 


ew, 
ch, 


ch 


les 
N). 
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wenn wir noch die Annahme machen, daß die sämtlichen vorkommenden Differenzen 
0, (2,p) — dı (2,0) usw. den zugehörigen Mitteltemperaturen proportional sind. So setzen 
wir: 
wı (2) — (2, 0) =Yı (2) (2) (2, p) = Yyı (2) (20). 
(2, p) (2, 0) (2) (2, pP) —d; (2, 0) Ö, "ws (2) 


Mit diesen, ihrer Definition nach sämtlich positiven Konstanten yı, dı, ergeben sich 
nun die folgenden Ausdrücke für die Koeffizienten cı, dı, €, da: 


u WE 7 Br 


(1 +7 8.) 
u 
a=(1- £)(- +1) 


Wir haben nun endlich noch einen letzten wesentlichen Punkt zu berück- 
sichtigen: 

In unserer bisherigen Ableitung war vorausgesetzt, daß die Wärmeübergangszahl « 
in der Gasperiode denselben Wert hat wie in der Windperiode. Dies ist in Wirklichkeit 
nicht der Fall, sondern beide Werte sind namentlich in den oberen Regionen des Wärme- 
speichers aus Gründen, die hier nicht näher zu erörtern sind, deutlich voneinander ver- 
schieden. Eine eingehende Ueberlegung, die im einzelnen durchzuführen uns zu weit 
führen würde, zeigt nun, daß dieser Tatsache in der einfachsten Weise in erster Annähe- 
rung dadurch Rechnung getragen werden kann, daß in der ersten unserer Gleichungen 
der Wert «,, in der zweiten der Wert «; eingesetzt wird. Aus den gleichen Gründen 


muß ferner der Wert von « in ıı und 4, getrennt werden, von denen der erste der Gas- 
periode, der zweite der Windperiode entspricht. 


Demnach entsteht jetzt definitiv das folgende System von Differentialgleichungen 
für die gesuchten Mitteltemperaturen: 


wı (2) — dı 
(21), 
30 
w (2) = — 7°? . [ey wı (2) — da (2)] \ 
wobei die Koeffizienten die Werte haben: 
dh= )ıı+3-8) 
(22) 
2 Lg" 
19 Ba ng „2 


Die Werte von I”, 2)”, &,, sind dabei aus den Formeln (16), (17), 


und (19) mit Hilfe von «, bzw. «, zu berechnen. 


Zur Erläuterung und Prüfung der Theorie soll diese auf die Beobachtungen in 
der Mitteilung Nr. 82 der Wärmestelle des Vereins Deutscher Eisenhüttenleute angewandt 


werden. Die folgenden Rechnungen sind unter Mitwirkung von Hrn. cand. ing. Gigling 
durchgeführt worden. 


Der dort berechnete Speicher bestand aus quadratischen Schächten von der Oefinung 
0,01 m*, deren Platten eine Dicke von 0,05 m besaßen. Da sich die Bestimmung der Z; und 
der C; in diesem Falle streng durchführen läßt, so konnten auch die Speicherkonstanten 


x, 2’, 2”, 2” genau berechnet werden; das Resultat war auf drei Stellen genau: (p= ?) 


26* 


| 

| 
n- 
ng 
rt, 
en 1 
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4 0,438 0,566 | 0,357 0,105 
6 0,385 0,449 0,213 0.035 
10 0,287 0,302 0,095 0,005 
14 0,221 | 0,225 0,055 0,001 
18 0,173 | 0,173 0,034 = 
22 0,139 | 0,140 0,023 _ 


Von Interesse ist dabei insbesondere der Gang der Hauptkonstanten I mit «. Es 
zeigt sich, daß I für «= 0 den Wert 1 annimmt, dann mit wachsendem « abnimmt und 
für @& >= gegen 0 strebt. Demnach wird das Produkt «Z für «— 0 ebenfalls 0, nähert 
sich dann mit wachsendem « einem Maximalwerte und strebt bei @—>» gegen (. 
Dieses Verhalten ist besonders deshalb von Interesse, weil die Konstante 2, wie wir 


gleich sehen werden, auch für die Steinoberflächentemperaturen von der größten Be- 
deutung ist. 


Die weiterhin in den Koeffizienten (22) auftretenden Konstanten müssen — dies 
ist eine Unvollkommenheit der Methode, die durch weitere Untersuchungen zu beheben sein 
wird — auf Grund der vorliegenden Erfahrungen über den Anstieg der Gastempe- 


raturen, den Abfall der Windtemperaturen und die Abweichung der mittleren von den 
Anfangs- und Endtemperaturen in beiden Halbperioden (vergl. die Formeln (20) geschätzt 
werden. Im vorliegenden Falle wurden die Werte yı = 0,2, ya = 0,145, dı = 0,203, 
d3 — 0,05 gewählt. 


Was nun endlich die Wahl von u anbetrifit, so ist zu beachten, daß infolge der 
Verschiedenheit von @ in den verschiedenen Schichten des Wärmespeichers eine Konstante 
in 2 kaum der Wirklichkeit entsprechen dürfte; es wurde im vorliegenden Falle der 
Quotient ” von z unabhängig gewählt, und zwar "' — 0,02, I = 0,03. Die übrigen in 

a ag 
den Formeln vorkommenden Konstanten sind der zitierten Mitteilung entnommen. 


Mit diesen Werten ergaben sich nun unter der Annahme einer Abhängigkeit der 
Wärmeübergangszahlen « von z, die ebenfalls im Rahmen der Beobachtungsfebler der 
Wirklichkeit entsprach, die folgenden Werte für die Koeffizienten: 


a, 2,0 0 ag 220 a9 La 
| dy ag ds 
6 0,216 0,202 4 0,142 0,124 
10 0,268 0,248 0,199 0,157 
14 0,294 0,267 10 0,260 0.168 
15 0,302 0,268 
22 0,304 0,263 


Wie man sieht, sind die Koeffizienten in Wirklichkeit nicht so stark veränderlich, 
wie es die durch die Beobachtungen nahegelegte starke Veränderlichkeit von « eigentlich 
erwarten lassen sollte; wenn es also nur auf die Bestimmung von w: (z) und ws (z) an- 
kommt, könnte man mit genügender Genauigkeit auch mit mittleren Werten rechnen und 
somit @, und s als Summe von zwei Exponentialfunktionen darstellen. 


Von Bedeutung werden diese Verschiedenheiten indessen, wenn man sich ferner noch 
für die mittleren Steinoberflächentemperaturen U, (z) und U; (z) interessiert. Diese ergeben 
sich nämlich aus unseren Gleichungen (6), (7) und (21) vermöge der Bezeichnungen 
c, und d, aus Gleichung (22) in der Form: 


— U, = (cı wı — dı we), 

und hier spielt die Veränderlichkeit von Y schon eine beträchtlichere Rolle. 
Rechnret man nun mit diesen Werten mit Hilfe der Runge-Kuttaschen Regel 
diejenigen Integrale unseres Systems von Differentialgleichungen aus, deren Anfangs- 


werte in zw, (2) für z=0 und in wm (2) für z=1! den Beobachtungen entsprechen, so er- 
gibt sich folgende Tabelle: 
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(2) w; (2) U] (2) UL (2) wa (z ws (2) Ug (2) Ug (2) 
beob. berechn. beob. berechn. beob. | berechn. beob. berechn. 
0 1125 1140 1100 1120 00 s20 950 980 . 
1 1000 1027 930 984 700 717 55 859 
3 800 730 775 530 565 690 701 
5 625 637 580 607 425 446 530 553 
7 530 507 465 474 340 351 420 437 | 
N) 440 410 375 376 280 279 335 350 ! 
11 355 334 300 296 215 216 270 276 
| 13. 280 269 220 228 160 156 200 210 | 
15 210 211 160 169 105 100 130 146 ! 
17 150 160 100 116 65 50 100 91 | 
» Angesichts der Unsicherheit der Beobachtungen kann hier von einer vollständigen | 
Uebereinstimmung nicht nur der mittleren Gas- und Windtemperaturen, sondern auch der | 
Steintemperaturen, und zwar sowohl in der Gas-, als auch in der Windperiode gesprochen | 
h werden. Diese ausgedehnte Probe aufs Exempel kann als eine starke Stütze für die 
Zulässigkeit unserer Vereinfachungen und damit des Gleichungssystems (21) an- | 
i gesehen werden. Es werden dabei auch die interessanten Beobachtungen der Wärme- | 
stelle über die Annäherung von Stein und Gastemperatur sowie über die starke Ver- 
schiedenheit von Stein- und Windtemperatur in den oberen Schichten des Speichers voll- | 
kommen bestätigt und durch den Gang der Hauptkonstanten I erklärt. 
Auch die gespeicherten bezw. in Verlust gegangenen Wärmemengen können aus 
unsern Formeln berechnet werden und ergeben im großen und ganzen Üebereinstimmung | 
mit den Beobachtungen. Im einzelnen sind diese Angaben noch einer Diskussion unter- 
worfen, weil die richtige Einschätzung der Verluste infolge der Mitwirkung des Brenn- 


schachtes und infolge unserer Vernachlässigung der Wärmeleitung durch den Speicher 
\ von oben nach unten schwierig ist. 

Es sei endlich bemerkt, daß der Vergleich unserer Formeln mit den Beobachtungen 
indirekt die Abbängigkeit der Wärmeübergangszahlen von zZ zu ermitteln gestattet, wenn 
auch in dieser Beziehung eine explizite Auflösung unserer Formeln nach & nicht vorliegt. 
Soviel kann jedenfalls auf Grund der Theorie und der Beobachtungen gesagt werden, 
daß die Annahme eines konstanten Wertes von « für die Gasperiode ausscheidet, daß 
vielmehr ein starkes Ansteiıgen nach der Höhe zu angenommen werden muß, um die 
erwähnten Erscheinungen bezüglich der Steinoberflächentemperaturen zu erklären. 866 | 

Breslau, den 23. Dezember 1927. ' 


Zur graphischen Bestimmung | 
der Knotenpunktsverschiebungen räumlicher Fachwerke. 
Von A.W,SOTOFF in Moskau. 


ie folgenden Betrachtungen stützen sich auf die Abhandlung von R. v. Mises 
»Graphische Statik räumlicher Kräftesysteme« (Zeitschr. f. Mathematik u. 


Physik Bd. 64 [1916]), deren Inhalt als bekannt vorausgesetzt wird. Die dort ein- 
| geführten Bezeichnungen und Terminologie werden ohne Erläuterung gebraucht'). | 


1. Die Aufgabe. Es sei P (Abb. 1) ein Knotenpunkt eines räumlichen Fach- 
| werkes, der durch die Stäbe (Aı) (As) (As) mit drei anderen Knotenpunkten (Pı) (P;) 
| verbunden ist. Sind die Verschiebungen der Punkte (P,ı)(P,)(P;) und die Längenände- 


!) Nachdem der vorliegende Aufsatz fertig war, ist mir die Arbeit des Hrn. W. Prager: >»Bei- 
trage zur Kinematik des Raumfachwerkes« [diese Zeitschr. Bd. 6 (1926), Heft 5] bekannt geworden. Ich 
behandle nur das, was in Heft 5, S. 345, von Hrn. Prager behandelt wird. Dabei verwende ich nur 
die ursprünglichen Mayor-Misesschen Abbildungsverfahren, während Hr. Prager neben diesen (B) 
Verfahren noch ein zweites (4) verwendet. Dazu möchte ich nur bemerken, daß das Wesentliche 
| im (A) Verfahren, der Bildpunkt (zo yo) Nr. 1, S. 343, nichts anderes als der Antipol des Bildstabes des 

(B) Verfahrens ist Deswegen scheint mir die von Hrn. Prager (S. 342) ausgesprochene Vermutung, 
daß Hr. B. Mayor die Tragweite dieses Elementes nicht erkannt habe, kaum gerechtfertigt, da er ge- 


rade die Antipole zur konstruktiven Bestimmung der Knotenpunktverschiebungen in der unten zitierten 
Arbeit verwendet. 
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rungen der Stäbe ermittelt, so lassen sich die Projektionen der Verschiebung des Punktes (P) 


auf die Stabrichtungen bestimmen. 


Unsere Aufgabe wäre demnach, den Verschiebungs- 


vektor (B) des Punktes (P), wenn die Projektionen desselben auf drei Richtungsgeraden 
(A) (As) (A;) [die nicht komplanar sein dürfen] bekannt sind, graphisch zu bestimmen. 


|RABSOZI 


Abb. 1. 


4 


2. Hilfssätze. Es sei (Abb. 2) « der Antipol der Geraden A in bezug auf einen 


Kreis von Halbmesser e und Mittelpunkt O, so daß 


ist. Dabei liegt « auf der Verlängerung des von O auf A gefällten Lotes Oa. Wir 
wollen die Gerade A »Antipolare« des Punktes « nennen. Es ist leicht, die folgenden 


Sätze zu beweisen: 


I. Die Antipolare \ des Punktes a geht durch «. 


2. Der Antipol # irgend 
einer durch a gezogenen 
Geraden B liegt auf N. 

3. Der Antipol y irgend 
einer durch 5 gezogenen 
Geraden € liegt auf B. 

Diese Sätze folgen di- 
rekt aus den Bezie- 


hungen 
Va: Ob-O0P 
=0c:07 (1). 


4. Wenn wir die Lote «a 
und @«b auf A bzw. B 
fällen, so wird 

ob 


Ba ad Oa 03 (2). 

(Ay) 
78) (C2) 
(RaB5023) (C,) 

Abb. 3, 


3a. 


a’ 
\ 
| (A;) \ 
\ 17 
\(A2) 
\ op 
\ 3 N 
\ N 
Abb. 2. 
B\ 
A, | 
1: e 
a N o 
\ 
A 
| | 
| \ A, 
(RA85023«] | \8 
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3. Auflösung. Erstes Verfahren. Es seien (Abb. 3): (B) der Verschiebungs- 
vektor; (A,), (4,) seine Projektionen auf zwei gegebene Richtungsgeraden (Aı) bzw. (4,); 
(Cı), (02) Vektoren, die die Endpunkte von (A,) und (B) bzw. (4,) und (B) so ver- 
binden, daß 2 

(B) = (41) + (Cı1) = (42) + 
ist. (Cı) und (C,) sind dann die Richtungsgeraden der Vektoren (C}) und (03). Abb. a 
zeigt nun die Abbildung auf der Ebene z, und zwar sind 
Aı, As, B die Bilder der Geraden (A,) (4;) (B), 
B das Bild des Vektors (B) (bekannt vorausgesetzt). 

Es mögen «,, &, die Antipole der Geraden B sein; der Durchschnitts- 
punkt Aı, As; die Antipolare von as, auf der «&, liegen; « der Durchschnittspunkt 
B, Yo; A die Antipolare von «, auf der ß liegt. 

Es ist jetzt leicht, die Bilder C;, ©, der Geraden (C}), (Ca) zu finden. Die Ge- 
raden (Cı), (C;) sind nämlich komplanar mit (Aı), (B) bzw. (As), (B) und senkrecht zu 


(41) bzw. (As). Ihre Bilder C,, © müssen dureh die Durchschnittspunkte Aı, B bzw. 
Ar, B und durch die Antipole «, bzw. «& gehen, wodurch ihre Lage bestimmt wird. Da 


(B) die Summe von (Aı)(C}) resp. (43) (5) ist, so muß, nach dem II. Hauptsatz von 
R. v. Mises, B die Resultante von Aı, Cı bzw. A:, C, sein, wenn Aı, Cı bzw. A., C, die 


Bilder der Vektoren (4,ı) bzw. (4A,) (C,) darstellen. 

Eine Anwendung des Parallelogrammsatzes führt zur Konstruktion der Vektoren 
Aı, Ci, Aa, Ca. Wir wollen nun den Momentensatz anwenden. Mit «, und dann «& für 
Bezugspunkt bekommen wir nacheinander 


B b' dı dı und (9 b,' = IAs| Ay (3), 


wo B, Aı, die Beträge der Vektoren und & aı, die Längen der 
Lote von «&, bzw. @ auf B, A, und A, bedeuten. Aus diesen Gleichungen finden wir 


A| ayaı 


(4). 
As ba Ay 
a 
Da aber 
a2 acag 


Daraus ersehen wir, daß, wenn A, und 4A, gegeben sind, so läßt sich der Durchschnitts- 
punkt « der gesuchten Geraden B mit VW konstruieren. 

Nun liegt aber eine Zweideutigkeit vor. Wir können nämlich auf \ zwei Punkte 
# finden, die der Gl. (4) genügen: « kann auch zwischen «, und @ liegen. Diese Zwei- 
deutigkeit wird durch die folgende Regel gehoben: Punkt « liegt außerhalb «,, «, 


wenn 4ı, A; beide zum Punkt a oder beide davon gerichtet sind; « liegt 


zwischen &,@«, wenn von 4ı, 4; eine zum Punkt m, die andere davon ge- 
richtet ist. Es ist leicht, die Richtigkeit dieser Regel zu prüfen. Wie ersichtlich, liegt 
der erste Fall für die angenommene Lage von Aı, 4, vor, und die Regel ist dabei be- 


stätigt. Wenn wir statt 4, die Gerade A; nehmen (mit «& als Antipol) und A;, C3 finden, 


so sehen wir, daß für A, und 4, der zweite Fall vorliegt, und die Regel ist dabei wieder 
bestätigt, da jetzt « zwischen «, und @ liegt; denn die Reihenfolge der Antipolen auf \, 
wenn eine Gerade um a, nach dem Uhrzeiger gedreht wird, ist &, &, &, ©, @, ds. 
4; ist hier bloß eine Hilfsgerade und muß nicht mit dem Bilde A, der dritten gegebenen 


Richtungsgeraden (4A;) (das nicht durch «, geht, weil (Aı) (Aa) (As) nicht komplanar sind) 
verwechselt werden.| Wenn wir jetzt die dritte Projektion (43) mit (A4,) in derselben 


Weise benutzen, so finden wir einen zweiten Punkt « der gesuchten Geraden B, und 
damit ist die Aufgabe gelöst. 
Zweites Verfahren. Aus (2), (3) und (4) finden wir 


aaı 02 Oay 


1. | 
| 
| 
1 
| 
| | 
ir 
n 
| 
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| ayay 
und Baı:Par' — |Aı : (6), 
Oaı Oag 


wo ?aı, Pay die Lote von $ auf Aı und 4; sind. Die letzte Gleichung gestattet, eine 
Gerade A zu konstruieren, auf der 2 liegen muß. Die entsprechende Zweideutigkeit wird 
gehoben, wenn wir bemerken, daß A außerhalb Aı, 4; und zwischen Aı, 4; liegt. Nach 
diesem Verfahren wird also der Antipol 5 der gesuchten Geraden B als Durchschniitts- 
punkt zweier Geraden gefunden. Es sei hier erwähnt, daß Hr. B. Mayor (C.R., Paris 
1915) auf analytischem Wege zur Bestimmung der Antipole $ gekommen ist. Dabei 
deutet er dieselben als Momentanpole ebener Scheiben (plaques), die den Raumknoten- 
punkten auf der z-Ebene entsprechen. Diese Deutung könnte aus (5) und (6) leicht ab- 
geleitet werden, aber wir sehen keinen Vorteil darin. 


Drittes Verfahren. Es seien (Abb. 1) 


(Aı') (Ay') Vektoren, von denen der 
erste auf (A,) und (A;), der zweite auf (A;) 
und (A,), der dritte auf (A,) und (43) senk- 
3 recht steht, und deren Beträge so bemessen 

sind, daß ihre Projektionen auf (Aı) (A,) (45) 


resp. gleich den gegebenen (A,), (A»), (4A;) 


sind. Dann wird (B) = (Aı') + + 
4 a, und sein Bild B zur Resultante der Bilder 
Aı', Ay’. Es handelt sich also um die 
KanN\ a, Konstruktion dieser drei Bilder. -- Es seien 
A; (Abb. 4) Aı, As, A;, &ı, die Bilder auf 
der z-Ebene, der drei Richtungsgeraden 
TEN (Aı) (43) (45) und ihre Antipole. — Die 
L \ 4 Bilder 4;' der Geraden (A;') (4;') 

| 


müssen resp. durch und @, und &; 
| “ und « gehen und werden dementspre- 
A, chend gezggen. 

Wenn nun (A), (Al), Vektoren 
sind, die die Endpunkte von (A,) und (Aı,') 
usw. so verbinden, daß 


(4)=(4)+ (U); (4)= (4) + U); 
(43) = + (ls), 
so schließen wir, genau so wie vorher, daß 
das Bild der Geraden (\I,) durch (den 
Abb. 4. Durchschnittspunkt Aı, Aı und durch «ı 
sehen muß, und dementsprechend |, und 5. 
Nachdem die Bilder Yı, U, WU; gezogen sind, führt die Anwendung des Parallelogramm- 
satzes zur Konstruktion der gesuchten Aı', Ay’, As’. Es bleibt ibre Resultante zu finden. 
850 
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Über die Durchrechnung der Größen zweiter Ordnung 


durch ein optisches System. 
Von M. HERZBERGER in Wetzlar. 


s ist das Verdienst Allvar Gullstrands, in allen seinen Arbeiten darauf hingewiesen 
zu haben, daß ein Fortschritt in der Theorie der optischen Instrumente nur dadurch 
zu erwarten ist, daß man die differentialgeometrischen Eigenschaften der in der Optik 

allein zu behandelnden Normalsysteme einer eingehenderen Betrachtung unterzieht. Es 
soll im folgenden gezeigt werden, wie man mit Hilfe vektoralgebraischer Methoden einige 
der für die Praxis wichtigsten Probleme Gullstrands ohne allzuviel Rechnung an- 
schaulich lösen kann, und es sollen die Ergebnisse so umgewandelt werden, daß sie für 
das praktische Rechnen Bedeutung gewinnen können, 


Az 
2, 
/ 
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Gegeben sei die brechende Fläche «a(o,?7) und eine zweifach unendliche Mannig- 

faltigkeit einfallender Strahlen 

a(o,r)-+A8(o,r). 

8(o,r) sei der Einheitsvektor entlang den Strahlen. Wir greifen einen Strahl heraus, 
den wir als Hauptstrahl unseres Bündels ansehen wollen, und fragen nun, wie ändern 
sich die Eigenschaften unseres Strahlenbündels entlang diesem Strahl, wenn die Fläche 
a(o,r) das einfallende Strahlenbüschel bricht. Als Gesetz nullter Ordnung bezeichnen 
wir das Brechungsgesetz. Im allgemeinen wird jeder Strahl eines Normalbündels zwei 
Flächen berühren, die sogenannten kaustischen Flächen des Strahlensystems (auch Brenn- 
flächen genannt, doch ist dieser Name in der Optik schon vergeben). Die Gesetze erster 
Ordnung geben uns die Berührungspunkte des Strahls mit diesen kaustischen Flächen, 
ausgedrückt durch die entsprechenden Größen für das einfallende Strahlenbündel und 
die brechende Fläche; die Gesetze zweiter Ordnung geben uns unter anderem die Gauß- 
sche Krümmung der beiden Brennflächenmäntel in den Berührungspunkten. 

Es ist nun die Aufgabe dieser Arbeit, für einen speziellen Fall, nämlich, wenn ein 
Hauptschnitt der brechenden Fläche mit einem Hauptschnitt des einfallenden Strahlen- 
bündels entlang unserem Hauptstrahl zusammenfällt, die Gesetze erster und zweiter Or- 
dnung vollständig aufzustellen. Dieser Fall ist praktisch von besonderem Interesse, da 
er sowohl den Fall behandelt, wo der Hauptstrahl in einem System von rotationssymmetrischen 
Flächen und gleicher Achse die Achse schneidet, als auch den Fall, wo in einem Prisma 


der betrachtete Hauptstrahl im sogenannten Hauptschnitt liegt, d.h. in einer Ebene senk- 
recht zur Prismenkante. 


1. Das Brechungsgeseiz (Gesetz nulliter Ordnung). Für die zur Ableitung 
unserer Sätze benötigte Betrachtung führen wir als Hilfskoordinaten die Bogenlängen © 
und 7 der beiden Scharen von Krümmungslinien auf der brechenden Fläche ein, und 


zwar sei a- der Tangentenvektor in der Einfallsebene, während a; auf der Einfallsebene 
senkrecht steht. 


Aus obigem folgt, daß 
ist. o sei der Einleitsvektor entlang der Flächennormale. Den Winkel © zwischen o 
und $ bezeichnet man als Einfallswinkel. 
Dann läßt sich das Berechnungsgesetz in der Form 


schreiben, wo das Zeichen -/ angibt, daß die Differenz des in der Klammer stehenden 
Ausdrucks vor und nach der Brechung zu bilden ist; also 


Gl. (2) kann man, da sie unabhängig von ©, 7 gilt, in die beiden folgenden Gleichungen 
zerlegen 


daraus folgt, daß der Vektor (An&) die Richtung von vo hat, also es ist 


wo I’ die sogenannte astigmatische Konstante, sich durch skalare Multiplikation von (3) 
mit 0 zu 


ergibt. 


2. Die Gesetze erster Ordnung. Aus Gl. (2) erhält man nun einfach durch 
Differenzieren die Gesetze höherer Ordnung. Gl. (2) bestimmt, wenn der einfallende Strahl 
und die brechende Fläche gegeben sind, den gebrochenen Strahl eindeutig; wir bekom- 
men also so ein vollständiges System von Beziehungen. Wir finden aus (2) 


oder unter Benutzung von. (4) 


Da (6) unabhängig von ©, 7 gilt, ergibt sich das Gleichungstripel: 
Ans-a-—= — And = (7). 


Wir brauchen nur noch die in (7) vorkommenden Größen erster Ordnung geometrisch 
zu deuten. Wir tun das, indem wir alle vorkommenden Vektoren in dem mit der Fläche 


| 
h. 
rd | 
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ei 
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| 
1) 
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3) 
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3) 
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verbundenen rechtwinkligen Koordinatensystem mit den Einheitsvektoren DO. d,, d- aus- 
drücken. Hier machen wir von unserer Voraussetzung Gebrauch, daß die Einfallsebene 
für unser einfallendes Strahlensystem wie für die brechende Fläche Hauptschnitt ist. 
Dann folgt aus Sätzen der Strahlengeometrie, wenn s, ? die beiden Hauptkrümmungsradien 
der einfallenden Fläche, vom Durchstoßpunkt an gemessen sind, analog r,, r, die Haupt- 
krümmungsradien der brechenden Fläche sind 


(1 X 0) (8) 
nun liegt aber vo. bzw. ®8- in der Einfallsebene (durch vo und a-) und zwar 0- senkrecht 


zu 0, also parallel a-; 0; bzw. &, liegt senkrecht dazu, also in der Richtung von \:. 
Wir können also setzen 


&-= un + Pı-, \ 
Hier soll der Index »Zs kennzeichnen, daß wir das zur brechenden Fläche gehörige 


Normalsystem meinen. Skalare Multiplikation mit o bzw. a-, ü; und vektorielle Multi- 
plikation mit 0 bzw. s gibt unter Beachtung von (8) und (1) 


I 1 1 
"t Ts 
t t 5 
Aus 
(12) 
folgt durch Differenzieren unter Berücksichtigung von (9) 
Aus (9), (10), (11) folgt 
g- = P, Is 0, = Y . 


Also gibt die zweite Gl. (7) eine Identität 0 = 0, während man aus den übrigen erhält 


An®=T Br, 
Das sind die bekannten Gleichungen zweiter Ordnung. In ihnen ist 
1 1 1 
P= . Y zu Pr = Ir = (15). 
t N: 


Die Durchrechnung geschieht, wie folgt. Gegeben ist /, s, 7. rs Einfallswinkel :. Aus 
(15) berechnet man fi, 7, Pr Yr. Aus (14) folgt dann 5, y' nach der Brechung. Die 
Lage der Berührungspunkte mit der Kaustik nach Brechung ergeben sich aus (15) zu 


Nun kann man weiter durch die nächste Fläche rechnen. 


3. Die Größen zweiter Ordnung. Man kann hier auch wieder von Gl. (2) 
ausgehen und durch Differentiation eine vollständige Reihe von vier Differentialglei- 
chungen erhalten. Die Rechnung lehrt, daß die Differentiation von (14) zu denselben 
Resultaten führt. Wir erhalten also | 

=Ilyr-+Il:yr, 

In den Formeln (17) soll der Index 0 bezw. z in bekannter Weise die Differentiation 
nach diesen Größen andeuten. 

Nun ist nach (5) 


Ancos’—- An (v$) 
Es ergibt sich aber 


dcos dcosi 
—= + 0,8 +0% 
unter Benutzurg von (10), (11) 
deosi Ira sin’ sin? cos? deosi 
= = + === () 
( rt t do 


ht 
Is. 


ge 


ält 


Jie 
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also T=4 —nsinid” 
t 
n sin? 
weil A = Ö 


ist (Brechungsgesetz). 
Das gibt in (17) eingesetzt 
t 
Die Gleichungen (18) enthalten eigentlich schon die Größen zweiter Ordnung. Es 
handelt sich nur noch darum, die vier Größen #-,ßs,y-,)’; in Beziehung zu den geometrischen 


Größen unseres Strahlensystems zu bringen, um Invarianten für die Weiterrechnung zu 
erhalten. 


Wir finden zunächst aus (11) 


t t dr t do s? s’ 
Nun ist nach der obigen Rechnung 
d cos dcosi sind cosi sin 
t | (20) 
cos? it- sin 27 fcosi cos?itz 8- \ 
und umgekehrt 


Es bleibt uns noch, die geometrischen Größen unseres Strahlensystems durch die 
vier Größen (21) auszudrücken. Zu dem Zweck verfahren wir, wie folgt: Wir betrachten 
die beiden kaustischen Flächen. Wir können schreiben 


. . 0.0.22). 


Dann wird bekanntlich die Gaußsche Krümmung der beiden Brennflächen 
gegeben durch 

= — 23). 
bi +? bes? — (ber bio)? — 5)? 3) 

Hier sind n, bezw. ı1, die Flächennormalen der Kaustischen Flächen im Berührungs-' 
punkt mit dem Hauptstrahl. 

Wir betrachten daneben noch die Krümmurgsradien der Kurve o — konst. auf der 
ersten Kaustischen Fläche, bezw. 7 = konst. auf der zweiten kaustischeu Fläche und finden 
für die Krümmung dieser Raumkurven, die bekanntlich geodätische Linien auf den beiden 
Brennflächen sind 

b; b, 

Diese vier Größen wollen wir nun zu den Größen (21) in Verbindung setzen. 


Wir betrachten zuerst die erste Brennfläche. Wir drücken a- und a; durch 8- und 8; aus 
und erhalten 


d- = sin . . . . . . . (25) 
also d,-= (sinö-+1.)$, — . 0. (86) 
daraus +? = (sind + t-)?, (sini+t:) t;, 


Aus (27) ergibt sich, unter Benutzung von (10) und (11), daß der Einheitsvektor 
in Richtung der Flächennormalen die Richtung 


hat. 
Aus (26) folgt weiter, unter Berücksichtigung von 


PW, j 
1S- 
ne 
st. ! 
pt- 
| 
| 
| 
| 
us 
| 
(2) 
ei- 
en 
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daraus b, gr = &, . . . . . . . . . (30) 
ts ts cos?i 8- 1 
(t—8) t— 8 8? (81), 


aus 6; d&-—= 0 und (31) folgt 


t; cos? i 8: 1 
Wir finden nun 
| t sind + t- t 
Di-- : bi == cos 2, = bı $- 
cos i t cos? (33) 
co8ı 8” cosı 
daraus 
(sin? -+t-) (sin! +)? (t- 9°? 
und endlich aus (23) und (24) 
1 8- 
K,= — 35 
(sind + t-) (35) 
( t- t 
Wir betrachten nun die zweite Brennfläche 
Es wird 
+ (s— = 8; 8 . . . . . . (38). 
Daraus 
9 ( N? g? 
83° (8 — cos?i 
b, >)2 . . . . . . . (40). 


Die Flächennormale der zweiten Brennfläche ergibt sich aus (38) zu 
und aus (38), (32) folgt 


tz s cos?i 


b, N, = (s t) = — N, = 0 b, N, . (42), 


8 
85 ) 


Damit haben wir das zur Durchrechnung nötige Formelmaterial gewonnen. Wir 
wollen jetzt einmal den Gang der Rechnung verfolgen. Gegeben sei für das eintreffende 
Strahlenbündel X‘, R,, K,, R, ebenso die entsprechenden Größen für die Normale der 
brechenden Fläche. (resucht die entsprechenden Größen für das gebrochene Bündel. 


Aus (35), (36), (44), (45) ergibt sich zuerst 


Rs cos i Ks (t — s) Rı cosi 
t t 
analog 
Rsr RıF Rır RırF 
Ts 
daraus aus (20) 
eos? i sin 27 /cosi 1 rt7 co3? rt3 
— \ Pr: = Ps 2 Pr: = > 
t t t rt (48) 


); 
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erhalten wird dann aus (18) Y-, durch 


rt t 


aus (21) 
— + = b = 8? . (50) 
und zum Schluß aus (35), (36), (44), (45) 
’ 


Damit ist die Durchrechnung der Größen zweiter Ordnung von einer Fläche zur nächsten 
beendet. Wir setzen noch zwei leicht verifizierbare Kontrollformeln her. Es gilt identisch 


Br — ‚cos isini 
Ps = B-y Prs, I = 
Man erhält diese Formeln, indem man mit Hilfe von Gl. (9) und (10) a;- im Koordinaten- 
system D, ü-, A, bildet, und die gefundenen Werte vergleicht. 
Man erhält aus (9) 
Prs Ir + Pr Is = 
also 


= — - Ns . . . (52). 
Aus (10) folgt 


also 
as 35 a yF 

— — — + I — - 

Aus (52) und (53) folgen die gesuchten Beziehungen 

ßs Fr YFt 

0 = = — 54). 
Br-yr’ PB—r Br-yr P-y (54) 


Zum Schluß wollen wir noch einige spezielle, besonders interessierende Fälle be- 


handeln. Gegeben sei ein System von Rotationsflächen mit gemeinsamer optischer Achse. 


Wir betrachten einen Hauptstrahl, der die optische Achse schneidet, und ein ihm ent- 
haltendes Normalenbündel, das irgendwann von einem festen Punkt herkam. Dann 


waren anfangs alle Größen zweiter Ordnung gleich null, also $0s,fos,to-,80-—=0. Da 
es sich um ein Rotationssystem handelt, ist R;r = 0, also 

daraus Psr=Ysr 0 und 
endlich 


Wir sehen also, unter den obigen Voraussetzungen schneidet die sagiıtale Bild- 
schale die Meridianebene stets in einer isolierten Kurve. Das Formelsystem reduziert 


sich von vier Formelgruppen auf deren zwei, die wir hier zum Schluß noch einmal zu- 
sammenschreiben. 


4. Formeln für die Umgebung der Meridianstrahlen eines rotations- 
symmetrischen Systems. Gegeben A,, K., gesucht 


Rı cosi 


. Kı s)? 
— sini, Recosi. . .. (1) 
sin?! ,cosi sin 
Ts 
„2 
’ Br t ’ cos? 1 ’ 92 
tz +2 t tg ( — ). = . . . . (4), 
(sin + 1 8-' 


| 
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Für Kugelflächen wird insbesondere nur in Formel (3) 

während sich die übrige Rechnung im allgemeinen nicht vereinfacht. 


Noch ein Fall ist von besonderem Interesse: Der Hauptstrahl bilde die Achse 
einer Rotationsfläche, dann wird auch 


d.h. die Meridiankaustik hat auf der Achse einen isolierten Punkt bzw. eine Spitze. 879 


(6), 


Über die Berechnung der Temperaturverteilung 


im Rotor einer Dynamomaschine. 
Von STEFAN BERGMANN in Berlin. 


(Aus dem Institut für angewandte Mathematik an der Universität Berlin.) 


n der vorliegenden Arbeit wird unter Zugrundelegung der in 1 besprochenen verein- 

fachenden Annahmen ein Verfahren zur Berechnung der Temperaturverteilung im 

Rotor einer Dynamomaschine angegeben. Wie in 1 gezeigt wird, reduziert sich infolge 
der angegebenen Vereinfachungen die Frage darauf, eine im Innern eines Kreiszylinders 
reguläre harmonische Funktion :w, die an der Berandung des Zylinders vorgegebene Werte 
annimmt, zu bestimmen. 

Unter Benutzung eines an einer anderen Stelle bewiesenen Satzes wird in 2 ge- 
zeigt, wie man eine durch die Randwerte gegebene harmonische Funktion in einem Kreis- 
zylinder durch Polynome approximieren kann. Die sich auf diese Weise ergebende 
Methode der Approximation ist rechnerisch verhältnismäßig einfach und eignet sich deshalb 
zur zahlenmäßigen Behandlung unserer Aufgabe. In 3 wird eine Majorante für den 


Fehler, den man begeht — falls man die Funktion durch ein entsprechendes Appro- 
ximationspolynom ersetzt — angegeben, und schließlich in 4 ein Beispiel zahlenmäßig 
durchgeführt!). 


1. Mathematische Formulierung der Aufgabe. Der Rotor einer Dynamo- 
maschine ist ein kreiszylindrischer Eisenkörper, in denen Mantel- und teilweise Stirn- 
flächen-Kupferdrähte eingebettet sind, die vom elektrischen Strom durchfilossen werden. 
Der Rotor befindet sich auf einer Welle, die den Rotationsantrieb übermittelt. Bei der 
Arbeit entwickelt sich sowohl in den Kupierdrähten, wie auch in dem Eisenkörper selbst 
eine beträchtliche Wärmemenge, die eine Temperaturerhöhung bewirkt. In der Regel 
sind wir im Stande, die Temperaturerhöhung an der Berandung des Rotors experimentell 


Abb, 1. Längsschnitt durch den Rotor einer Dynamomaschine, 


zu bestimmen (oder zu berechnen), ferner läßt sich (verhältnismäßig genau) diejenige 
Wärmemenge berechnen, die im Innern infolge der Eisen- und Hpysteresis-Verluste, wie 
auch infolge des wechselnden magnetischen Feldes entwickelt wird. Es entsteht die 
Aufgabe, aus diesen Angaben die Temperaturverteilung im Innern des Rotors zu berechnen. 


I) [eh benutze die Gelegenheit, Herrn Dr. M. Pin! für die Hilfe bei der Durchführung der 
numerischen Rechnungen meinen besten Dank auszusprechen. 


eh. 


se 
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Wir sind aber hierbei gezwungen, eine Reihe von vereinfachenden Annahmen zu 
machen: 


I. Das Material, aus dem der Rotor hergestellt ist, ist homogen. Insbesondere die 
Welle, auf der der Rotor aufgekeilt ist, besitzt dieselbe Wärmeleitfähigkeits- und Kapa- 
een wie der Rotor selbst. Diese Koeffizienten sind überall im Rotor 
konstant ''). 

II. Die Wärmeverteilung an der Berandung des Rotors ist uns vollkommen bekannt. 
Insbesondere soll uns die Temperatur längs des Schnittes des Rotors mit der Welle (Teile 
der Schnitte Bı B, und 2, B, in der Abb. 2) bekannt sein. Die Wärmeerzeugenden 
Faktoren im Innern des Rotors sind uns ebenfalls quantitativ bekannt. 

III. Der Wärmezustand ist stationär. 


IV. Die Wärmequellen im Innern und auf der Berandung des Rotors sind rotations- 
symmetrisch angeordnet’). 


Bezeichnen wir mit x, y,2 die Kartesischen Koordinaten des Raumes, mit 7’ (x, y, 2) 
die Temperatur, mit P(x, y, 2)°), diejenige Wärmemenge, die im Innern des Rotors [an 


der Stelle (x, y, z2)] entwickelt wird, so besteht — was wir aus der Physik als bekannt 
übernehmen — die partielle Differentialgleichung 


wobei 7 (x, y,z) an der Berandung des Kreiszylinders (Rotors) die vorgeschriebenen 
Werte annehmen soll und die Konstante e gleich dem Produkt der spezifischen Masse 
mal spezifischer Wärme dividiert durch den Wärmeleitfähigkeitskoeffizienten des Ma- 
terials ist. 

Wir führen folgende Bezeichnungen ein: 4? sei der Laplacesche Operator; B— unser 
Kreiszylinder; 2A — seine Höhe; 20 — sein Durchmesser; € sei die Berandung von B (die 
Stirn- und Mantelflächen); dr — das Volumenelement von B; do — das Flächenelement von ('; 
n — die Normale (nach innen gerichtet) von C. Mit x, r =|\y:-+2°), g bezeichnen wir 
zylindrische Koordinaten des Raumes. Wir nehmen ferner an, daß die x Achse mit der 


Achse des Zylinders, die yz Ebene mit der Mittelebene COC (Abb. 2 S. 112) des 


Zylinders zusammenfällt. (Nach IV sind somit die sämtlichen auftretenden Funktionen 
von unabhängig.) 


Die Gl. (1) können wir auch in der Form 


schreiben. Um 7 aus (la) zu finden (vergl. Fußnote 5), spaltet man bekanntlich die 


Funktion 7 
T(®, Y, z) —u Y, z) +w (x, z) (2), 


"+ y—-n)'+(2— 5%” 


wobei 


ist, und w diejenige in B regulärharmonische Funktion ist, die an der Berandung ( die 
vorgeschriebenen Werte 


annimmt. 

Unsere ganze Aufgabe hat sich somit auf die Auffindung einer in BD regulärhar- 
monischen Funktion w reduziert, die auf € die durch (4) gegebenen Werte annimmt. 
In 2 werden wir zeigen, wie man (verhältnismäßig einfach) eine durch die Randwerte 
gegebene Potentialfunktion mit Hilfe von Polynomen approximieren kann, wobei mit 
wachsendem Grad des Polynomes, der Polynom in B gegen die gesuchte Funktion w strebt. 


2. Approximation der rotationssymmetrischen Potentialfunktionen in einen 
Kreiszylinder durch Polynome. In den Arbeiten »Ueber die Entwicklung der har- 
monischen Funktionen der Ebene und des Raumes nach Orthogonalfunktionen« und »Ueber 


') Die Lüftungskanäle im Rotor lassen sich bis zu einem gewissen Grade berücksichtigen, indem 
wir sie als negative Wärmequellen betrachten. 


2, Diese Vereinfachung ist nicht wesentlich. 


3) Wir können annehmen, daß Pix, y,z) eine stetig differenriierbare Funktion des Ortes "ist, was 
für die folgenden Ueberlegungen von Bedeutung ist. 
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die Bestimmung der elastischen Spannungen und Verschiebungen«!) wurde folgender Satz 
bewiesen: 


Zu jedem einfach zusammenhängenden ganz im Endlichen gelegenen konvexen Bereich 
B des dreidimensionalen Raumes läßt sich eine unendliche dreieckige Zahlenmatrix 


0 0 

0 


angeben, mit Hilfe deren man jede im abgeschlossenen Bereiche B reguläre Potentialfunktion ww 
in der Form 


v=m 
w = limes w„ = — limes 2 K,(x,y,z) do - (6) 


On 


darstellen kann. Dabei bedeuten K,(x,y,z) die Kugelfunktionen, 4A, eine Konstante und 


yP m @ =m 


p p=v v 

\ 

P=v p=v+1 p=v+2 

F =. m 

77 


Es läßt sich weiter zeigen, daß eine Funktion w, die in einen abgeschlossenen 
B beschränkt und in jedem ganz im Innern vom B gelegenen Teilbereiche von B regulär 


(8) 
| | Ik 
ist, sich ebenfalls in der Form (6) darstellen läßt. = WITT besitzt den Wert 
Kıl Ka! IK, K,-ı! 


I Kıt {Ki+ı Kat +1 Ks-ıt 


(1) — 
/ IK, Kat | IKı Ka! IKı Kst 
IKs-ı Ka} IKs-ı Ks-ı} | | ıK: Ki! IK Ka} IK: 
wo 
(5) (5, Iy ‘2 
ist. 


Der oben angegebene Satz läßt sich dazu verwenden, die durch die Randwerte 
gerebene Potentialfunktion w durch Polynome zu approximieren, indem für m endliche 


!) Mathematische Annalen 86 (1922) S. 238—271 und 98 (1927) S. 248-263. Ich möchte bei 
der Gelegenheit einige Druckfehler in der letzten Arbeit richtig stellen! Es soll S. 251, zweite Zeile von 


unten 3 (4) und nicht »R P, sy«; S.252. dritte Zeile von unten » Verschiebungs- 
| 


vektoren gr, r,r « und nicht »Versehiebungsvektoren Kr. r *; Seite 263, zehnte Zeile von unten: 
sals (41)2I)« und nicht »als (40)?"«, ferner »durch (38) und (41)« und nicht »durch (38) und (40)« 
heißen, 


nn | 


ir 


rt 
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Werte eingesetzt werden!)., Die Berechnung der Ausdrücke (8) ist im allgemeinen 
rechnerisch zeitraubend. 


Handelt es sich um Bereiche, die durch eine endliche Anzahl von Parametern charak- 


(s) 
9% 
terisiert sind, so liegt es nahe, die y =(— 1)"*‘ y ArEg ein für allemal als Funktionen 
95 9s+2 


dieser Parameter auszurechnen, um in der Zukunft unangenehme Determinäntenberech- 


(s) 
Ik 
nungen zu vermeiden. In der vorliegenden Arbeit werden wir die y” = (— 1)**‘ y 0) (+) 
9; +2 


für Kreiszylinder für den Fall berechnen, das w eine kreissymetrische Funktion ist. Die 
Folge der Kugelfunktionen reduziert sich dann auf die kreissymetrischen Kugelfunktionen ?), 
die in Zylinderkoordinaten lauten: 


> 9 3 3 
(10). 
21 5 105 
— — 
4 128 
IT (m) 
Das EN: w Ay do ist in diesem Falle: 
n 
m m) 
) 0) | )e az! (1i). 
)r )r 
Unsere ganze Aufgabe besteht darin, die in (8) durch die Determinanten ange- 
9% 
gebenen Werte der =, Fre für den Fall, daß B ein Kreiszylinder von der Höhe 
I, 9 +2 


2h und vom Durchmesser 2e ist, zu berechnen. gu” sind in A und 0 homogene Polynome, 


wie man leicht übersieht. Wir übergehen die rein formale Berechnung dieser Polynome 
und stellen sofort das Resultat auf. 


Es ist: —0 (n,m— ganz). IK, K, } ist für die übrige v, «dem Aus- 
2 


druck 4rho”*v | > a (.) | gleich. In der Zahlentafel I wurden die Werte 


Nn=( 


für kleinere zusammengestellt. 


!), Das auf diese Weise erhaltene wm stellt die Lösung der Minimumaufgabe dar, unter allen 


Polynomen mten Grades Gm (@,y,2) von der Form p=m 
p—=0 


dasjenige zu finden, das das folgende Integral zu einem Minimum macht. 


(w— Gm) Iw— Em)\? — Gm)\? 
(w — Em) — + )-+( ) dt 
Oz 


?) Die kreissymetrischen X, lassen sich bekanntlich in der Form 
Vy'+2? 
K,= P,, (cos #) E arctang — 
x 
schreiben, wo P(eos# die Legendreschen Polynome erster Art hedeuten. Die P,(eos 9) sind tabuliert. 
Vgl. das Verzeichnis berechneter Funktionentafeln, Berlin (1928), II, S. 21. 
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Zahlentafel I. 


(v, u) (v, (v, (v, v.) (v, (v, | 
A, A, A, A, Ay 
1 
IK; Kıt 
2 
Ka} 
1 
Kı, - 
3 
ı ha Kat 0 10 
N 4 
4 3 
\ 5 5 1 
\ 
16 4 2 
45 5 15 
Kı 
ET 128 16 14 
15 25 25 25 
Kı K 0 
128 14 18 
15 b 
Kr 0 | 
A 64 | 4 7 
9 15 | 9 33 f 
12 4 14 3 
75 45 | 45 60 15 es 
Ka Ko 
256 2 16 7 N 11 
35 35 21 1 
Kı\ 
j 128 16 8 2 
Ss 21 5 21 15 7 
Ka, 64 32 16 4 6 
at 175 35 35 15 35 35 
12 128 8 8 22 
175 1029 245 245 294 49 
IK; 0 — 
512 128 16 12 44 2. 
ı ha Ka; 29 48 g 
5 7 91 16 
128 16 16 11 
75 175 189 285 115 81 24 
usa 256 128 392 16 6 11 13 
ix x. 1225 175 245 161 3955 476 140 32 
BEE 4096 96 39 4 72 11 13 15 


. . (2x 
Die 9, dagegen besitzen die Werte g,, Ä (s, — ganz), 


pzn p= » » 


\ 


Die Werte der a‘). wurden für kleinere s und p in der Zahlentafel II angegeben. 


an 
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— 
— 


08098208 0 
181290°0 186618°'61 ‚819009°70 FE — — — — 02091700 — 02998000°0 — (0 
—D 
807 89897 809 691 col (2) 
sccHl 66 96 _ _ 960F 16 
8r07 96% 9887 L9L care (2) 
col 9 
r991 68 (9) 
F701 167 6 (9) 
N 
9] 9 (2) 
L 91 
L + r 91 15 
sr 81 :) 
(0 
| | Pr 13 
| 
TI 
| I (8)0 
| | | I 
| | | (6 
U®Z 


ew, 
) 
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Für manche Rechnungen ist es bequemer, wenn man die 9,” in der Form 
n=——— —) 
2 2 s2 v 2 2 h 5) n gerade) 
g, (4n) h » ) 


— (4) N’ Ip,2n ( ) 


darstellt. In der Zahlentafel III haben wir die Werte für q,),, für kleinere s,p (vor- 
wiegend in Form der Dezimalbrüche) angegeben. 


Um die g und or zu berechnen, benutzt man zweckmäßig die Formeln: 
=—ıK, + ıKs 4 Ks} iKs Ks | 


Damit ist unsere eigentliche Aufgabe gelöst. Es bleibt nur übrig das kons'ante 
Glied A„, das bei unserer bisherigen Rechnung noch unbestimmt geblieben war, zu be- 
rechnen. Am zweckmäßigsten ist es, A, derart zu bestimmen, daß 


w.)? do oder W— Um do oder do . (13) 


Zahlen 
| 
0,375 0,5 
0,5 
gu) 0,25 0,5 | 
0,25 0,66 0,5 
0,014648437 0,17578125 0.109375 0,16071428  0,08571428 
(9) | — 0,05859375 — 0,46875 0,25 0,2857142 | 
yo) 0,046875 0,375 0,2 | 
9) | 0,00343945312 0,140625 0,4546875 0,45535714 0,43809523 0,538 
gie) 0,0117185 — 0,609375 — 0,8125 0 + 0,9904761 |  0,20317460 
| 0,0078125 0,1875 0,666 0,6857142 0,12190476 
0,00005007 0,0036048889| — 0,041389465 — 0,15609783 — 0,12483629 — 0,02047821 
gi 0,00040055 0,02403259 0.23391723  0.,09521484 0.30338535  0,53503787 
0,0008010864 0.04165649 | — 0,37109375 — 0,9179687 0,15625 0,86363636 
gi 0,0004577636 0.021972656 0,20703125 0,6473214 0,59523 0,19047619 
| — 0,000005721875 — 0,00074387 0,01795955 0,1668222 — 0,6281225 — 1,0136123 
0,0000286090 0.00328066 0,07173153 0,58349623 1,88886621  2,5702121 
6 65536 32768 40960 5376 1680 
3 183 172 


2p+ 1 +) 
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Minimum wird. In unserem Beispiel (vgl. unter 4) haben wir A, mit Hilfe des letzten Inte- 
grals bestimmt. 


3. Fehlerabschätzung. Da eine harmonische Funktion das Maximum ihres Betrages 
stets am Rande annimmt, so erhalten wir als rohe Fehlerabschätzung die Tatsache, daß 
in jedem Punkte von B 


ist, wobei © das Maximum des Unterschiedes der Werte von w und w„ an der Be- 
randung ist. 


Unsere Methode bietet die Möglichkeit, eine bessere Majorante anzugeben. Es sei 


Der Wert %,(x, y, z) an der Berandung C ist bekannt, und wir wollen eine genauere 
Abschätzung für (x, z) im Innern von B finden. 
Bezeichnen wir mit yo: 


7 
Klar). 


die Orthogonalfunktionen, die zu e gehören (vgl. die erste ER 2 genannte Arbeit), 
so ist L 


(x, y, 2) 


s=m-+ On 
und Im (x, x,Y, 2) y, 2) 
s—m-+1 
tafel Ill. 


"0,08707482 


0,16968839 0,15800857 0,0410404 0,00219885 


| 
0,65728342 0,22510822 | 0,01539202 | 
| 
0,35670995  0,02770562 | 
| 
0,014512471 | 
—0,5351659 0,3534460 0,9089564  0,926533 0,3018747 | 0,0354281 0,00104974 
1,65922445  1,9601721 2,6089412 .1,01675040 | 0,13816930 | 0,00472364 
| 876544 16220160 
1,0604159 2,9474533 | 1,26924194 
4729725 | 2341213875 
937472 | 2048 
0,6174397 - '0,0032166170 
517440 24255 | 


| 
Ä 
N 
| 
| 
r- 
| 
te | 
e- 
| 
| 
1: 
| 
| 
| | 
| 
| 
| | 
| | 
| 
21 
87 
36 1 
19 | 
3 | 
1 
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Entsprechende Gleichungen lassen sich für die Ableitungen nach y bzw. z schreiben. 
Nach der Schwarzschen Ungleichung folgt: 


m- m—+1 


m 


Wir wollen nunmehr für die beiden unter der Quadratwurzel auftretenden Ausdrücke 
Majoranten einsetzen. Nach Formel (27) S. 263 der ersten am Beginn von 2 zitierten 
Arbeit ist 


„+1 4nA’(r,y, 2) 4n a? y,2) | 


wobei A (x,y,2) den Abstand des Punktes (x, y,:z) von der bedeutet. 
Um eine Majorante für den zweiten Ausdruck in (20) zu finden, ersetzen wir in 


0, durch einen anderen zweimal differentiierbaren Ausdruck ©, (&,y,2), der auf C dieselben 
Werte wie %,(x, y,:) annimmt. [©,„(x,%y,:2) braucht die Laplacesche Gleichung nicht 
zu befriedigen]!),. Da die harmonische Funktion die Eigenschaft besitzt, daß sie das 
Integral (22) zu einem Minimum macht, wenn zum Vergleich alle in B zweimal stetig 
differentiierbare Funktionen, die auf C mit 9,, (X, y, :) übereinstimmen herangezogen werden, 


so ist 


(' 


Ist w(x,r) auf der Berandung eine viermal differentierbare Funktion, so kann man 
in unserem Falle als ©, (x, y,:) den Ausdruck 


!) Die Heranziehung der Minimaleigenschaft des Diriehletschen Integrals zu Fehlerabschätzung 
wurde in letzter Zeit von E, Trefftz (über Theorie und Praxis «des Ritzscehen Verfahrens, Vortrag 
in der Gesellschaft für angew. Mathematik und Mechanik in Berlin am 13. Mai 1927, erscheint dem- 
nächst in den Mathem. Annalen unter dem Titel! »Konvergenzsätze und Fehlerschätzung beim Ritzschen 
Verfahren«) benutzt. Vergl. ferner die Arbeit von E. Trefftz, »Ein Gegenstück zum Ritzschen Ver- 
fahren«, Verhandlungen des 2, internationalen Kongresses für technische Mechanik, Zürich und Leipzig 
1927), S. 131 bis 137. 


Zahlen 
| 

r 0 0,1 0,2 0,3 0,4 0,5 0,6 | 0,7 

x 0,5 0,5 0,5 0,5 0,5 0,5 0,5 | 0,5 

Ka 0,25 0,245 0,24 0,205 0,17 0,125 0,07 0,005 
Kr 0,0625 0,055 0,0331 — 0,002 — 0,0479 . —0,1015 — 0,1589 — 0,2147 
0,015625 0,01108 —0,000875 — 0,01536 | — 0,02458 | —0,018378 0,01435 0,0880 
Ks 0,003906 000188 —0,00234 0,00403 0,0022 | 0,0187 0,0403 0,0488 
Wy 0,405 0,402 0,399 0,377 0,349 0,329 0,294 0,253 
u 0,794 0,780 0,755 0,682 | 0,594 0,487 0,363 0,262 
ws 1,1778 1,144 1,082 0,91 0,731 0,513 0,3028 0,131 
08 1,1004 1,347 1,228 0,978 0,731 0,469 0,212 0,051 


ke 
en 
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"9 (—h,r) (—h, r) 
(sin nr| dr +cos ur | dr) 


—2h(r—o) 
0 


2 


h —+h 
‘ . Im Ian ( 


n=( 
Im (h,r) 9m (h,r) 
+ (e-+h)(r—o) 2 (sin nr | nr | Ar dr) 
0 


benutzen. Ersetzt man in (20) den ersten Faktor durch (21) und den zweiten Faktor 
durch den in (23) auf der rechten Seite angegebenen Ausdruck, so erhält man eine 
Majorante, die man leicht angeben kann. 


4. Ein Beispiel, Im vorliegenden Abschnitt werden wir (unter Zugrundlegung 
der in 1 besprochenen vereinfachenden Annahmen) die Wärmeverbältnisse in einem 
kleinen Kurzschlußrotor betrachten, der keine Luftkanäle hat. Es sei A='/; 0 —=1. 
(Vergl. Abb. 2 und 3, S. 412.) Der Durchmesser der Welle sei ebenfalls 1. Wir nehmen 
ferner an, daß das magnetische Feld und die Eisen- und Hysteresis-Verluste so gering 
sind, daß sie keine Wärmequellen verursachen. Dagegen entstehe infolge der Kupfer- 
verluste längs der Mantelfläche B, cB, die Temperatur 50°; auf den Kreisscheiben vom 
Radius '/; um die Achse; auf den beiden Stirnflächen sei die Temperatur 20°, auf den 
übrig gebliebenen Kreisringen auf beiden Stirnflächen sei die Temperatur 50° Unsere 
Frage besteht darin, die Temperatur 7’ (x y,:z) im Innern des Zylinders (Rotors) aus diesen 
Angaben zu berechnen. Stellen wir 7’(x&,y,:) in der Form 


— 0 w@,y,2) - - (8) 


auf, so reduziert sich unsere Aufgabe darauf, eine in B reguläre harmonische Funktion 
zu berechnen, die auf den Kreisscheiben vom Radius '/s um die Achse auf den Stirnflächen 
den Wert 1 besitzt und auf den übrigen Teil der Berandung den Wert 0. 


Wir berechnen mit Hilfe der Zahlentafel II die zu unseren Zylinder (k=!/,0=|) 
(s) 
gehörigen g\” und daraus die Pe(—ı)t* '). Die Matrix (5) lautet: 


V 
(s) 


') Wir haben die y,’ für die ungerade 8 nicht ausgerechnet, da sie für unser Beispiel ohne 


Interesse sind. 


tafel IV 

0.8 0,9 | 1,0 1,0 1,0 1,0 1,0 1,0 
0,5 0,5 0,5 0,4 0,3 0,2 | 0,1 0 
0,011 — 0,155 — 0,25 — 0,34 — 0,41 — 0,46 — 0,49 |-05 

— 0,2643 — 0,298 — 0,3125 — 0,0794 0,1131 0,2566 0,345 | 0875 
0,2108 0,392 0,64068 0,3951 0,1338 | — 0,0995 — 0,25658 | — 0,3125 

— 0,00822 | — 0,139 — 0,487 — 0,510 | —0,811 | — 0,0345 0,188 0,273437 
0,2065 | 0,151 | 0,09% 0,038 — 0,005 — 0,035 — 0,055 — 0,06 
0,100 | — 0,049 — 0,183 — 0,097 -- 0,028 0,037 0,077 0,041 

— 0,013 — 0,055 — 0,019 | 0,062 0,049 — 0,047 — 0,023 — 0,052 


— 0,027 — 0,012 0,048 — 0,003 — 0,010 0,004 0,015 0,023 


| 

sch. | 

r | | 

EN. n | 

| 

| 

} 

in | 

| 

en | 

‚ht | 

as 

ig | 

| 

| 

| 
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0 0 0 0 0 0 0 
0 0,628 0 0 0 0 0 0 
0 0 0 0 0 
0 0,613 0 0,785 0 0 0 0 . 
0 0 0 0 0 
(26). 
0 0,7223 0 1,2572 0 1,0045 0 0 
0 0 0 0 
0 0,5682 0 1,5942 0 2,2451 0 1,2906 


Abb. 4b. 


Die Werte der Funktionen 
wy„n und w längs der Mantel- 
linie Bxc. 


Abb. 4a. 


Die Werte der Funktionen 
an und w längs des Ra- 
diusses ABı. 


[RA8237 4a u 
\n, 
B; 5, 
c 
A| 
Ir 
Abb. 2. Aufriß des Zylinders. Abb. 3.  Grundriß des Zylinders. 


Mit Hilfe von (7) berechnen wir 7” (für ein bestimmtes m). Die Werte von 


K: 

“dw lauten: 


| N 
m \ 
| 
— 
\ 
2 
| 


w, 
>h. 
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v 1 2 3 | 4 


0 1,57 0 0,196 
On 


| (27). 
0 — 0,258 0 — 0,0413 


OT (m) 


Auf Grund von (11) und (27) berechnen wir / J w 540. Abgesehen von der Be- 
n 


stimmung der Konstanten A,„, erhalten wir mit Hilfe von (6) die gesuchten Approximations- 
polynome. Unsere Rechnungen haben wir für m = 2, 4, 6 und 8 durchgeführt. Mit 
Hilfe der letzten in (13) angegebenen Nebenbedingung wurde A„ bestimmt. 


Die Werte von A, lauten: 


m | 0 | 2 4 6 8 (28) 
Am 0,125 0,255 0,414 | 0,503 0.522 
Wir erhalten mithin: 
= 0,125 
ws = 0,25 + 0,620 Ka 
ws = 0,414 + 1,30 Ka + 0,872 K, + IM). 


we = 0,501 2,11 2,28 Kı + 1,125 K; 
w = 0,522 + 2,470 Ka —+ 3,48 Kı + 2,546 Ks + 0,518 K; 


An der Berandung C nehmen XA;,„ (wie aus den Tabellen zu entnehmen ist) und 
Wg„ (wie wir berechnet haben) die Werte, die in der Zahlentafel IV (S. 410/411) ver- 


zeichnet sind, an. 


Die Werte der Approximationspolynome an der Berandung C' wurde in den Abb. 4a 
und 4b aufgetragen. Die Basis der Abb. 4a ist der Radius ABı, die Basis von 4b ist 
die halbe Leitlinie Bi c der Mantelfläche des Zylinders (vergl. Abb. 2 und 3). 823 


KLEINE 


Zur Ermittlung der theoretischen Leistung 
von Turbomaschinen. Da in vielen Fällen 
erhebliche Unterschiede zwischen der nach der 
Eulerschen Grundgleichung berechneten theo- 
retischen Leistung und der tatsächlichen Lei- 
stung sich ergeben, erschien es notwendig, in 
die Untersuchung der Bewegungsvorgänge in 
den Kreiselrädern tiefer einzudringen. In den 
veröffentlichten hydrodynamischen Arbei- 
ten, die diese Frage behandeln 
(z. B. den Untersuchungen von 
W. Spannhake!) und von E. 
Sörensen?), erfolgt die Bestim- 


mung des wirklichen Wertes 
der theoretischen Leistung auf 4 


Grund der bekannten Föttin- 
gerschen Gleichung. Im Fol- 


genden gestatte ich mir nachzu- 7 


) W. Spannhake, Anwendung 
der konformen Abbildung auf die 
Berechnung von Strömungen in Krei- 
selrädern. Hydraulische Probleme, 
S. 180 bis 200 (1926). £ 
9) E. Sörensen, Potentialströ- Pz 
mungen durch rotierende Kreisel- 
räder, diese Zeitschr. Bd. 7 (1927), 
Heft 2. 


weisen, daß die Ermittlung der theoretischen 
Leistung nicht auf Gıund der Föttinger schen 
Gleichung, sondern aus einer andern geschehen 
muß. | 

Aus dem Momentensatz, ergibt sich, wie be- 
kannt, daß für eine ebene Strömung in einem 
Zentrifugalpumpenrad (Abb. 1). 


(2) (1) 

W 


Abb. 1. 


| 

| 

| 

V 
N 
| 
| 
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wobei 1. — Leistung, M — Drehmoment, — Win- 
kelgeschwindigkeit, y — spezifisches Gewicht, 
9 — Erdbeschleunigung, u — Umfangsgeschwin- 
digkeit, c, Umfangskomponente der absoluten 
(reschwindigkeit c, dgq — Flüssigkeitsmenge, die 
durch den Umfangselement d’! fließt, bezeich- 
nen; die Indizes 1 wnd 2 entsprechen dem Ein- 
tritt- und Austrittumfang. 

Ist c„ die Normalkomponente der Absolutge- 
schwindigkeit, so ergibt sich 

und folglich 


g I 
(2) 1) 
oder 
2r 


L == (Cu Cn)2 dy ur, (Cu (1 a). 
q 
0 
Aus der Gleichung (1) oder (la) sieht man 
also, daß Z vom Produkt (eu €n), foglich von 
der gestrichenen Fläche der Abb. 1, abhängig 
ist. 


Die Gl. (1) kann in folgenden Formen darge- 
stellt werden: 
L= sın 2 a)g dis 
Eur 
(2) 
1 
(1) 
q 
(2) 
g 
(1) 
1 
(ug? — u?) — ug sin 2 3)a dia 
(2) 
+ u | w* sin 
(1) 


wobei @«— den veränderlichen Winkel (c, u), ıw,„ 
und ww, die Normal- und Umfangskomponen- 
ten der Relativgeschwindigkeit w, B — den ver- 
aänderlichen Winkel (w, u) bezeichnen. 

Aus der Gl. (1) ergibt sich 


L= ug du, —? u du 
(2) 


un (ug Cug — Cut) 
4 


wobei c„ı und cug unbekannte mittlere Werle 
von Cu; und cu» sind. Die Gl. (2) ist streng 
korrekt. Führt man die Annahme, daß c„ den 
Wert ctzß’” (3’-Winkel, unter dem die 
2ır 

Schaufel den Umfang schneidet) hat, so erhal- 
ten wir den von Euler gefundenen Wert von 

Da diese Annahme nicht immer der Wirk- 
lichkeit entspricht, so unterscheidet sich L von 
der nach Euler berechneten Leistung Z.. 
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Aus der Gl. (1) folgt auch, daß 


g g 
(2) (1) 


=" 
g 
wobei cy„ı und Cya unbekannte mittlere 
Werte von und c„9, und 79 — die Zir- 
kulationen um die Umfänge 1 und2 bezeichnen. 
Die Gl. (3) ist ebenfalls streng korrekt. 
gleich sei, so ergibt sich die Föttingersche 
Gleichung 


Führt man die Annahme, daß c„der Größe 


oder 
wobei z—Schaufelzahl, 7%, Zirkulalion um 


eine einzelne Schaufel. 
Aber es gibt keinen Grund anzunehmen, daß 


C„ der Größe En gleich sei; im Allgemeinen ist 
C„ veränderlich an den Umfängen 1 und 2, und 


es gibt keinen Grund dafür, daß die Größen 


[ 


denselben Wert haben. 

Demzufolge ist die Föttingersche Glei- 
chung wie die Eulersche nur eine ange- 
näherte, — um die Eulersche Gleichung 
zu bekommen wird ein hypolhetischer 
Wert für c„ eingeführt, und die Föltinger- 
sche Gleichung ergibt sich unter hypotheti- 
scher Annahme für den Wert von c,„ 

Danach muß man bei hydrodynami- 
schen Untersuchungen die Gl. (1) (oder eine 
von den Gleichungen 1a, 1b, Ic, 1d) und nicht 
die Gl. (4a) benutzen. In den oben erwähnten 
wichtigen Untersuchungen von W. Spann- 
hake und von E. Sörensen sind aber die 
Gl. (4) und (4a) verwendet für Berechnung der 
Leistung ZL:; demzufolge müßte man in die Er- 
sebnisse dieser Arbeiten entsprechende Berichti- 
sungen einbringen. 

Unsere Ableitungen gelten auch für die übri- 
gen Arten der Turbomaschinen (axiale Pumpen, 
Turbinen, Turbokompressoren usw.). 

B. Meisel. 897 

Ueber die Approximation konvexer ver- 
mittels linearer Funktionen. Will man eine 
vorgegebene Funktion f(x) im Intervall 
—a,b> durch eine lineare =pxr--gq er- 
setzen, so daß der durchschnittliche Fehler 


Charkow. 


— la) 


möglichst klein wird, so führt die Differential- 
rechnung auch bei einfachen f(x) auf unüber- 
sichtliche Bedingungen zur Bestimmung von 
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p und qg. Nun kann man im Falle einer kon- 
vexen (oder konkaven) Funktion f(x) diese 
Aufgabe graphisch und rechnerisch unmittelbar 
lösen. Es gilt nämlich folgender 

Satz. Ist f(x) in <a,b>> stetig und kon- 
vex (konkav) und sind x,. die Abszissen 
der Punkte die <a, b> in vier gleiche Teile 
zerlegen, so ist die durch /(x,. f(x,)) und 
(X, f(X3)) gelegte Gerade das Bild der gesuch- 
ten Funktion /(x). 

Beweis. 1. Es gibt eine Gerade die J/(l) 
minimisiert, denn F(p,qg) ist nicht negativ, 
stetig und für große 'p) + |g beliebig groß. 

2. Die minimisierende Gerade schneidet das 
Bild von y=f(x), denn sonst könnte man 
durch Parallelverschiebung der Geraden J(I) 
vermindern. 

3. Es muß zwei Schniltpunkte 5” und 5” 
geben, denn sonst würde eine Drehung um den 
einzigen J/(!) vermindern. 

4. Es seien /,. I, I, die Längen der drei 
Strecken, in welche die Strecke y--l(x), 
a—x—_b5b, durch 5 und 5” zerlegt wird; es 
sei also == Bei der Parallelverschiebung 
der Geraden um dq muß das Differential dJ 


verschwinden, wenn J minimal sein soll. Ein 
Blick auf die Zeichnung ergibt aber: 
wo 3 den Neigungswinkel der Geraden y = (x) 
gegen die x-Achse bedeutet. Daraus folgt 
+1; 
>». Bei der Drehung der Geraden um den 
Winkel d$ um 8° muß das Differential dJ 
ebenfalls verschwinden. Sind aber x,. %,, 7, die 
llächen der drei von den Kurvenbögen einer- 
seits und den Strecken /,. begrenzten li- 
suren, so ist 
dJ=dn +drne +dnz 
— 


— 
+2 — 1? +l2=0. 


Führt man (1) ein, so folgt 
I, l;, 


also 


also 1; und die Projektion der 
2 2 

Strecke 7, —-1,--1, auf die x-Achse teilt (a, b) 

im Verhältnis 1:2:1 w. z. b. w. 


l.wöw. II.Steinhaus. 


BUCHBESPRECHUNGEN 


(Die hier angezeigten Bücher sind dureh die VDI-Buchhandlung, Berlin NW 7, Ingenieurhaus, zu beziehen.) 


Dr.-Ing. OTTO MOHR, Wirklicher Geheimer 
Rat und Professor. Abhandlungen aus 
dem Gebiete der Technischen Mecha- 
nik. Dritte erweiterle Auflage zur Jahrhun- 
dertfeier der Technischen Hochschule Dresden. 
llerausgegeben von Dr.-Ing. K, Bever, o. Pro- 
fessor a. d. Techn. Hochschule Dresden, u. NH. 
Spangenberg, 0. Professor a. d. Techn. 
Hochschule München. Mit 528 Texlabbildungen 
Verl. Wilhelm Ernst & Sohn, Berlin 1928. 
XV --622S. Preis geh. 30M, geb. 33M. 

Eine schöne und würdige  Festgabe zur 
Jahrhundertlfeier der Technischen Hochschule 
Dresden bildet die Neuherausgabe der klassi- 
schen Abhandlungen Mohrs über Technische 
Mechanik. Man kennt den Inhalt der Ar- 
beiten, die zu ihrer Zeit großenteils grund- 
legende Lösungen für Probleme der Statik 
und der Festigkeilslehre gebracht haben. Die 
Mohrsche Methode für die Spannungsberech- 
nung an Fachwerken, für die Konstruktion der 
elastischen Linie des Balkens, für die geo- 
melrische Darstellung des Trägheils- und Span- 
nungstensors sowie vieles andere sind längst 
Gemeingut aller wissenschaftlich arbeitenden 
Ingenieure. In der vorliegenden Neuausgabe 
ist eine Abhandlung über das slalisch unbe- 
stimmte Fachwerk und ein kurzer Aufsalz 
unter dem Titel »Beitrag zur Berechnung der 
Rahmenträger« zum ersten Mal aufgenommen. 
Wenn diese Arbeiten auch bei dem heuligen 
Stand der technischen Mechanik kaum etwas 
wesentlich Neues bringen, so erfreut doch die 
bekannt klare Durchdringung der Aufgabe jeden 
sachverständigen Leser. Im ganzen genommen 


hat das vor mehr als zwei Jahrzehnten zum 

erstenmal erschienene Buch noch nichls von 

seiner Jugend und Aktualität verloren. 
Mises 


J. L. COOLIDGE, Professor der Mathematik 
an der Harvard University. Einführung 
in die Wahrscheinlichkeilsrech- 
nung. Deutsche Ausgabe von Dr. Friedrich 
M. Urban, Sekretär. der »Phönix-Wiener«, 
Vereinigte Versicherungs-A.-G. Sammlung ma- 
Ihematisch-physikalischer Lehrbücher. Heraus- 
segeben von E. Trefftz. Nr.24. Verlag von 
B. G. Teubner, Leipzig und Berlin 1927. 2128. 
Preis geb. 10M. 

Das vor kurzem erschienene englische Ori- 
vsinal dieses Buches ist in dieser Zeitschrift an- 
sezeigt worden). Das Werk gehört zu den 
verständigsten und  verständlichsten  Einfüh- 
rungen in die Elemente der Wahrscheinlich- 
keitsrechnung, die in den lelzten Jahren er- 
schienen sind. Die Ueberselzung, die der Ver- 
lag in handlicher Form herausgebracht hat, ist 
durchaus zu begrüßen. Gelegentliche Ueber- 
sehen des Ueberselzers, wie auf S. 11, werden 
hoffentlich bei einer Neuauflage in Ordnung 
sebracht werden. Mises. 9% 


Dr.”’Dr.-Ing. AUG. FÖPPL +, Professor an 
der Techn. Hochschule in München, Geh. 
Hofrat, und Dr. LUDWIG FÖPPL, o. Professor 
an der Techn. Hochschule in München. Dran« 
und Zwang. Eine höhere Festigkeilslchre Fur 


I) Diese Zeitschrift Bd. 6 (1926), S. 82. 


1- 
| 
_ 
| 
| 
), | 
1. 
) | 
j 
| 
/ | 
| 
| 


116 Buchbesprechungen 


Ingenieure. Mit 79 Abb. im Text. Zweiter Band, 
Zweite Auflage. R. Oldenbourg, München und 
Berlin 1928. VIIT --382S. Preis geh. 16M, geb. 
17 50M. 

Die erste Auflage dieses Werkes und der 
erste Band der zweilen sind hier schon ange- 
zeigt worden. Es mag genügen, zu der jelzt er- 
schienenen Fortsetzung zu sagen, daß der Ver- 
fasser Ludwig Föpp| viele neuere Fortschritte 
der Festigkeitslehre berücksichtigt und in den 
Zusammenhang des Buches glücklich hineinge- 
arbeitet hat. Dabei ist, um PJalz zu sparen, 
manches Ueberflüssige der ersten Auflage aus- 
seschieden worden. Mises. 95 


E. HELLINGER und 0. TOEPLITZ, Inte- 
sralgleichungen und Gleichungen 
mit unendlich vielen Unbekannten 
Sonderausgabe aus der Eneyklopädie der ma- 
Ihemalischen Wissenschaften. Mit einem Vor- 
wort von E. Hilb. Verlag B. @. Teubner, 
Leipzig und Berlin 1928. 1355 bis 1616 S. Preis 
geb. 16M. 

Der Sonderabdruck aus der Enzyklopädie der 
mathematischen Wissenschaften gibt eine außer- 
ordentlich sorgfältig bearbeitete Darstellung 
der Theorie der Integralgleichungen, wie sie 
ausgehend von dem lredholmschen Ge- 
danken hauptsächlich durch Hilbert und 
seine Schüler, durch Poincare, Erhard 
Schmidt und andere entwickelt worden ist. 
Die Verfasser waren bemüht, die grundlegenden 
Tatsachen in Zusammenhang mit den entschei- 
denden Beweisgedanken darzustellen, so daß 
ihre Arbeit mehr ist als sonst ein Enzyklo- 
pädie-Referat: Wer sich für irgend etwas inter- 
essiert, was mit Integralgleichungen cder mit 
der Idee der unendlich vielen Variablen zu- 
sammenhängt, auch für die prinzipiellen An- 
wendungen dieser Theorie, wird kaum an 
anderer Stelle so gute Auskunft finden wie in 
dem vorliegenden Buch. Mises. 925 


Dr. W. LIETZMANN, Göttingen. Ueber 
die Beurteilung der Leistungen in 
der Schule, Mathematisches, Psychologisches, 
Pädagogisches. Verlag von B. G. Teubner, Berlin 
1927. 1168. Preis kart. 6M. 

»Die Schrift kann als ein Versuch angesehen 
werden, der praktischen Pädagogik mit Psy- 
chologie und Mathematik zu Hilfe zu kommen.« 
Sie versucht, »die Leistungsbeobachtung psy- 
chologisch und das Psychologische in ihr ma- 
Ihematisch zu erfassen«. Das verwendete ma- 
Ihematische Rüstzeug ist kein schweres: Dia- 
gramme werden weitgehend benutzt, die Gauß- 
sche Kurve und der Korrelationsbegriff enl- 
wickelt, ohne daß viel mathematisch abgeleitet 
würde. Zur Erläuterung diene ein Beispiel: 
Man fragt sich, wie die verschiedenen Zen- 
surennummern 1, 2,3, 4,5 auf die Gesamtheil 
der Schüler verteilt sein mögen. Zunächst 
könnte man sagen, die normale Verteilung sei 
die, bei der mit jeder Nummer gleich viele 
Schüler bedacht werden; das entspricht aber 
offenbar in keiner Weise der Wirklichkeit. Dem 
entgegen steht das Postulat, die Verteilung der 
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Zensuren richte sich im Allgemeinen — und 
solle sich richten — nach der Gaußschen 


Kurve, wobei erklärt wird, auf welche Weise 
jeder Lehrer seine »persönliche Gleichung« — 
seine Glockenkurve — aufzustellen habe; das 
kommt natürlich auf die Bestimmung der zwei 
charakteristischen Konstanten hinaus. Die Frage 
des Zusammenhanges verschiedener »Rangord- 
nungen« führt zur Benutzung der Korrela- 
tionsformeln usf. Im ganzen stellt das Buch 
tatsächlich einen recht anregenden Versuch 
dar, der auch für den angewandten Mathema- 
tiker von Interesse ist. 

Berlin. 


Dr.-Ing. P. WERKMEISTER, ord. Professor 
a. d. Technischen Hochschule Dresden. Ein- 
führung in die Ausgleichungsrech- 
nung nach der Methode der kleinsten Qua- 
drate. Mit 34 Figuren im Text. Wittwers Tech- 
nische Hilfsbücher, Bd. 9. Verlag Konrad Wiltl- 
wer, Stuttgart 1928. VI--198$S. Preis geb. 
10,50 M. 

Das Büchlein enthält eine recht gelungene. 
elementar gehaltene Einführung in die Aus- 
sleichsrechnung, ohne auf die wahrscheinlich- 
keitstheoretischen Grundlagen einzugehen. Viele 
Beispiele vorwiegend aus dem Gebiete der 
Vermessungskunde, an denen die Fragestellung 
erläutert wird, erleichtern das Studium. — 
Nach einem Abschnitt, der die wichtigsten 
Termini erklärt, wird das Prinzip der Methode 
der kleinsten Quadrate dargelegt und schließlich 
die Aufgaben der Ausgleichung von direkten, 
von vermittelnden und von bedingten Beobach- 
tungen behandelt. 


Berlin. H.Pollaczek-Geiringer. 924 


H.Pollaczek-Geiringer. 921 


CARL FRIEDRICH GAUS$S, Bestimmung 
der Anziehung eines elliptischen 
Ringes. Nachlaß zur Theorie des Arithme- 
tisch-geometrischen Mittels und der Modulfunk- 
tion. Uebersetzt und herausgegeben von Dr. 
Harald Geppert, Privatdozent an der Uni- 
versität Gießen. Mit 4 Fig. Ostwalds Klassiker 
der exakten Wissenschaften Nr. 225. Akadem. 
Verlagsges. m.b. H., Leipzig 1927. 208S. Preis 
11,60 M. 

Der Verfasser hat es unternommen, in die- 
sem Bändchen zum erstenmal eine deutsche 
Ausgabe des Gaußschen Nachlasses zur Theorie 
der Modulfunktion zu verar ‘ten. 

Ueber diesen Gegenstand ha. Gauß selbst nur 
eine kleine Mitteilung über die Anziehung des 
elliptischen Ringes veröffentlicht, wobei er das 
arithmetisch-geomeltrische Mittel zur Berech- 
nung eines elliptischen Integrales verwendet. 
Alle anderen Arbeiten von Gauß, von denen 
einige tief in die Theorie der Modulfunktion 
eindringen, haben sich erst im Nachlaß von 
Gauß gefunden und wurden hauptsächlich 
durch die Bearbeitung von L. Schlesinger im 
10. Bande der Gesamtausgabe der Gaußschen 
Werke erst dem mathematischen Publikum zu- 
sänglich gemacht. 

Der Herausgeber hat nun auf Anregung und 
unter Beihilfe Schlesingers die Gaußschen Auf- 
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zeichnungen geordnet und hat es verstanden, 
indem er verbindenden Text sowie erläuternde 
Anmerkungen hinzufügte, einen eigenartigen, 
der Gaußschen Gedankenwelt entsprechenden 
l,ehrgang zur Einführung in die Theorie der 
Modulfunktion zu schaffen. Jedem, der sich 
eine konkrete Vorstellung über diese groß- 
artigen Schöpfungen von Gauß erwerben will, 
wird durch diese Art der Darstellung der Weg 
in wohl durchdachter Weise sehr geebnelt. 


Prag. P. Funk. 907 


Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 


Dr. ERICH MOSCH, Öberstudienrat am 
Mommsen - Gymnasium in Charlottenburg. 
Lehrbuch der Physik, Oberstufe. Viertes 
Heft: Von der Elektrizität und dem Magne- 
tismus. Mit 199 Abb. Freytag A.-G., Leipzig 
1927. X+2298S. Preis Hiblwd. 4,80 M. 


W.W.ROUSE BALL, Fellow and Tutor of 
Trinity Colledge (Cambridge). Histoire des 
Mathömatiques. Edition francaise revue et 
augmentee. Traduite sur la troisieme edition 
anglaise par L. Freund, Lieutenant de 
vaisseau. Tome Premier. Librairie Scienti- 
figque J. Hermann, Paris 1927. VII -+ 326 S. 
Preis 40 Fr. 


Dr.-Ing. OTTO KIRSCHMER, Vorstand des 
Forschungsinstituts für Wasserbau und Wasser- 
kraft e. V., München. Untersuchung der 
Ueberfallkoeffizienten und der Kolk- 
bildungen am Absturzbauwerk I im Sempt- 
flutkanal der »Mittleren Isar«, Heft 1 der 
Mitteilungen des Forschungsinstituts für 
Wasserbau und Wasserkraft e. V., München. 
Mit 44 Abb. R. Oldenbourg, München und 
Berlin 1928. 39 S., 1 Tafel. Preis geh. 4,50 M. 

Dr. med. JOHANNES LUDWIG SCHMITT, 
München. Kosmologie. Geheimnisse und 
Erkenntnisse. Dom-Verlag M. Seitz & Co., 
Augsburg. 34 S. Preis 1,50 M. 

Dr.-Ing. EWALD PERTZ, Die Bestim- 
mung der Baustoffdämpfung nach dem 
Verdrehungsausschwingverfahren. Mit 


42 Abb. und mit einem Vorwort von Prof. 
Dr.-Ing. ©. Föppl, Braunschweig. Sammlung 
Vieweg, H. 91. Verlag Friedr. Vieweg & Sohn 
A.-G., Braunschweig 1928. 62 S. Preis geh. 
3,60 M. 


Dr.-Ing. MAX MAYER, Beratender Ingenieur 
und Professor an der Staatl. Bau-Hochschule 
in Weimar. Nomographie des Bauinge- 
nieurs. Mit 47 Fig. Sammlung Göschen 
959. 111 S. Verlag von Walter de Gruyter 
& Co., Berlin und Leipzig 1927. Preis 1,50 M. 


Professor E. SCHULTZ, Oberlehrer an der 
Staatl. Maschinenbau- und Hüttenschule Duis- 
burg, und Gewerbeschulrat E. DIECKMANN, 
Oberstudiendirektor i. R. Mathematische 
und technische Tabellen für Bau- 
gewerkschulen und die Praxis. I. Zahlen- 
tafeln für den Hochbau. 239 — 250. Tausend 
der Schultzschen Tabellen. Vierzehnte, neu- 
bearbeitete Auflage von E. Dieckmann. 
Baedeker-Verlag, Essen 1928. XIV + 306 S. 


Dr. G. JUNGE, Professor, Bln.-Lichterfelde. 
Einführung in Wesen und Wert der 
Mathematik. Wissen und Wirken, 56. Bd. 
G. Braun, Karlsruhe 1928. 92 S. Preis 3M. 


GUSTAV RICHTER, Bozen. Die Kraft 
als fünfte Dimension. 1. Teil der Philo- 
sophie der Einmaligkeit. Wilhelm Braumüller 
Universitäts-Verlagsbuchh. Ges. m. b. H., Wien- 
Leipzig 1928. Preis 3,50 M. 


Dr.-Ing. RICHARD WAGNER, Die Be- 
stimmung der Dauerfestigkeit der 
knetbaren, veredelbarenLeichtmetall- 
legierungen. Berichte aus dem Institut für 
Mechanische Technologie und Materialkunde 
der Technischen Hochschule zu Berlin. Her- 
ausgegeben von Professor Dr.-Ing. P. Rieben- 
sahm. Erstes Heft. Julius Springer, Berlin 
1928. 64 S. Preis 6M. 


H. B. LÜBSEN, Ausführliches Lehr- 
buch der analytischen oder höheren 
GeometriezumSelbstunterricht. 17. Auf- 
lage, völlig neu bearbeitet von Prof. A. Do- 
nadt. Mit 129 Figuren im Text. Verlag 
Friedrich Brandstetter, Leipzig 1928. VII + 
291 S. Preis 6,50 M. 


NACHRICHTEN 


Gesellschaft für angewandte Mathematik 
und Mechanik. Die Hauptversammlung fand, 
wie angekündigt, in den Tagen vom 19. bis 
21. September in Hamburg statt. Gegenüber 
dem ursprünglichen Programm fielen die Vor- 
träge der Herren Burrau, Malkin und 
Schwerin wegen Verhinderung der Vortra- 
genden aus. Dafür trat Hr. L. Schiller- 
leipzig mit einem Vortrage über. das Problem 
der Wärmeübertragung ein. Hr. Kempf, der 
durch Krankheit gezwungen war, der Veranstal- 
tung fernzubleiben, wurde durch Hrn. Lerbs 
sowohl hinsichtlich des Vortrags wie in der 
Führung durch die Hamburger Versuchsanstalt 
ersetzt. 


Die Fachsilzungen der Gesellschaft waren von 
etwa 60 bis SO Mitgliedern und Gästen, darunter 


auch einer Reihe von Ausländern, besucht. Die 
folgende Uebersicht gibt die Titel der gehalte- 
nen Vorträge und die Namen der Diskussions- 
redner. Ein ausführlicher Bericht mit kurzer 
Inhaltsangabe sämtlicher Vorträge wird im 
nächsten Heft veröffentlicht. 


Erste Sitzung. 
Mittwoch, den 19. September, nachm.: 
Hydromechanik. 


H. Lerbs-Hamburg: Reibungsmessungen bei 
hohen Reynoldsschen Zahlen. — Aussprache: 
Reißner, v. Mises, Spannhake. 

A. Busemann-Göltingen: Zeichnerische Er- 
mittlung von ebenen Strömungen mit Ueber- 
schallgeschwindigkeit. — Aussprache: Reiß- 
ner, Konr. Müller. 
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J. Nikuradse-Göltingen: Turbulente Strö- 
mungen in keilförmisen Kanälen. Aus- 
sprache: Schiller. 

I. Hopf-Aachen: Ueber die Geschwindigkeils- 
verteilung an einer Plalte und an einem 
kohr .— Aussprache: Schiller, Reißner. 

Th. Troller-Aachen: Zur Wirbeltheorie der 
L.uftschrauben. Aussprache: Grammel, 
Reißner, Spannhake, Burgers. 

". Eisner-Berlin: Physikalisches zur Vereini- 
sung der Grenzschichttheorie mit der asym- 
p:olischen Widerstandstheorie. Aussprache: 
Zeilon, Reißner, Hamel, Burgers. 

M. Lagally-Dresden: Strömung im Außen- 
raum zweier Kreise. — Aussprache: Gram- 
mel, Busemann, Wilh. Müller, 
Bergmann. 


Zweite Silzung. 
Donnerslag, den 20. September, vorm.: 
Statistik und Versicherung. 


Dobbernack-Berlin: Die Einwirkungen 
von Krieg und Inflation auf die malhemati- 
schen Grundlagen der deuischen Sozialver- 
sicherung. -— Aussprache: Böhm, v. Mises, 
Riebesell. 

Die versicherungsmalthemalischen Grund- 
lagen der deutschen Arbeitslosenversicherung, 
Aussprache: Plaut, Riebesell, 

v. Mises. 

Basch-Wien: TFehlertensoren und Fehler- 
übertragung. Aussprache: v. Mises., 
Winkelmann. 

l. Baur-Berlin: Probleme der Mehrfachkorre- 
lation. Aussprache: v. Mises, Riebe- 
sell. Deuehler. 

Rh. Böhm - München: Einige Bemerkungen über 
die Theorie des Preises und über die damit 
zusammenhängenden ragen des  Grenz- 
nulzens in der theoretischen Nationalökono 
mie. Aussprache: Riebesell, Lorev 


Dritte Silzung. 
Donnerstag, den 20. September, nachm.: 
Mechanik. 

II. Ueber die Zugdiasonalen- 
felder in dünnen Blechen. Aussprache: 
HHamel, Trefftz, Reißner. 

II. Pollaezek-Geiringer-Berlin: Zur 
Praxis der Lösung linearer Gleichungen in der 


Stalik. Aussprache: Hamel,v. Mises, 
Trefftz. 

W. Müller-Hannover: Zur hydrodvnamischen 
Deutung der elliplischen Funktionen. Aus- 


sprache: Burgers. 

\M. Herzberger-Jena: Geometrische Optik 
und differentielle Liniengeomeirie. Aus 
sprache: v. Mises. 

St. Bergmann-Berlin: Leber die Berechnung 
des magnelischen Feldes in einem Transfor- 
malor. Aussprache: Emde, Stein, 
ineser. 

Sehiller-Leipzig: Untersuchungen zum 
Wärmeübergangsproblem. Aussprache: 
Busemann, Körner. 
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Vierte Sitzung. 

"reilag, den 21. September, nachm.: Gemeinsam 
mit der Deutschen Gesellschaft für technische 
Physik. 

(. Reutlinger - Darmsladt: Mechanische 
Schwingungsmesser hoher Empfindlichkeit. 
II. Backhaus-Berlin: Ueber Strahlungs- und 
Richtwirkungseigenschaften von Schallstrahlen. 
J. Trendelenburg-Berlin: Ueber Herztöne 
und lHerzgeräusche. 
Gauer-Göllingen: 

Siebkelten. 

K. Pohlhausen-Berlin: Fragen aus der Phy- 
sik der Ilochspannungsventilröhren. 

Die Geschäftssitzung der Gesellschaft am 
Freitag, den 21. September, fand infolge Ver- 
hinderung des Vorsitzenden Hrn. Prandtl 
unter Leilung des stellvertretenden Vorsitzenden, 
Hrn. Reißner, statt. Er gedachte zunächst 
der im abgelaufenen Jahre verstorbenen Mit- 
glieder, der Herren Krey-Berlin, Schachen- 
meier-München und A.Schoenflies-Frank- 
furt a. M. Nach dem Bericht des Geschäfts- 
führers ist die Mitgliederzahl im abgelaufenen 
Jahr von 322 auf 394 gestiegen. Neue Orls- 
gruppen wurden nicht errichtet. Der Kassen- 
abschluß zum 1. Oktober 1927 mit einem Bar- 
bestand von 513,88 M ist von den Kassenprüfern 
richtig befunden worden. Die Versammlung 
erteilt Entlastung. Für die satzungsgemäß aus- 
scheidenden Ausschußmitglieder Prof. Dr. E. 
Meißner und Geh. Reg.-Rat Prof. Dr. H. 
l.orenz werden von der Versammlung die 
Herren Prof. Ing. K. Körner-Prag und Prof. 
Dr. G. Hamel-Berlin gewählt Zu Kassen- 
prüfern werden Hr. Prof. Dr. E. Schwerin 
und Prof, Dr. E. Trefftz wiedergewählt. Der 
Mitgliedsbeitrag wird in der bisherigen llöhe 
beibehalten. Als Ori der nächsten Jahresver- 
sammlung wird nach längerer Aussprache mit 
allen Stimmen bei zwei Stimmenthaltungen 
Prag gewählt. Hr. Körner begrüßt im 
Namen der Prager Mitglieder «diesen Entschluß. 


Ueber Vierpole und 


Anwendung der Statistik in der Technik. 
Das Außeninslitut der Technischen Hochschule 
Berlin veranstaltet in den Monaten November 
und Dezember eine Vortragsreihe »Fabrikalions- 
kontrolle auf statistischer Grundlage«. Nähere 
Angaben hierzu beim Sekretariat der Hochschule. 


Persönliches. Hr. Dr.-Ing. C. Weber in 
Dortmund ist zum o. Prof. der Mechanik (Nach- 
folger von E. Trefftz) an der Technischen 
llochschule in Dresden ernannt worden. Hr. 
Trefftz hat den Lehrstuhl von M. Grübler 
übernommen. 

Ir. Dr. A. Berger, Direktor der Versiche- 
rungsgesellschaft Phönix in Wien hat sich an 
der Universität Wien für Versicherungsmathe- 
maltik habilitiert. 

Ir. a.o. Professor Dr. J. Lense an der 
Technischen IHlochschule München ist zum ord. 
Professor der angewandten Mathematik ernannt 
worden. 

Ir. Prvdzt. Dr. B.Schilling in Dresden ist 
zum nicht beamleten a. 0. Prof. ernannt worden. 


W. 
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ZUSCHRIFTEN AN DEN HERAUSGEBER 


Ueberschallgeschwindigkeit in zylindri- 
schen Rohren. Hr. E. A. Wedemeyer be- 
handelt in Heft 3, S. 238 die scheinbare Un- 
stimmigkeit der beiden experimentellen Tat- 
sachen, daß der Ausfluß eines komprimierten 
Gases durch ein zylindrisches Rohr im Rohr 
höchstens mit Schallgeschwindigkeit vor sich 
seht, während man ein in demselben Rohr be- 
findliches Geschoß mit einer höheren Geschwin- 
digkeit als Schallgeschwindigkeit hinaustreiben 


Jla 
Abb. 1. 


kann (Abb. 1). Zu der dort angegebenen Rech- 
nung ist folgendes zu bemerken: 


Zunächst ist das benutzte Integral der Schwin- 
gungsgleichung genau der Energiesalz. Außer- 
dem ist die Expansionsarbeit pro kg eines ab- 


Pg 
geschlossenen Gases A =/p -dv nicht gleich 
der Indikatorarbeit pro kg verarbeiteles Gas 


Pa 
L= fo . dp, denn beide unterscheiden sich 


p 
durch die Differenz der Ansaug- und Ausschub- 
arbeit: L=A--(p,v, — Beim Ausfluß 
eines Gases im Beharrungszustand bleibt die 
ganze Energie im Gase selbst, und zwar han- 
delt es sich in diesem Falle um die Arbeit L. 
Neben der Kontinuitätsgleichung für die durch- 
strömende Masse besteht bei verlustloser Strö- 
mung die Kontinuitätsgleichung der hindurch- 
strömenden Energie. Daraus ergibt. sich die 
2 

Gleichung L=”_, die besagt, daß die pro kg 
freigewordene Arbeit Z sofort in die kinetische 
Energie desselben kg übergeht. Für eine nicht- 
stationäre Strömung, wie sie das Austreiben 
eines Geschosses aus einem Rohre ist, kann 
diese Energiegleichung nicht genau so lauten, 
weil im allgemeinsten Falle auch noch andere 
Inergieabgaben möglich sind. 


l:s sind die beiden experimentellen Tatsachen 
also nicht ohne weiteres miteinander vergleich- 
bar. Da das Problem der Ausströmung eines 
(sases mit Schall- und Ueberschallgeschwindig- 
keil eine eigenartige Schwierigkeit aufweist, soll 
jedoch der einzig vergleichbare Spezialfall lang- 
samer Beschleunigung bei unendlich großem 
Behälter näher betrachtet werden. 


Für den stationären Ausfluß durch ein zylin- 

drisches Rohr berechnet man aus dem Energie- 
2 

u 

salz RZ und der Adiabatengleichung die 


Durchflußmenge pro qem Rohrquerschn 
Zu diesem Zweck ist in Abb. 2 die Dichte des 


Gases p=—- und die Geschwindig- 
v pi 

keit w aus dem Energiesatz dargestellt. Das 


Produkt von p und w ist die Durchflußmenge 
». An dieser Kurve, die auch bei Wede- 
meyer in Abb. 2 angegeben ist, interessiert 


| 
| 
| \ 
/ | \ 
/ | \ 
0 7% 
[RAsıZ2} 
2. 


nur die Tatsache, daß es eine maximale Durch- 
flußmenge Vmax bei dem sogenannten kritischen 
Druck p, und der Schallgeschwindigkeit für den 
kritischen Zustand w, gibt. Bei größeren oder 
kleineren Drücken geht daher eine geringere 
Menge durch das Rohr. Betrachtet man nun 
den Ausfluß durch das Rohr, so wird am Ein- 
lauf der kritische Druck px herrschen und es 
wird in das Rohr die Durchflußmenge max 
einströmen. Aus reinen Platzgründen muß nun 
im ganzen Rohr Schallgeschwindigkeit herr- 
schen, weil sonst die eingeströmte Menge nicht 
weilerfließen kann. Durch eine Einschnürung 
am Einlauf (Lavaldüse, Abb. 3) kann man es 


Abh. 3. 


dagegen erreichen, daß die Zuflußmenge soweit 
verringert wird, daß nunmehr genügend Platz 
dafür vorhanden ist, daß das Gas in dem zy- 
lindrischen Rohr auch mit höherer als Schall- 
seschwindigkeit strömt. Der Anschluß an den 
Außendruck erfolgt aber in beiden Fällen erst 
nach der Mündung des Rohres, wo für das Gas 
genügend Platz vorhanden ist, um mit jeder 
beliebigen Geschwindigkeit strömen zu können. 
Die eigenartige Tatsache, daß im Beharrungs- 
zusland das Gas in einem zylindrischen Rohr 
ohne Einlaufsverengung höchstens mit Schall- 
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geschwindigkeit strömt, erklärt sich also ledig- 
lich aus Platzgründen. 

Betrachtet man ein langsam beschleunigtes 
(eschoß in diesem Rohre, so werden die Platz- 
fragen wesentlich anders. Das Geschoß wird 
solange beschleunigt, als der Druck an der 
Rückseite größer ist als der Druck p.„' vor dem 
(reschosse. pa ist wegen des Aufstaues der vor 
dem Geschoß befindlichen Luft beim Austreiben 
größer als der Außendruck pa. Fliegt das 
(eschoß im betrachteten Augenblick mit der 
Schallgeschwindigkeit ıw,, so herrscht hinter 
dem Geschoß der kritische Druck p;,, der größer 
sein möge als der Druck pa und die Reibung. 
Das Geschoß wird daher noch weiter beschleu- 
nigt. Da das Geschoß aber dann mit größerer 
(Geschwindigkeit als Schallgeschwindigkeit fliegt, 
während die weiter hinten im Rohre nach- 
strömende Luft nur Schallgeschwindigkeit hat. 
wird hinter dem Geschoß Platz gewonnen, in 
den ein Teil der nachströmenden Luft ein- 
dringen kann und daher ohne weiteres Ueber- 
schallgeschwindigkeit annehmen kann. in Abb. 4 


WR 
[RAsız 2«] 
Abb. 4. 


sei dies noch näher erklärt. Bis zu einer Länge / 
habe das nachströmende Gas Schallgeschwindig- 
keit!). Es müssen dann alle Querschnitte A in 
diesem Bereich die maximale Durchflußmenge 
hindurchlassen. Für einen Querschnitt B hinter 
der Länge ! besteht dagegen dieser Zwang noch 
nicht und hier ist es daher möglich, eine Strö- 
nung mit Ueberschallgeschwindigkeit zu haben. 

Es ist also für den Fall langsamer Beschleu- 
nigung bestätigt, daß das Geschoß über die 
Schallgeschwindigkeit hinaus beschleunigt wer- 
den kann und daß das nachströmende Gas auch 
dann noch weiter einen Druck auf das Geschoß 
zu seiner weiteren Beschleunigung ausüben 
kann. Wenn auch die nichtstationären Vor- 
sänge beim Abschießen einer Kanone (wegen 
der großen rechnerischen Schwierigkeiten bei 

!), 7 ist bei idealer Strömung zeitlich konstant 
und stellt den Punkt dar, an dem das Geschoß 
die Schallgeschwindigkeit überschritt. Das Gas 
hinter dem Geschoß außerhalb ? stellt eine Ver- 
dünnungswelle in einem mit der Geschwindigkeit 
ws, mit bewegtem Koordinatensystem dar, die gegen 
die Strömung anläuft und deren Kopf dauernd in 
! bleibt. Vergl. Riemann-Weber, Die part. 
Differentialgleichungen, Bd. II, Braunschweig 1919, 
Ss. 503 oder Ackeret, Handbuch der Physik, 
Bd. VII, Berlin 1927, S. 322. 


der Behandlung nichtstationärer Vorgänge) 
nicht direkt als Beweis herangezogen werden 
können, so mag doch aus dem angegebene 
Beispiel hervorgehen, daß die maximale Ge- 
schoßgeschwindigkeit nicht die Schallgeschwin- 
digkeit ist. Es ist vielmehr so: Die Masse 
gleich hinter der Kopfwelle beim plötzlichen 
Oeffnen einer zylindrischen Rohrleitung kann 
Ueberschallgeschwindigkeit besitzen und diese 
wird ausgenutzt bei der Beschleunigung des Ge- 
schosses. Würde man dagegen Erbsen während 
der stationären Strömung in das Rohr hinein- 
werfen, so würden diese innerhalb des Rohres 
höchstens mit Schallgeschwindigkeit bewegt 
werden können. A.Busemann. 912 


* 


Erwiderung. Zweck meines Aufsatzes war, 
zu zeigen, daß der vielfach als Dogma aufge- 
faßle Satz, „Die Höchstgeschwindigkeit im zy- 
lindrischen Rohr ist die Schallgeschwindigkeit‘, 
nicht immer zutrifft, daß es vielmehr Fälle 
gibt, bei denen die Schallgeschwindigkeit ohne 
Anwendung einer Düse überschritten wird. In 
diesem Sinne freut es mich, daß der Herr Ein- 
sender bei anderer Darstellung und unter Her- 
anziehung eines anderen Beispiels zu demselben 
Ergebnis kommt, was den Wert der Arbeit 
nur erhöhen kann, ihr aber nicht widerspricht. 

Der Unterschied der Integrale Z und A ist 
allgemein bekannt, wie auch, daß für Pulver 
I,.=A mit guter Näherung gesetzt werden kann. 
Ich hatte geglaubt, deshalb nicht näher darauf 


eingehen zu sollen. Fernerist/p dv =/v dp=-I. 


p 

Ebenso glaubte ich, um die Arbeit nicht un- 
nülz lang zu gestalten, auf die Erwähnung der 
Sekundärwellen (Viellesche Wellen) verzich- 
ten zu sollen, da sie zwar stören, aber nicht 
das Wesen der betrachteten Erscheinung ver- 
ändern. 

Es sollte im Aufsatze nur gezeigt werden, 
daß nicht immer die Einfügung einer Laval- 
düse in eine zylindrische Leitung nötig ist, 
um Ueberschallgeschwindigkeit zu erzielen, daß 
vielmehr bei sinkendem Anfangsdruck im Be- 
hälter diese Erscheinung eintritt. Eine Erkennt- 
nis, die vielleicht für Auspufferscheinungen 
bei Maschinenzylindern von Wert sein kann 
(Stumpf, Die Gleichstromdampfmaschinen), 
da hier die Ausströmung wesentlich dem Aus- 
treiben eines masselosen Geschosses ähnelt. 

Zu der Unterdruckwelle muß gesagt werden, 
daß sie ihren Kopf nur solange in / haben 
wird, als das Geschoß noch beschleunigt wird. 
Würde es sich in dem Versuchsgefäß (Abb. 35, 
S.239) befinden, so würde es nach Ueberschrei- 
ten von dessen Mitte verzögert werden, wobei 
die Unterdruckwelle sich der Rückwand nähert 
mit zeitlich wachsender Geschwindigkeit. 


E.A.Wedemeyer. 912a 


(Redaktionsschluß 10. Oktober 1928.) 
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